Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,579,823	A→G	N63S (AAC→AGC)	tktB →	transketolase 2, thiamine triphosphate‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,579,823	0	A	G	100.0%	157.9 / NA	45	N63S (AAC→AGC)	tktB	transketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (29/16);  total (29/16)

ATGCGGTACAAAAAGCCAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCA  >  NC_000913/2579682‑2579966                                                                                                                                              |                                                                                                                                                aTGCGGTACAAAAAGCCAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCAcgc                                                                                                                                        >  2:755068/1‑151 (MQ=255) aTGCGGTACAAAAAGCCAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGCATGATCGCGACCGCTTTATTCTTTCCAGCggcc                                                                                                                                            >  2:315743/1‑145 (MQ=255)       tACAAAAAGCCAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGAt                                                                                                                                  >  2:356159/1‑151 (MQ=255)        aCAAAAAGCCAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGAt                                                                                                                                  <  1:682456/150‑1 (MQ=255)           aaaaGCCAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGAtgctg                                                                                                                              >  2:698541/1‑151 (MQ=255)                 cAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTAc                                                                                                                        >  2:658511/1‑151 (MQ=255)                        ggTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGc                                                                                                                 <  2:421890/151‑1 (MQ=255)                               ccGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCt                                                                                                          >  2:833074/1‑151 (MQ=255)                               ccGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCt                                                                                                          <  1:781476/151‑1 (MQ=255)                                         gATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTAc                                                                                                >  1:226135/1‑151 (MQ=255)                                              gCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCt                                                                                           <  1:699954/151‑1 (MQ=255)                                                   gCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGc                                                                                      <  2:174271/151‑1 (MQ=255)                                                    cTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGcc                                                                                        <  1:504897/148‑1 (MQ=255)                                                    cTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCt                                                                                     >  1:658645/1‑151 (MQ=255)                                                       atatTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTgg                                                                                   >  1:471114/1‑150 (MQ=255)                                                            gCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGgaagaac                                                                             >  2:473678/1‑151 (MQ=255)                                                             ccGAAGTGCTGTGGAACGATTTTCTTAAACATAGCCCTACCGTCCCAACCTGGTATGATCGCCCCTGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACATTTTTGTACATTTGCCCGGTTACGACTTGCCGCTGgaagatc                                                                             <  1:548701/150‑3 (MQ=255)                                                                      tgtgGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACtt                                                                   >  2:667798/1‑151 (MQ=255)                                                                        tgGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTcc                                                                 >  2:284598/1‑151 (MQ=255)                                                                         ggAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCg                                                                <  2:479232/151‑1 (MQ=255)                                                                             cGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAg                                                            <  1:857998/151‑1 (MQ=255)                                                                              gATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAg                                                            >  2:307772/1‑150 (MQ=255)                                                                              gATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGt                                                           >  2:338351/1‑151 (MQ=255)                                                                               aTTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGtt                                                          >  2:765131/1‑151 (MQ=255)                                                                                  ttCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGc                                                        <  1:770640/150‑1 (MQ=255)                                                                                       aaaCATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCg                                                   >  2:27067/1‑150 (MQ=255)                                                                                        aaCATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGaa                                                 <  2:512712/151‑1 (MQ=255)                                                                                         aCATAACCCTACCGACCCAACCTGGGATGATCGCGACCGCTTTATTCTTTCTAGCGGTCACGAGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGCTACGACCTGGCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGaa                                                 >  1:381965/1‑150 (MQ=255)                                                                                          cttagcgcgaCCGACCCAACCTGGCATCATTCCAACCGATGAATTCCATCCGGCTGTCACGCGTCGATGCTGCTCTCCAGATTGCGACAACTGCCCGTTTCCGACCTGCCGCGGGCAGAACTGGGGAACTTCCGTCATTTTCATCCGaaa                                                >  1:402170/10‑150 (MQ=255)                                                                                                   aCCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGcc                                      >  1:106299/1‑151 (MQ=255)                                                                                                    ccGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCa                                     <  2:300902/151‑1 (MQ=255)                                                                                                      gACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCAcc                                   >  2:582019/1‑151 (MQ=255)                                                                                                             ccTGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCgg                                >  1:543855/1‑147 (MQ=255)                                                                                                               tGGTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTg                          >  2:183469/1‑151 (MQ=255)                                                                                                                 gTATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGc                        >  2:476569/1‑151 (MQ=255)                                                                                                                  tATGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCt                       >  2:320039/1‑151 (MQ=255)                                                                                                                   aTGATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCt                       >  2:689176/1‑150 (MQ=255)                                                                                                                     gATCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCtat                     >  2:748615/1‑150 (MQ=255)                                                                                                                      aTCGCGACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTCCAGTTTGCTCCATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCtat                     <  1:92970/149‑1 (MQ=255)                                                                                                                         gcgACCGCTTTATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGcc                >  1:548690/1‑151 (MQ=255)                                                                                                                                   tATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGCCCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCtata                    <  1:571219/137‑1 (MQ=255)                                                                                                                                   tATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCtata                    >  2:571257/1‑137 (MQ=255)                                                                                                                                   tATTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGa       >  1:280677/1‑150 (MQ=255)                                                                                                                                    aTTCTTTCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGaa      <  2:825901/150‑1 (MQ=255)                                                                                                                                          tCCAGCGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAacca  <  1:65928/148‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                ATGCGGTACAAAAAGCCAACTCTGGTCATCCCGGCGCGCCGATGGGCATGGCTGATATTGCCGAAGTGCTGTGGAACGATTTTCTTAAACATAACCCTACCGACCCAACCTGGTATGATCGCGACCGCTTTATTCTTTCCAACGGTCACGCGTCGATGCTGCTCTACAGTTTGCTACATCTGACCGGTTACGACCTGCCGCTGGAAGAACTGAAGAACTTCCGTCAGTTGCATTCGAAAACCCCAGGCCACCCGGAGATTGGCTATACGCCAGGCGTTGAAACCA  >  NC_000913/2579682‑2579966

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