Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,579,943	G→A	G103D (GGC→GAC)	tktB →	transketolase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,943	0	G	A	100.0%	64.5 / NA	20	G103D (GGC→GAC)	tktB	transketolase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (12/8);  total (12/8)

ATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG  >  NC_000913/2579709‑2580171                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     aTCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGtt                                                                                                                                                                                                                      >  1:199088/1‑251 (MQ=255)                  ggCATGGCTGATATTGCCGAAGTGCTGTGGACCGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTACGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATCCGCCAGGCGTTGAAACCACCACCGGccc                                                                                                                                                                                                     >  2:19318/1‑251 (MQ=255)                              atTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAgg                                                                                                                                                                                         >  1:144829/1‑251 (MQ=255)                                    gAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGc                                                                                                                                                                                   >  1:16243/1‑251 (MQ=255)                                     aaGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCg                                                                                                                                                                                  <  1:16891/251‑1 (MQ=255)                                         gCTGTGGAACGATTTTCTTAAACATAACCCTCCCGCCCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCTGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAAc                                                                                                                                                                               <  1:258965/250‑1 (MQ=255)                                                                                  ccTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAg                                                                                                                                     <  1:113222/251‑1 (MQ=255)                                                                                                                   aCGGTAACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCAc                                                                                                    >  1:195068/1‑251 (MQ=255)                                                                                                                         aCGCGTCGATGCTGATCTACAGTTTGCTACATCTGACAGGTAACGACCTGCCGATGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTCTACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCCGACCATGAGATCGTCGATCACTTCAcc                                                                                              >  2:112647/1‑251 (MQ=255)                                                                                                                           gcgTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTa                                                                                            <  1:40868/251‑1 (MQ=255)                                                                                                                             gTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGCCCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTAtg                                                                                          <  2:141853/251‑1 (MQ=255)                                                                                                                                  tgctgcTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGttt                                                                                     >  1:168892/1‑251 (MQ=255)                                                                                                                                     tgctCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATg                                                                                  >  1:267545/1‑251 (MQ=255)                                                                                                                                                                          gCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGa                                               >  1:37216/1‑249 (MQ=255)                                                                                                                                                                            cGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGCAGGTATTTCCCACGAAGTc                                           >  2:197183/1‑251 (MQ=255)                                                                                                                                                                                         gaagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcag                              <  2:194573/251‑1 (MQ=255)                                                                                                                                                                                          aagaactTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTggcagg                             >  2:113974/1‑251 (MQ=255)                                                                                                                                                                                                      cAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACt                 >  1:145057/1‑251 (MQ=255)                                                                                                                                                                                                                    agccccAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATTGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTg   <  1:240905/249‑1 (MQ=255)                                                                                                                                                                                                                     aCCCCAGGCCACCCGGAGATTGACTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTgg  <  1:157508/251‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                     ATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCACCACCGGCCCGCTTGGACAAGGTTTGGCGAACGCCGTCGGGCTGGCGATAGCAGAGCGTACACTGGCGGCGCAGTTTAACCAGCCAGACCATGAGATCGTCGATCACTTCACCTATGTGTTTATGGGCGACGGCTGCCTGATGGAAGGTATTTCCCACGAAGTCTGTTCGCTGGCAGGCACGCTGGGACTGGGCAAGCTGATTGG  >  NC_000913/2579709‑2580171

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