Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	72,517	T→G	R137R (CGA→CGC)	thiQ ←	thiamine ABC transporter ATP binding subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	72,517	0	T	G	100.0%	188.3 / NA	54	R137R (CGA→CGC)	thiQ	thiamine ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (34/20);  total (34/20)

GTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGA  >  NC_000913/72222‑72794                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                      gTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcgct                                                                                                                                                                                                                                                                                   >  1:599132/1‑300 (MQ=255) gTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcgct                                                                                                                                                                                                                                                                                   >  2:578712/1‑300 (MQ=255) gTGCATACTAACCCGTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACgcgct                                                                                                                                                                                                                                                                                   >  2:483914/1‑300 (MQ=255)               gTAATCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAGTGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGAccgc                                                                                                                                                                                                                                                                       <  2:405792/298‑1 (MQ=255)                  aTCCCCAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAANGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGaa                                                                                                                                                                                                                                                                  >  1:145352/1‑300 (MQ=255)                    ccccAATAGTGCTGAANCACTCTCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCCCACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTAGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACTAGACAACNCGCTAACGCCACGCGCTGNCGCTGACCGCCGGaaa                                                                                                                                                                                                                                                                 >  1:581160/1‑299 (MQ=255)                    ccccAATAGTGCTGAANCACTCTCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACNCGCTAACGCCACGCGCTGNCGCTGACCGCCGGaaa                                                                                                                                                                                                                                                                 >  1:581092/1‑299 (MQ=255)                    ccccAATAGTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCAACAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGaaa                                                                                                                                                                                                                                                                 >  1:240087/1‑299 (MQ=255)                            gTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGNATCTTCCACACTGTGTGAAACCATCAATAACGNCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTAAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGAAAACGAGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGcca                                                                                                                                                                                                                                                         >  2:513315/1‑298 (MQ=255)                            gTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCAAACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATTGGCTGTTCGCGCACCAGCCAACGCGCTAACGCAACGTGCTGTCGCTGACCGCCGGAAAGCTCGcca                                                                                                                                                                                                                                                         >  2:60019/1‑298 (MQ=255)                            gTGCTGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCTTAACCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAACAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGcca                                                                                                                                                                                                                                                         >  2:59394/1‑298 (MQ=255)                                           ccTTACCGCTCAACAACTCATTGGTCTTACCCTGCCGGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGAAGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGAGCTAACGCCACGCGCTGTAGCTGACCGCAGGAAAGCTCGCCAGGTAACCGTGCCATTa                                                                                                                                                                                                                                         >  2:401310/1‑300 (MQ=255)                                             ttACCGCTCAACAACTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTaaa                                                                                                                                                                                                                                       >  1:423995/1‑300 (MQ=255)                                                      aacaacTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAAt                                                                                                                                                                                                                               <  1:267130/299‑1 (MQ=255)                                                           cTCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCTCGCCATTAAATTATCAATCCCCAt                                                                                                                                                                                                                         >  1:230412/1‑300 (MQ=255)                                                            tCATTGGTCTTACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATc                                                                                                                                                                                                                        >  1:394166/1‑300 (MQ=255)                                                                  gtcttACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTggc                                                                                                                                                                                                                    >  1:194273/5‑298 (MQ=255)                                                                  gtcttACCCTGCCAGGCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTggc                                                                                                                                                                                                                    >  1:194311/5‑298 (MQ=255)                                                                                 gCGATGCGCCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCAACAGCGTCAGCATCTACTGACGTNACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCAtt                                                                                                                                                                                                   >  2:504134/1‑300 (MQ=255)                                                                                      gcgcCCGTCTGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGTCGGGCGATAGCGTGCATTTTctc                                                                                                                                                                                              <  2:394166/300‑3 (MQ=255)                                                                                        gcCCGTCGGCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCct                                                                                                                                                                                            >  1:384812/1‑300 (MQ=255)                                                                                            gTCGGCGACCACCACCGAGCGCGTGGCAATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACTTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAAGAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTCCTCctgc                                                                                                                                                                                          <  2:470104/298‑1 (MQ=255)                                                                                               ggCGACCACCACCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGtgc                                                                                                                                                                                     <  2:563704/300‑3 (MQ=255)                                                                                                  gcccacccacGAGCGCGTGGCGATCCGCGCGGCATCTTCCACCCTGTGTGACACCATCAATAACGTCCTTTTTTGCTGCTGGCCGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGACGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTCACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGAAAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCATGCTGTACcgc                                                                                                                                                                                  >  2:427590/10‑300 (MQ=255)                                                                                                  gaccaccacCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACTTCATTTTTTGCTGCTTGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCcg                                                                                                                                                                                   <  1:343824/292‑1 (MQ=255)                                                                                                  gaccaccacCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCcg                                                                                                                                                                                   <  1:246201/292‑1 (MQ=255)                                                                                                  gaccaccacCGAGCGCGTGGCGATCCGCGCGGCATCTTCCACACCGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCctg                                                                                                                                                                                           >  2:427584/8‑291 (MQ=255)                                                                                                                                   aTCTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGATAACGCCAAGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCAATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCAGGATTTAGCCCCAGca                                                                                                                                                  >  2:35968/1‑298 (MQ=255)                                                                                                                                     cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGa                                                                                                                                               >  1:169332/1‑300 (MQ=255)                                                                                                                                     cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTACTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGa                                                                                                                                               >  2:399109/1‑300 (MQ=255)                                                                                                                                     cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGACAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGa                                                                                                                                               >  2:163959/1‑300 (MQ=255)                                                                                                                                     cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGa                                                                                                                                               >  1:366539/1‑300 (MQ=255)                                                                                                                                     cTTCCACACTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCAATCTGACGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCCGccc                                                                                                                                                 >  2:240445/1‑298 (MQ=255)                                                                                                                                          acacTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGTCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTc                                                                                                                                          <  1:536562/297‑1 (MQ=255)                                                                                                                                          acacTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTc                                                                                                                                          <  1:536572/297‑1 (MQ=255)                                                                                                                                             cTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCtgc                                                                                                                                       >  1:341154/1‑299 (MQ=255)                                                                                                                                             cTGTGTGACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCTCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCtgc                                                                                                                                       >  1:434651/1‑299 (MQ=255)                                                                                                                                               gtgtgACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACTCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCtgcgc                                                                                                                                     >  1:513717/5‑297 (MQ=255)                                                                                                                                               gtgtgACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCtgcgc                                                                                                                                     >  1:484993/5‑297 (MQ=255)                                                                                                                                               gtgtgACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCAGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCCGTTTCAATCCCGGATTTAGCCCCAGCCCGATGttatgcgc                                                                                                                                     >  2:59492/5‑294 (MQ=255)                                                                                                                                               gtgtgACACCATCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTACTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCtgcgc                                                                                                                                     >  1:686702/5‑297 (MQ=255)                                                                                                                                                         aTCAATAACGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAgg                                                                                                                           <  2:666543/300‑1 (MQ=255)                                                                                                                                                                acGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGNCAGGTAACCGCGCCATTAAATTATCAATCCCCNTCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGa                                                                                                                    <  2:513717/299‑1 (MQ=255)                                                                                                                                                                     aTTTTTTGCTGCTTGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCNGGTAACCGCGCCATTAAATTATCAATCCCCATNTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAAc                                                                                                                  <  2:517335/297‑1 (MQ=255)                                                                                                                                                                              tgctgGCAGCTCGTGCTCACCAGCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTcc                                                                                                      <  2:555664/300‑1 (MQ=255)                                                                                                                                                                                                    gCGTCAGCATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGAc                                                                                <  1:163959/300‑1 (MQ=255)                                                                                                                                                                                                               tcCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGAcggc                                                                     >  2:433021/1‑300 (MQ=255)                                                                                                                                                                                                                      gTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcg                                                              <  1:629659/300‑5 (MQ=255)                                                                                                                                                                                                                           gCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtggt                                                          >  1:477098/1‑291 (MQ=255)                                                                                                                                                                                                                           gCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCACACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACggcggcgtcgtggt                                                          >  1:477128/1‑291 (MQ=255)                                                                                                                                                                                                                                                        gggCAATAAAATCGACTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAAGTATCAATCCCCATCTGGCGGGCGAGAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGGCAGCGAAc                             <  2:430546/297‑1 (MQ=255)                                                                                                                                                                                                                                                        gAGCAATAAAATCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAAc                             <  2:430551/299‑1 (MQ=255)                                                                                                                                                                                                                                                               aaaaTCGGCTGTTCGCGCACCAGACAACGCTCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCgctggct                     <  1:272352/300‑1 (MQ=255)                                                                                                                                                                                                                                                               aaaaTCGGCTGTTCGCGCACCAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCgctggct                     <  1:271928/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                  ccAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAg                                      <  1:28072/264‑1 (MQ=255)                                                                                                                                                                                                                                                                                  ccAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTCCTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAACGCCGTCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCAGCGa  >  2:305190/1‑300 (MQ=255)                                                                                                                                                                                                                                                                                  ccAGACAACGCGCTAACGCCACGCGCTGTCGCTGACCGCCGGAAAGCTCGCCAGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCTACTGCTGTGCCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGCGCGACCGTCAGGTGACTGAACAGGTTGTTTTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCGTCGTGGTGTGATCAATGCCGTCGATAg                                      >  2:28072/1‑264 (MQ=255)                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                      GTGCATACTAACCCGTAATCCCCAATAGTGCCGAAGCACTCGCCTTACCGCTCAACAACTCATTGGTCATACCCTGCCAGGCGATGCGCCCGTCGGCGACTACTACCGAGCGCGTGGCGATCCGCGCCGCATCTTCCACGCTGTGCGACACCATCAATAGCGTCATTTTTTGCTGCTGGCAGCTCGTGCTCACCAGCGTCAACATCTCCTGACGTAACGCCGGATCGAGCGCAGAGAACGGTTCATCGAGCAATAAAATCGGCTGTTCGCGTACCAGACAACGCGCTAACGCCACTCGCTGTCGCTGACCGCCGGAAAGCTCGCCCGGTAACCGCGCCATTAAATTATCAATCCCCATCTGGCGGGCGATAGCGTGCATTTTCCCCTGCTGTACCGCGTTCAGTTTCAATCCCGGATTTAGCCCCAGCCCGATGTTCTGTGCGACCGTCAGGTGGCTGAACAGGTTGTTCTCCTGAAACAGCATCGACACCGGACGGCGTGACGGCGGCATAGTTGTGTGATCTACGCCATCGATAGTCAGCGAACCGCTGGCTGGCGTCAGAAAACCGGCGA  >  NC_000913/72222‑72794

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