Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,644,343	A→T	5.8%	F280L (TTT→TTA)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,644,343	0	A	T	5.8%	221.9 / 2.7	82	F280L (TTT→TTA)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base A (32/45);  new base T (3/2);  total (35/47)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.78e-01

GCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACATCTCTTCACGCA  >  NC_000913/3644195‑3644491                                                                                                                                                     |                                                                                                                                                     gCAACTCATCGACGCCAAATTCGTCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTa                                                                                                                                                      >  2:608168/1‑149 (MQ=255)  cAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTaa                                                                                                                                                     >  2:313477/1‑149 (MQ=255)     cTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCTTTTTGCCAGATAGGTTAAGAAATACAGCACCTGCTGCGGGTTTTCTGCCATTTTTGTGGCAAGGGATTTAAAgg                                                                                                                                                  >  1:794996/1‑149 (MQ=255)      tCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGc                                                                                                                                                 <  2:434942/149‑1 (MQ=255)      tCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGGGTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGc                                                                                                                                                 <  1:346287/149‑1 (MQ=255)      acATCGCCGCCAAATTCGGCTTTGGCAAACGCACGCAATTGGGCCAGCTCTTTTTCGCCTTGGGGACGCGCGCGTTTTGCCAGAGCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGc                                                                                                                                                 <  2:710070/148‑1 (MQ=255)           gACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTTAAGGCGTAGt                                                                                                                                            >  2:438571/1‑149 (MQ=255)           gACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTTAAGGCGTAGt                                                                                                                                            >  2:433326/1‑149 (MQ=255)            aCGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGCCGCGCGCGTTTTGCCAGAGCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGtt                                                                                                                                           >  1:251018/1‑149 (MQ=255)              gCCAAATTCGGCTTTGTCAAAGTCACTCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTGTGCCAGATCGTTTAAGAAATCCAGCACCTTCTGCGGGTTTTCTGCCATTTTTGTGTCAAGGGATTTAAAGGCGTAGtttt                                                                                                                                         <  2:82934/149‑1 (MQ=255)              gCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTTAAGGCGTAGtttt                                                                                                                                         <  1:438571/149‑1 (MQ=255)              gCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTTAAGGCGTAGtttt                                                                                                                                         <  1:433326/149‑1 (MQ=255)               ccAAATTCGGCTTTGGCAAAGGCCCGCAATTGCGCCAGCTGTTTTTCGCCTGGTGGACGCGCGGGTTTTGCCAGATCGGTTAAGAATTCCAGCTCCTGCTGCGGGTTTTCTGCGATTTTAGTGGCTCGGGCTTTAAAGGCGTAGttttg                                                                                                                                        >  2:667864/1‑148 (MQ=255)               ccAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTc                                                                                                                                        <  1:1033751/149‑1 (MQ=255)                 aaaTTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCGTTCTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATGTAAAGGCGTAGTTTTCaa                                                                                                                                      >  2:1118845/1‑149 (MQ=255)                  aaTTCGGCTTTGGCAAATGCATGCAATTGCGGCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCtaa                                                                                                                                     <  2:624464/149‑1 (MQ=255)                     tCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGcc                                                                                                                                  <  1:1075023/149‑1 (MQ=255)                     tCGGCTTTGACAACGGCACGCAATTGCGCCAGCTCTTTTTAGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGcc                                                                                                                                  <  2:130824/149‑1 (MQ=255)                       ggCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGGCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAACTCCAGCCCCAGCTGCGGGTTTTCTGCGATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGcccc                                                                                                                                >  1:1014043/1‑148 (MQ=255)                       cgCTTTGCCAAAGGCACGCAATTGCGCCAGCGCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGACATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCca                                                                                                                                <  2:328027/148‑1 (MQ=255)                         cTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCcagc                                                                                                                              >  1:356520/1‑149 (MQ=255)                                aaaGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTgcg                                                                                                                       >  1:470529/1‑149 (MQ=255)                                aaaGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGATgcg                                                                                                                       >  2:189048/1‑149 (MQ=255)                                    gCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACCCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGAGATTTTATTGGCAAGGGATTTTAAGGCGTAGTTTTCAAAGCCCAGCAGTTACGCCAg                                                                                                                   >  1:894702/1‑149 (MQ=255)                                        gCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCg                                                                                                               >  1:622157/1‑149 (MQ=255)                                         cAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGt                                                                                                              <  1:1127624/149‑1 (MQ=255)                                           aTTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGa                                                                                                            <  1:474463/149‑1 (MQ=255)                                              gcgcCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGc                                                                                                         <  1:415410/149‑1 (MQ=255)                                              gcgcCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTAAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCAATTTAGAGGCAAGCGTTTTAAAGGCGTAGTTATCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGc                                                                                                         >  2:929042/1‑149 (MQ=255)                                              gcgcCAGCGCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTAAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGAAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGc                                                                                                         <  1:342908/149‑1 (MQ=255)                                                     ctctTTTTCGCCTTGTGGACGCGCGCGGTGTGCCAGATCGGTTACGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGa                                                                                                  <  1:1056960/149‑1 (MQ=255)                                                      tctTTTTCGCCTTGTGGACGAGCTCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTCCCATTTTCGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAg                                                                                                 <  2:231469/149‑1 (MQ=255)                                                           ttCGGCTTGTGGACGCGCGCGTTTTGCCAGATAGGTTAAGAAATCAAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGATTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGAtct                                                                                            <  2:1037260/149‑1 (MQ=255)                                                           ttCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGAtct                                                                                            >  2:311966/1‑149 (MQ=255)                                                            tCGCCTTGTGGACGCGCGCGATTTGCCAGTTCGGTTAAGAAAACCAGCACCTGCTGAGTGTTTACTGCCAGTTTAGTGGCAAGTGATTTAAAGGAGTAGTTTTCAAAGCACACCAGAAGCGCAAGTTCTTGACGCAGCGCGAGGAtctc                                                                                           >  2:514264/1‑149 (MQ=255)                                                              gCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCtt                                                                                         <  1:560723/149‑1 (MQ=255)                                                               ccTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTc                                                                                        <  1:494287/149‑1 (MQ=255)                                                               ccTTGTGGACGCGCGCGTTTTGCCAGAGCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTc                                                                                        <  1:175965/149‑1 (MQ=255)                                                                     ggACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCAc                                                                                  >  2:365667/1‑149 (MQ=255)                                                                      gACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCAcc                                                                                 <  2:353989/149‑1 (MQ=255)                                                                         gcgcgcgTTTTGCCATATCGGTTAACTAATCCAGCAACTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTg                                                                              <  2:1076514/149‑1 (MQ=255)                                                                         gcgcgcgTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTg                                                                              <  1:526113/149‑1 (MQ=255)                                                                         gcgcgcgTTTTGCCAGAGCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTg                                                                              <  1:431595/149‑1 (MQ=255)                                                                           gcgcgTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCGGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCt                                                                            >  1:525351/1‑149 (MQ=255)                                                                           gcgcgGTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCt                                                                            <  1:799219/149‑1 (MQ=255)                                                                                   tGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATccc                                                                    <  1:576917/149‑1 (MQ=255)                                                                                    gCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGATCGTGACGCAGCGCGAGTATCTCTTCCATCACCTAGCTGTTATCCCa                                                                   >  2:1133221/1‑149 (MQ=255)                                                                                    gCCAGATCGGTTAAGAAATCCAGCACCTCCTGCGGGTTTTCTGCCATTGTAGTGGCAAGGGAGTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGATCGTGACGCAGCGCGAGTATGTCTTCCATCACCTTTCTGTTATCCCa                                                                   >  2:249581/1‑149 (MQ=255)                                                                                     ccAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGGTTTCTGCCATGTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCCCAGCCCAGCAGTTGCGCCAGGTCGTGCCGCAGCGCGAGGATCTCTTCCATCCCCTTGCTGTTATCCCAt                                                                  >  2:297199/1‑149 (MQ=255)                                                                                         aTCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTAc                                                              >  1:831158/1‑149 (MQ=255)                                                                                               tAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGt                                                        >  2:1029128/1‑149 (MQ=255)                                                                                                aaGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGtt                                                       >  1:192697/1‑149 (MQ=255)                                                                                                aaGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGtt                                                       >  1:125272/1‑149 (MQ=255)                                                                                                aaGAAATCCAGCACCTGCTGCGGGATTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTATCAAAGCCCAGCAGTTGCGCCAGATCGTGACGCAGCGCGATGATCTCTTCCATCACCTTGCTGGTATCCCATTTACAGGCGtt                                                       >  1:234717/1‑149 (MQ=255)                                                                                                 tgataTCCAGCAGCTGCTGCGGGTTTTCTGCCATTTTAGTGGCATGGGATTTAACGGCGTAGTTTTCATCGCCCAGCAGTTTCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCATGCTGTTATCCCTTTTACCGGCGttt                                                      <  2:1014043/145‑1 (MQ=255)                                                                                                    aaTCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTggg                                                   <  1:313477/149‑1 (MQ=255)                                                                                                     aTCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGc                                                  <  2:356520/149‑1 (MQ=255)                                                                                                       ccAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCt                                                <  1:248326/149‑1 (MQ=255)                                                                                                             ccTGCTGCGGGTATTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCCTCGTTTTCAACGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTACTGTGAACCCCTTTACAGACGTTTGGGCCTTGAtag                                          >  1:773193/1‑147 (MQ=255)                                                                                                              ctgctgCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGACAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCgg                                         >  1:60020/1‑149 (MQ=255)                                                                                                               tgctgcGGGTTTTCTGCCATTTTTGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGAGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGa                                        <  2:1023119/149‑1 (MQ=255)                                                                                                               tgctgcGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGa                                        <  1:578790/149‑1 (MQ=255)                                                                                                               tgctgcGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGa                                        <  1:608168/149‑1 (MQ=255)                                                                                                                   gcGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGgcg                                    >  2:210504/1‑149 (MQ=255)                                                                                                                   gcGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGgcg                                    >  1:60463/1‑149 (MQ=255)                                                                                                                      ggTTTTCTGCCATTTTAGTGGCTAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCgg                                 <  2:978400/149‑1 (MQ=255)                                                                                                                       gTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCggg                                >  1:582903/1‑149 (MQ=255)                                                                                                                       gTTTTCTGCCATTTTAGTGGCAAGGGAGTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTGCGTGACGCAGCGCGAGGATCGCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCggg                                <  2:62451/149‑1 (MQ=255)                                                                                                                          ttCTGCCATTTTCGTGGCATGGGAGGTAAAGGCGTAGTGTTCACAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCGCATTTACCGGCGTTTTGGCCTTGATCGGAGGCGCTGGTGc                             <  1:500647/149‑1 (MQ=255)                                                                                                                             tGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCACAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGt                          <  1:1137073/149‑1 (MQ=255)                                                                                                                               ccATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAACGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGATGCTGTaa                        >  1:924596/1‑149 (MQ=255)                                                                                                                                   tttAGTGGCAAGGGCTTTAAAGGCGTAGTTTTCAAAGCCCAGCCGTTGCGCCAGTTCGTGACGCAGCGCGAGTATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAgcgc                    <  2:794996/149‑1 (MQ=255)                                                                                                                                   tttAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAgcgc                    >  2:506325/1‑149 (MQ=255)                                                                                                                                     tAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGTGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGa                  >  1:947163/1‑149 (MQ=255)                                                                                                                                     tAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGa                  <  1:925636/149‑1 (MQ=255)                                                                                                                                          gCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACAt             <  2:1179375/149‑1 (MQ=255)                                                                                                                                          ccAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGCTCGTGACGCAGCGCGAGGATCTCATCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACAt             <  2:1174856/148‑1 (MQ=255)                                                                                                                                           cAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACAtc            >  2:1153849/1‑149 (MQ=255)                                                                                                                                            aaGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACAtct           <  2:582903/149‑1 (MQ=255)                                                                                                                                              gggATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACATctct         <  2:125272/149‑1 (MQ=255)                                                                                                                                                gATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGAGCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTGTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACATCTCTTc       <  2:218007/149‑1 (MQ=255)                                                                                                                                                 aTTTAAAGGCGTAGTTATCAAAGCCCAGCAGATCCGCCAGGTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACATCTCTTCa      <  1:823118/149‑1 (MQ=255)                                                                                                                                                   ttAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTAACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACATCTCTTCACg    >  1:1042443/1‑149 (MQ=255)                                                                                                                                                    tAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACATCTCTTCACGc   <  2:622157/149‑1 (MQ=255)                                                                                                                                                     aaaGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGGGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACCTCTCTTCACGCa  >  2:33308/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCAACTCATCGACGCCAAATTCGGCTTTGGCAAAGGCACGCAATTGCGCCAGCTCTTTTTCGCCTTGTGGACGCGCGCGTTTTGCCAGATCGGTTAAGAAATCCAGCACCTGCTGCGGGTTTTCTGCCATTTTAGTGGCAAGGGATTTAAAGGCGTAGTTTTCAAAGCCCAGCAGTTGCGCCAGTTCGTGACGCAGCGCGAGGATCTCTTCCATCACCTTGCTGTTATCCCATTTACCGGCGTTTGGGCCTTGATCGGAGGCGCGGGTGCTGTAAGCGCGATACATCTCTTCACGCA  >  NC_000913/3644195‑3644491

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑