Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	125,227	T→C	R737R (CGT→CGC)	aceE →	pyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	125,227	0	T	C	100.0%	146.5 / NA	41	R737R (CGT→CGC)	aceE	pyruvate dehydrogenase, decarboxylase component E1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (20/21);  total (20/21)

CTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGC  >  NC_000913/125093‑125357                                                                                                                                       |                                                                                                                                   cTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACg                                                                                                                                >  2:159061/1‑139 (MQ=255) cTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACg                                                                                                                                >  1:538768/1‑139 (MQ=255)    aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTcc                                                                                                                             >  2:43465/1‑139 (MQ=255)    aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTcc                                                                                                                             >  2:577505/1‑139 (MQ=255)    aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTcc                                                                                                                             >  2:696521/1‑139 (MQ=255)     aCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCg                                                                                                                            <  2:300642/139‑1 (MQ=255)     aCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCg                                                                                                                            >  1:397434/1‑139 (MQ=255)       gAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTg                                                                                                                          <  1:543619/139‑1 (MQ=255)            cTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagca                                                                                                                     <  1:257632/139‑1 (MQ=255)                   aTGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTgaga                                                                                                              >  2:676648/1‑139 (MQ=255)                           aaTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                      <  2:399526/139‑1 (MQ=255)                           aaTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                      <  1:696521/139‑1 (MQ=255)                                    aGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                             <  2:397434/139‑1 (MQ=255)                                       tGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGc                                                                                          >  2:62783/1‑139 (MQ=255)                                         cTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGt                                                                                        >  1:726815/1‑139 (MQ=255)                                         cTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGt                                                                                        <  1:319921/139‑1 (MQ=255)                                                    aTCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACg                                                                             >  2:526992/1‑139 (MQ=255)                                                      ccGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGtt                                                                           >  2:675432/1‑139 (MQ=255)                                                        gTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTta                                                                         >  1:548862/1‑139 (MQ=255)                                                             ggTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCg                                                                    <  1:62783/139‑1 (MQ=255)                                                             ggTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCg                                                                    <  1:252531/139‑1 (MQ=255)                                                               tATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTg                                                                  <  1:675432/139‑1 (MQ=255)                                                                aTCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGa                                                                 >  2:768982/1‑139 (MQ=255)                                                                   tACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcct                                                              >  1:302449/1‑139 (MQ=255)                                                                       aaCTCGAAACTATTGAAGGTAGCAAAGGTCAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                     <  2:486215/128‑1 (MQ=255)                                                                       aaCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGc                                                                     >  1:486215/1‑128 (MQ=255)                                                                                        ggTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                         <  1:82785/139‑1 (MQ=255)                                                                                           aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                                 <  2:641665/96‑1 (MQ=255)                                                                                           aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                                 >  1:641665/1‑96 (MQ=255)                                                                                           aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                         <  2:396369/136‑1 (MQ=255)                                                                                           aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                         >  1:396369/1‑136 (MQ=255)                                                                                           aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATg                                      >  2:644891/1‑139 (MQ=255)                                                                                                    aaaGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGAtt                             <  1:526992/139‑1 (MQ=255)                                                                                                       gTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGCCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgtg                          <  2:622899/139‑1 (MQ=255)                                                                                                            gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGa                                                                               >  1:89353/1‑81 (MQ=255)                                                                                                            gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGa                                                                               <  2:89353/81‑1 (MQ=255)                                                                                                              tgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTg                   <  2:371666/139‑1 (MQ=255)                                                                                                                         ggTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACAtgctgc        <  1:210430/139‑1 (MQ=255)                                                                                                                         ggTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACAtgctgc        <  1:644891/139‑1 (MQ=255)                                                                                                                             cTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCAccc    <  1:760688/139‑1 (MQ=255)                                                                                                                               aTCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGc  >  1:662015/1‑139 (MQ=255)                                                                                                                                       |                                                                                                                                   CTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGC  >  NC_000913/125093‑125357

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