Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,581,601	G→A	A656T (GCA→ACA)	tktB →	transketolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,581,601	0	G	A	100.0%	63.5 / NA	20	A656T (GCA→ACA)	tktB	transketolase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/11);  total (9/11)

GGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACA  >  NC_000913/2581465‑2581717                                                                                                                                         |                                                                                                                     ggAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACaa                                                                                                                   <  2:400148/140‑1 (MQ=255)  gAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACaaa                                                                                                                  <  2:154418/140‑1 (MQ=255)  gAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACaaa                                                                                                                  <  1:470633/140‑1 (MQ=255)  aaaGCAGGTATTGCCGACTACTCGTACAAGTTTGTCGGTTTGAAAGGGTCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACaaa                                                                                                                  <  2:178167/139‑1 (MQ=255)        ggTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCATTGGCTTAACCGCCCAGAATAATGTGACAAAAGCGCa                                                                                                            >  1:247888/1‑140 (MQ=39)           aTTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATaa                                                                                                         <  2:384681/140‑1 (MQ=255)             tGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCGCGTTCTTTGGCTGTACAGCCGAGAAGATTGTGACAAAAGCGCATAAgg                                                                                                       >  2:188479/1‑140 (MQ=255)                         gTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGa                                                                                           >  1:612511/1‑140 (MQ=255)                         gTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGa                                                                                           >  2:604984/1‑140 (MQ=255)                         gTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGGGACAAAGGCGCATGAGGTGCTGGGAGTGa                                                                                           >  1:303469/1‑140 (MQ=255)                          tACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGaa                                                                                          <  1:457038/140‑1 (MQ=255)                                    ttGGTCTGAAAGGGGCAATTGTAGGGATGACGGGTTACGGGGAATCTGCTACGGCGGAAAAGCTGTTCACGATATTTGGCTTTACAGCCGCGAATACTGTACAAAAAGCGCATAAGGTGCTGGGCGCGAACGGTACCTgc                                                                                >  2:15911/1‑139 (MQ=255)                                                      ttGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtgattgccggatgct                                                              <  2:521055/140‑1 (MQ=255)                                                                  cGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtgattgccggatgctgattgccggatg                                                  <  1:433743/140‑1 (MQ=255)                                                                  cGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGtgattgccggatgctgattgccggatg                                                  <  2:612511/140‑1 (MQ=255)                                                                       tACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACg                                             >  1:577439/1‑140 (MQ=255)                                                                           gggAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGa                                         <  2:577439/140‑1 (MQ=255)                                                                                tCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGAGTGCCGGAGGCTAATGGCCGGAAGCGACGCTGACGCGt                                    >  1:198483/1‑140 (MQ=255)                                                                                                       cccGTTCTTTGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTa             >  2:113965/1‑140 (MQ=255)                                                                                                               ttGGCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCt     >  2:496965/1‑140 (MQ=255)                                                                                                                  gCTTTACCGCCGAGAATATTGTGACAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACa  <  1:496965/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                     GGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAAGCGCATAAGGTGCTGGGAGTGAAAGGTGCCTGATGGTGATTGCCGGATGCTGATTGCCGGATGCGACGCTGACGCGTCTTATCCGGCCTACAGCGTCTTATCCGGCCTACA  >  NC_000913/2581465‑2581717

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