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%	188.9 / NA	52	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 (28/24);  total (28/24)

AACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATG  >  NC_000913/125102‑125347                                                                                                                              |                                                                                                                         aaCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                             <  2:285696/139‑1 (MQ=255)     aCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGa                                                                                                               <  1:687981/133‑1 (MQ=255)     aCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGa                                                                                                               >  2:687981/1‑133 (MQ=255)     aCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagcagc                                                                                                         >  1:991461/1‑139 (MQ=255)            gCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                  >  2:872727/1‑139 (MQ=255)              cattgccggcGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCt                                                                                                >  1:360685/10‑139 (MQ=255)              cGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCt                                                                                                >  1:309189/1‑139 (MQ=255)                  aaTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                            <  2:991461/139‑1 (MQ=255)                  aaTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                            <  2:309189/139‑1 (MQ=255)                    tGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGaa                                                                                          >  1:1254732/1‑139 (MQ=255)                              tGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGc                                                                                >  2:950569/1‑139 (MQ=255)                                cTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATcc                                                                                                 <  2:89479/120‑1 (MQ=255)                                cTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATcc                                                                                                 >  1:89479/1‑120 (MQ=255)                                       aGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                       <  2:158471/139‑1 (MQ=255)                                        ggTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTg                                                                      >  1:1349759/1‑139 (MQ=255)                                          tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                         <  1:101118/134‑1 (MQ=255)                                          tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                         >  2:101118/1‑134 (MQ=255)                                             ccGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTgg                                                                                              >  1:509387/1‑110 (MQ=255)                                             ccGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTgg                                                                                              <  2:509387/110‑1 (MQ=255)                                               gTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTta                                                               >  1:708398/1‑139 (MQ=255)                                                    ggTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCg                                                          <  2:708398/139‑1 (MQ=255)                                                       aTCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGa                                                       <  1:872727/139‑1 (MQ=255)                                                               aCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTc                                               >  2:728237/1‑139 (MQ=255)                                                                   gAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCg                                           >  2:58575/1‑139 (MQ=255)                                                                             aaGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGcgcg                                 <  2:1254732/139‑1 (MQ=255)                                                                             aaGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGcgcg                                 <  2:1349759/139‑1 (MQ=255)                                                                                  aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                             >  2:315820/1‑58 (MQ=255)                                                                                  aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                             <  1:315820/58‑1 (MQ=255)                                                                                    cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                             <  2:671422/56‑1 (MQ=255)                                                                                    cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                             >  1:671422/1‑56 (MQ=255)                                                                                    cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                   >  2:113094/1‑82 (MQ=255)                                                                                    cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                   <  1:113094/82‑1 (MQ=255)                                                                                    cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                          >  1:221075/1‑139 (MQ=255)                                                                                    cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                          >  1:479202/1‑139 (MQ=255)                                                                                    cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                          <  2:479202/139‑1 (MQ=255)                                                                                      aaGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCa                        <  1:785641/139‑1 (MQ=255)                                                                                       aGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                             <  2:231758/53‑1 (MQ=255)                                                                                       aGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                             >  1:231758/1‑53 (MQ=255)                                                                                       aGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                  <  2:1093981/64‑1 (MQ=255)                                                                                       aGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                  >  1:1093981/1‑64 (MQ=255)                                                                                         gTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATcc                                                                                                 <  1:762552/63‑1 (MQ=255)                                                                                         gTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATcc                                                                                                 >  2:762552/1‑63 (MQ=255)                                                                                          tAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGAt                    >  2:754536/1‑139 (MQ=255)                                                                                            aaGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCa                                                                                                                        >  2:121291/1‑37 (MQ=255)                                                                                            aaGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCa                                                                                                                        <  1:121291/37‑1 (MQ=255)                                                                                              gTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgtg                <  1:950569/139‑1 (MQ=255)                                                                                                   gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                         <  1:817086/77‑1 (MQ=255)                                                                                                   gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                         >  2:817086/1‑77 (MQ=255)                                                                                                   gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGc           >  1:1037963/1‑139 (MQ=255)                                                                                                        tCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGaa      <  2:221075/139‑1 (MQ=255)                                                                                                        tCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGaa      >  2:896547/1‑139 (MQ=255)                                                                                                            cTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACAtg  >  2:566999/1‑139 (MQ=255)                                                                                                                              |                                                                                                                         AACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATG  >  NC_000913/125102‑125347

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