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%	221.1 / NA	61	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/32);  total (29/32)

GCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCT  >  NC_000913/125092‑125358                                                                                                                                        |                                                                                                                                    gCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCAc                                                                                                                                  <  2:1501292/139‑1 (MQ=255) gCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCAc                                                                                                                                  <  2:141148/139‑1 (MQ=255)   tGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGt                                                                                                                                >  2:636894/1‑139 (MQ=255)     aaCGAAAAGTACCACATGCCGGCACTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAGGGTATCTACAAACTCGAAACTATTGAAGGTAGCTAAGGTAAAGTTCAGCTGCTCGACTCCGTTTCTATCCTGCGCCACGTcc                                                                                                                              >  2:258178/1‑139 (MQ=255)     aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTcc                                                                                                                              >  2:826829/1‑139 (MQ=255)     aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTcc                                                                                                                              >  1:361704/1‑139 (MQ=255)     aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTc                                                                                                                               <  2:392920/138‑1 (MQ=255)     aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTc                                                                                                                               >  1:392920/1‑138 (MQ=255)     aaCGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTG‑AGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCg                                                                                                                             >  1:1455767/1‑139 (MQ=255)        gAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTg                                                                                                                           <  1:826829/139‑1 (MQ=255)           aaCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                        <  1:636894/139‑1 (MQ=255)           aaCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                        <  1:820224/139‑1 (MQ=255)            aCTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagc                                                                                                                       <  2:1468693/139‑1 (MQ=255)                 cacaTGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagc                                                                                                                       <  2:91877/134‑1 (MQ=255)                 cacaTGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAagc                                                                                                                       >  1:91877/1‑134 (MQ=255)                      gCCGGCAATGCCGGAAGGTTCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                             >  1:1454542/1‑139 (MQ=255)                          gCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTgg                                                                                                         <  1:43479/139‑1 (MQ=255)                            aaTGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                       <  2:361704/139‑1 (MQ=255)                                     aGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                              <  2:1454542/139‑1 (MQ=255)                                          cTGAGGAAGGTATCCGTAAAGGTATCTAGAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGt                                                                                         >  1:951595/1‑139 (MQ=255)                                           tGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCAc                                                                                                                                  >  1:1009043/1‑97 (MQ=255)                                           tGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCAc                                                                                                                                  <  2:1009043/97‑1 (MQ=255)                                                 aGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCt                                                                                  <  1:337246/139‑1 (MQ=255)                                                  ggTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCa                                                                                                                                   <  1:239139/89‑1 (MQ=255)                                                  ggTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCa                                                                                                                                   >  2:239139/1‑89 (MQ=255)                                                  ggTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTACGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTg                                                                                 >  1:816864/1‑139 (MQ=255)                                                    tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                        <  1:1458979/98‑1 (MQ=255)                                                    tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                        >  2:1458979/1‑98 (MQ=255)                                                    tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                             >  1:567596/1‑109 (MQ=255)                                                    tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                             <  2:567596/109‑1 (MQ=255)                                                    tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                       >  1:626136/1‑115 (MQ=255)                                                    tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                       <  2:626136/115‑1 (MQ=255)                                                    tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAgg                                                                                      <  2:302690/132‑1 (MQ=255)                                                    tATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAgg                                                                                      >  1:302690/1‑132 (MQ=255)                                                     aTCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTg                                                                                                                           <  1:395123/94‑1 (MQ=255)                                                     aTCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTg                                                                                                                           >  2:395123/1‑94 (MQ=255)                                                                tATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTg                                                                   <  2:611440/139‑1 (MQ=255)                                                                         aCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcctc                                                              >  1:1324187/1‑135 (MQ=255)                                                                         aCTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGAcctc                                                              <  2:1324187/135‑1 (MQ=255)                                                                            cGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                                    <  2:209840/110‑1 (MQ=255)                                                                            cGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGtt                                                                                    >  1:209840/1‑110 (MQ=255)                                                                                    ttGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                             >  2:24004/1‑77 (MQ=255)                                                                                    ttGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATc                                                                                                             <  1:24004/77‑1 (MQ=255)                                                                                     tGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtata                                                                        >  1:291639/1‑113 (MQ=255)                                                                                     tGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTtata                                                                        <  2:291639/113‑1 (MQ=255)                                                                                     tGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGgcg                                              <  2:657446/139‑1 (MQ=255)                                                                                         ggTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGa                                                                                                      <  1:629457/79‑1 (MQ=255)                                                                                         ggTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGa                                                                                                      >  2:629457/1‑79 (MQ=255)                                                                                         ggTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTg                                          <  2:1455767/139‑1 (MQ=255)                                                                                           tAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGAt                                        <  2:816864/139‑1 (MQ=255)                                                                                            aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGt                                                                                         <  1:360415/89‑1 (MQ=255)                                                                                            aGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGt                                                                                         >  2:360415/1‑89 (MQ=255)                                                                                                 aGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                              >  1:704966/1‑79 (MQ=255)                                                                                                 aGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTAc                                                                                              <  2:704966/79‑1 (MQ=255)                                                                                                        gTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                       <  1:143747/63‑1 (MQ=255)                                                                                                        gTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCg                                                                                                       >  2:143747/1‑63 (MQ=255)                                                                                                             gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                        >  2:453408/1‑41 (MQ=255)                                                                                                             gctgctCGGCTCCGGTTCTATCCTGCGCCACGTCCGTGAag                                                                                                                        <  1:453408/41‑1 (MQ=255)                                                                                                                                aTCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGc   >  1:142258/1‑139 (MQ=255)                                                                                                                                aTCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGc   >  2:857291/1‑139 (MQ=255)                                                                                                                                 tCCTGCGCCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCt  >  1:612103/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                    GCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCTGCGTCACGTCCGTGAAGCAGCTGAGATCCTGGCGAAAGATTACGGCGTAGGTTCTGACGTTTATAGCGTGACCTCCTTCACCGAGCTGGCGCGTGATGGTCAGGATTGTGAACGCTGGAACATGCTGCACCCGCT  >  NC_000913/125092‑125358

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