Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	73,319	T→C	P401P (CCA→CCG)	thiP ←	thiamine ABC transporter membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	73,319	0	T	C	95.7%	212.6 / ‑1.4	70	P401P (CCA→CCG)	thiP	thiamine ABC transporter membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (51/15);  total (54/16)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

AACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTC  >  NC_000913/73074‑73564                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          aacaccGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAg                                                                                                                                                                                                                                                      <  2:469040/251‑1 (MQ=255) aaaaacGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTAGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCATATTGCGGCAg                                                                                                                                                                                                                                                      <  2:301157/248‑1 (MQ=255)    aacGCCACCACGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGcc                                                                                                                                                                                                                                                   <  1:600197/251‑1 (MQ=255)        ccaccaCGCCAAAATCCCCAATCGCCAGCCAGCAGGCACAAGCTAGCGCGTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCTCTTTTAAGCGCGACCAGCCTTCAATGACCCGCGACTGACACAACATACTGTAGCGGGCGGTGATGTCGCGCATCGGGTTTTCCAGTTCTTGAAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGAGTGCGGCAGGccaatg                                                                                                                                                                                                                                               <  2:648462/251‑5 (MQ=255)          accacGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGccaatggt                                                                                                                                                                                                                                             >  1:171170/1‑245 (MQ=255)              cGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtg‑ttg                                                                                                                                                                                                                                        >  1:547762/1‑251 (MQ=255)              cGCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtg‑ttg                                                                                                                                                                                                                                        >  1:628762/1‑251 (MQ=255)               gCCAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtg‑ttgt                                                                                                                                                                                                                                       >  1:239667/1‑251 (MQ=255)                 cAAAATCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑ttgttg                                                                                                                                                                                                                                     >  1:257383/1‑251 (MQ=255)                      tCACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGc                                                                                                                                                                                                                                  >  1:489217/1‑249 (MQ=255)                      ccACCAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAg                                                                                                                                                                                                                                <  2:702875/250‑1 (MQ=255)                         ccAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTaa                                                                                                                                                                                                                             <  2:595239/251‑1 (MQ=255)                          cAATCGACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTaaa                                                                                                                                                                                                                            <  1:574164/251‑1 (MQ=255)                               gACAGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCAGCAGGCCACTGGGG‑TTGTTGAGCAGTAAAaaga                                                                                                                                                                                                                        >  2:358489/1‑250 (MQ=255)                                  aGCACGCAGGCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑Gcc                                                                                                                                                                                                                   >  1:238228/1‑251 (MQ=255)                                         aggCAAAAGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGcc                                                                                                                                                                                                            >  1:497060/3‑251 (MQ=255)                                                   cAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑Tcc                                                                                                                                                                                                 >  1:582967/1‑251 (MQ=255)                                                             gcCAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGAAGATTGCGGCAGGCCAATGGGG‑GTTATGAGCAGTAAAAAGAA‑GCAGGTACCCAGCACAA‑TCCCCGGAATGGc                                                                                                                                                                                       >  2:479102/1‑251 (MQ=255)                                                               cAGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCg                                                                                                                                                                                      >  1:235300/1‑250 (MQ=255)                                                                aGTGGACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCCTGGCg                                                                                                                                                                                      >  1:444623/1‑249 (MQ=255)                                                                   ggACGTTTCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGAACCA‑TCCCCGGCATGGCGAggtt                                                                                                                                                                                 >  2:409091/1‑249 (MQ=255)                                                                          tCAGGGCGCGCAGCTCGACCACTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATg                                                                                                                                                                          >  1:648279/1‑251 (MQ=255)                                                                                         cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACCTACTGTAGCGGGCGGTGATCTCGCGCCTCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGAATTTGTGAAAATCACAATGACGTCATCAGATTGCGGCAGGCCAATGGTG‑TTGTTGCGAAGTAAAAAGCA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATAGCTAGGATCAACGTGCCgatcctttcaacc                                                                                                                                                           >  2:124544/3‑239 (MQ=255)                                                                                         cgaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGc                                                                                                                                                           <  1:532623/249‑1 (MQ=255)                                                                                          gaccacTTTTAAGCGCGACCAGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTCCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACCATGCCGTCAGCAGATTGCGGCAGGCCGATGGTG‑TTGATGAGCAGGAAAAAGAA‑GCCGGTAGGCAGCACAA‑TCCCCGGCATGGCGAGGATACACGTGCCGCTCATCTccggcg                                                                                                                                                          >  2:30381/2‑247 (MQ=255)                                                                                                      gcgcgACCAGCCTTCAATGCCCAGCGACTGACACAACATACGGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCa                                                                                                                                              <  1:409091/247‑1 (MQ=255)                                                                                                       cgcgACCAGCCTTTAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATAGCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAg                                                                                                                                             <  1:487243/248‑1 (MQ=255)                                                                                                              aGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTc                                                                                                                                      >  1:701057/1‑251 (MQ=255)                                                                                                              aGCCTTCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTc                                                                                                                                      >  1:15299/1‑251 (MQ=255)                                                                                                                  ttCAATGCCCAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACAAGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGcc                                                                                                                                  >  2:675116/1‑251 (MQ=255)                                                                                                                     aatgcccAGCGACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCcgcc                                                                                                                               >  1:180463/5‑250 (MQ=255)                                                                                                                          ccAGCGACTGACACAACATACTGTAGCGGGCGGTGATCTCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGCAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAACA‑GCCGGCAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCACATCTCAAGCGCCTGACCCGCCCGCATTTTCTGcccaccccc                                                                                                                           >  2:655913/1‑243 (MQ=255)                                                                                                                               gACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAg                                                                                                                                                                                                                                   >  2:522224/1‑143 (MQ=255)                                                                                                                               gACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAg                                                                                                                                                                                                                                   <  1:522224/143‑1 (MQ=255)                                                                                                                               gACTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTc                                                                                                                     >  1:628589/1‑251 (MQ=255)                                                                                                                                aCTGACACAACATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCCCCCGAAGTTcg                                                                                                                    >  2:322371/1‑251 (MQ=255)                                                                                                                                        aaCATACTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCCGCAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCCTCTCCAGCCCCTGACCCCACAGCCTTTTCTGCCGCGccctaaattctcgactgcc                                                                                                            >  2:142582/1‑235 (MQ=255)                                                                                                                                           atacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCCTCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGCCTGCTcc                                                                                                          >  1:130781/4‑250 (MQ=255)                                                                                                                                            tacTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATAAACCTGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCCCGCCCGAAGTTCGCGACTGCTccct                                                                                                        >  2:71379/3‑249 (MQ=255)                                                                                                                                                gTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGGTCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCACTGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGACCAA‑TCCCCGGCATGGCGAGGTTCAACATGCCGCTCATCTCAAGCGCCTGACCCGCCAGCATTTTCAGCCGCGCCCGAAGTTCGCGACGGCTACAtcac                                                                                                     >  2:27616/1‑247 (MQ=255)                                                                                                                                                     gggcggTGATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGACCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGAAGATGGCGGCAGGCCAATGGTG‑TTGTTGAGCATAAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGcatc                                                                                               <  1:251620/251‑1 (MQ=255)                                                                                                                                                            gATATCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAAGGATAGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAAGGGTG‑TTGGTGAGCAGTAAAAAGAA‑GCCGGTCGCCAGGACAA‑TCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCCGACCAGCCAGCATTTTCTACAGAGCCAGGAGTACGCGACGGCTCCATAGCAGAATCATGGTc                                                                                         >  2:368314/1‑250 (MQ=255)                                                                                                                                                               atCGCGCATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGc                                                                                      >  1:108139/1‑250 (MQ=255)                                                                                                                                                                  gcgcATCGGGTTTTCCAGTACTTTCAGCGCCTAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCCCAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGTAAAACGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGCGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCCTTTTCTGCCGCGCCCGCACTTCGCGCCCGCTCCCTAGCCGCCTCCTGGTCAGccac                                                                                   >  2:77243/1‑247 (MQ=255)                                                                                                                                                                   cgcATCGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTAc                                                                                >  1:168520/1‑251 (MQ=255)                                                                                                                                                                        cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGtg‑tt                                                                                                                                                                                                                                         <  1:180645/96‑1 (MQ=255)                                                                                                                                                                        cGGGTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGAGGCAGGCCAATGGtg‑tt                                                                                                                                                                                                                                         >  2:180645/1‑96 (MQ=255)                                                                                                                                                                        cGGGGTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGGG‑TTGTTGAGCAGAAAAAAGAA‑GCCGGGAGCCAGCACAA‑TCCCCGGCATGGCGGGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CCACGcac                                                                            >  1:124368/1‑250 (MQ=255)                                                                                                                                                                           gTTTTCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCCGCA‑CTACGCACAAt                                                                         >  1:492711/1‑250 (MQ=255)                                                                                                                                                                             tttCCAGTACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCCGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGATCCCTAGCAGCATCATGGTCAGCA‑CTACGCCCAATac                                                                       >  1:710645/1‑250 (MQ=255)                                                                                                                                                                             ttcccagttcttacggcgcATAAGGGATCGC‑ATATAACGTATGGGTAAAAATCACAGTGCCTGCAGCAGATGGCGGCAGCCCAATGGTG‑TTGTTGAGCATAAAAAAGAG‑GCCGGTAGCCAG‑ACAACTCCCCGGCATGGCGAGTATCAACATGCCTCCCCTCTCCAGCGCCTGCCCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATac                                                                       <  2:560012/235‑1 (MQ=255)                                                                                                                                                                                 cagtaCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCAGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCAGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTATGCCGCGCCCGAAGTTCGGGACTGCTCCATAGCAGCATCATGGACAGCA‑CTACGCACAATACACCtt                                                                  >  2:498160/5‑250 (MQ=255)                                                                                                                                                                                    taCTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTgccc                                                               >  1:170675/2‑250 (MQ=255)                                                                                                                                                                                          cAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGCG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAgcc                                                         >  2:110044/1‑250 (MQ=255)                                                                                                                                                                                             cgaaTAAGGGATCCCCAT‑TACTCCATTGGTAAAAATCACAATGCCGCCAGCAGCTTGCGGC‑GGGCAATGGGGGTGGTTGAACAG‑AAAAAGAAAGCGGGTGCCCACCACAA‑CCCCCGGCATGGCGAGACCCAACATGCCGCTCATCTCCAGCGCTGGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CGACGCACAATACACCTGCCGCCAGCGCaa                                                      <  2:31981/248‑1 (MQ=255)                                                                                                                                                                                              gcATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTTGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCCCAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAAt                                                     >  1:315422/1‑251 (MQ=255)                                                                                                                                                                                              gcATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAAt                                                     >  1:150998/1‑251 (MQ=255)                                                                                                                                                                                                aTAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATAc                                                   >  1:534879/1‑251 (MQ=255)                                                                                                                                                                                                      gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                             >  1:378357/1‑251 (MQ=255)                                                                                                                                                                                                      gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                             >  1:361226/1‑251 (MQ=255)                                                                                                                                                                                                      gATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACg                                             >  1:172415/1‑251 (MQ=255)                                                                                                                                                                                                                  aCGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAgc                                  >  1:91573/1‑251 (MQ=255)                                                                                                                                                                                                                    gCATTGGTGAAAATAACTATGCCGTCAGCAGATTGCGGCAGGCCCATGGTG‑TTGTTGCGCAGTAAACAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCACTCATCTCCACCGCCTGACCCGCCAGCAATTTATGCCGCGCCCGCAGCTCGCGCATGCTCCCTAACCGCATCATGGTCACCAACAACAC‑CAACACCCCTGACGCCAGCGCCATCCCCAACGAGGCCcacacccc                                >  2:618329/1‑247 (MQ=255)                                                                                                                                                                                                                      aTTGGTGAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCt                              >  1:78648/1‑251 (MQ=255)                                                                                                                                                                                                                         ggTGAAAATCACAATGCCGTCAGCAGATTGCGGGAGGCCAATGGTG‑TTGTAGAGCAGAAAAAAGAA‑GCGGGTAGCCAGCACAA‑CCCCCGGCAGGGCGGGGATCAACAGGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGcc                           <  1:289090/251‑1 (MQ=255)                                                                                                                                                                                                                          gagAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCCTGGGCAGCA‑CTACGCACAATACACCTCACGCCCGCGCCACACGCCACGCGGTCCACAgcccccgcc                           >  2:28280/3‑243 (MQ=255)                                                                                                                                                                                                                            gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCCTAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCACCGAGGGCCACAGCGCCTGCcacc                        >  2:18106/1‑250 (MQ=255)                                                                                                                                                                                                                            gAAAATCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTg                                                                  >  1:292238/1‑209 (MQ=255)                                                                                                                                                                                                                             aaaaTCACAATGCCGTCAGCAGATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGGCCAAAGCGCCTGCcac                         >  1:39739/1‑249 (MQ=255)                                                                                                                                                                                                                                                  gATTGCGGCAGGCCAATGGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCGCCTGACCCGCCAGCATTTTCTGCCGCGCCCGAAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCACTACGCAACGAGGTCCACAGCGCCTGCCACAGCAACGGATGTGCCCGCACTTc  >  2:555823/1‑251 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                          AACAACGCCACCACGCCAAAATCACCAATCGACAGCACGCATGCAAAGGCCAGCGCCTGCGCCAGTGGACGTTTCAGGGCGCGCAGCTCCACCACTTTTAAGCGTGACCAGCCTTCAATCCCCAGCGACTGACATAACATGCTGTAGCGGGCGGTGATATCGCGCATCGGGTTTTCCAGCACTTTCAGCGCATAAGGGATCGCCAT‑TAACGCATTGGTGAAAATCACAATGCCGTCAGCAGATTGTGGCAGGCCGATAGTG‑TTGTTGAGCAGTAAAAAGAA‑GCCGGTAGCCAGCACAA‑TCCCCGGCATGGCGAGGATCAACATGCCGCTCATCTCCAGCACCTGACCCGCCAGCATTTTCTGCCGCGCCCGCAGTTCGCGACTGCTCCATAGCAGCATCATGGTCAGCA‑CTACGCACAATACACCTGCCGCCAGCGCAATACGCAACGAGGTCCACAGCGCCTGCCACAGCACCGGTTGTGCCAGCACTTC  >  NC_000913/73074‑73564

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