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%	226.8 / NA	63	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 (29/34);  total (29/34)

CGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCC  >  NC_000913/125091‑125355                                                                                                                                         |                                                                                                                                 cGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCa                                                                                                                                <  1:837800/139‑1 (MQ=255) cGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCa                                                                                                                                >  1:1440629/1‑139 (MQ=255)  gCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCAc                                                                                                                               <  1:1497223/139‑1 (MQ=255)           aaaCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACATCCGTGaa                                                                                                                      >  2:465173/1‑139 (MQ=255)                aCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                     <  2:318924/135‑1 (MQ=255)                aCCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                     >  1:318924/1‑135 (MQ=255)                               tGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                     >  1:1100453/1‑120 (MQ=255)                               tGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                     <  2:1100453/120‑1 (MQ=255)                                                gAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                                 <  1:1403920/139‑1 (MQ=255)                                                  aGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGt                                                                                      <  2:1654024/132‑1 (MQ=255)                                                  aGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGt                                                                                      >  1:1654024/1‑132 (MQ=255)                                                  aGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                               <  2:947538/139‑1 (MQ=255)                                                        ccGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGa                                                                                                                       <  2:532631/93‑1 (MQ=255)                                                        ccGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGa                                                                                                                       >  1:532631/1‑93 (MQ=255)                                                             aaGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGa                                                                                                                       <  2:764112/88‑1 (MQ=255)                                                             aaGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGa                                                                                                                       >  1:764112/1‑88 (MQ=255)                                                                   tCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                     <  1:1407268/84‑1 (MQ=255)                                                                   tCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                     >  2:1407268/1‑84 (MQ=255)                                                                    cTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCg                                                                                                                          >  1:1228232/1‑78 (MQ=255)                                                                    cTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCg                                                                                                                          <  2:1228232/78‑1 (MQ=255)                                                                     tACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcct                                                            >  2:644293/1‑139 (MQ=255)                                                                          aCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTc                                                       >  2:589883/1‑139 (MQ=255)                                                                            tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtat                                                                      >  1:111014/1‑122 (MQ=255)                                                                            tCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtat                                                                      <  2:111014/122‑1 (MQ=255)                                                                               aaaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATcc                                                                                                         >  2:863305/1‑84 (MQ=255)                                                                               aaaCTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATcc                                                                                                         <  1:863305/84‑1 (MQ=255)                                                                                     ttGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGgc                                            <  1:644293/139‑1 (MQ=255)                                                                                          ggTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                       <  1:465173/139‑1 (MQ=255)                                                                                             aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCAc                                                                                                                               >  1:146250/1‑48 (MQ=255)                                                                                             aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCAc                                                                                                                               <  2:146250/48‑1 (MQ=255)                                                                                             aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                       >  2:433939/1‑136 (MQ=255)                                                                                             aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                       <  1:433939/136‑1 (MQ=255)                                                                                               cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTgg                                                                                                      <  2:1206572/71‑1 (MQ=255)                                                                                               cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTgg                                                                                                      >  1:1206572/1‑71 (MQ=255)                                                                                               cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                    <  1:88468/73‑1 (MQ=255)                                                                                               cAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                    >  2:88468/1‑73 (MQ=255)                                                                                                aaaGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTc                                 <  2:1440629/139‑1 (MQ=255)                                                                                                 aaGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCa                                <  1:589883/139‑1 (MQ=255)                                                                                                    gTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                                  <  1:1389836/134‑1 (MQ=255)                                                                                                    gTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGt                                  >  2:1389836/1‑134 (MQ=255)                                                                                                     tAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtata                                                                     >  1:1142279/1‑98 (MQ=255)                                                                                                     tAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtata                                                                     <  2:1142279/98‑1 (MQ=255)                                                                                                        aGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATtgt                         <  1:1305971/139‑1 (MQ=255)                                                                                                             aGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                               <  2:1351874/80‑1 (MQ=255)                                                                                                             aGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                               >  1:1351874/1‑80 (MQ=255)                                                                                                             aGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtat                                                                      <  1:666248/89‑1 (MQ=255)                                                                                                             aGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtat                                                                      >  2:666248/1‑89 (MQ=255)                                                                                                             aGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGACCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                               <  2:199081/80‑1 (MQ=255)                                                                                                             aGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGACCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                               >  1:199081/1‑80 (MQ=255)                                                                                                              gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCg                                                                  >  1:1046871/1‑92 (MQ=255)                                                                                                              gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCg                                                                  <  2:1046871/92‑1 (MQ=255)                                                                                                              gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGc                   >  2:1229235/1‑139 (MQ=255)                                                                                                               ctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                    <  2:513144/57‑1 (MQ=255)                                                                                                               ctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                    >  1:513144/1‑57 (MQ=255)                                                                                                               ctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGa                                                                                               >  1:693171/1‑62 (MQ=255)                                                                                                               ctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGa                                                                                               <  2:693171/62‑1 (MQ=255)                                                                                                               ctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGTCGTAGGTTCTGACGTGTATAGCGTGACCTCCTTCACCGAGCTGgc                                            <  1:963057/113‑1 (MQ=255)                                                                                                               ctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGgc                                            >  2:963057/1‑113 (MQ=255)                                                                                                                  ctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                           >  2:130072/1‑63 (MQ=255)                                                                                                                  ctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                           <  1:130072/63‑1 (MQ=255)                                                                                                                            gTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                                 <  2:6517/63‑1 (MQ=255)                                                                                                                            gTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                                 >  1:6517/1‑63 (MQ=255)                                                                                                                               cTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCAccc  <  1:1023298/139‑1 (MQ=255)                                                                                                                                         |                                                                                                                                 CGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCC  >  NC_000913/125091‑125355

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