Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	126,432	T→C	A246A (GCT→GCC)	aceF →	pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	126,432	0	T	C	100.0%	218.2 / NA	57	A246A (GCT→GCC)	aceF	pyruvate dehydrogenase, dihydrolipoyltransacetylase component E2
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/26);  total (31/26)

GCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGCAGCTCCTG  >  NC_000913/126301‑126562                                                                                                                                    |                                                                                                                                   gCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATgg                                                                                                                             >  2:1288274/1‑139 (MQ=255) gCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATgg                                                                                                                             >  2:1135619/1‑139 (MQ=255)    gCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAg                                                                                                                          <  2:159813/139‑1 (MQ=255)            gAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGgcg                                                                                                                  <  2:1041485/139‑1 (MQ=255)                      aaCGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAg                                                                                                        <  2:1064858/139‑1 (MQ=255)                           tCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGcgtc                                                                                                   >  1:672561/1‑139 (MQ=255)                            ccGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGcgtcg                                                                                                  >  1:126161/1‑139 (MQ=255)                                          tGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGaaa                                                                                    >  1:936654/1‑139 (MQ=255)                                              gAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCa                                                                                <  2:672561/139‑1 (MQ=255)                                                  ttGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGt                                                                            <  2:929649/139‑1 (MQ=255)                                                         gACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGAt                                                                     <  1:1288274/139‑1 (MQ=255)                                                          aCTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGc                                                                                                      <  1:1008651/105‑1 (MQ=255)                                                          aCTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGc                                                                                                      >  2:1008651/1‑105 (MQ=255)                                                          aCTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATa                                                                    >  2:1294899/1‑139 (MQ=255)                                                                  gATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAACGCCTCTATGGAAGTTCCGGCACCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGaaa                                                            >  1:569509/1‑139 (MQ=255)                                                                   aTGGTGAAAGTGGGCGATAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGaaaa                                                           <  2:685840/139‑1 (MQ=255)                                                                    tGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGgcg                                                                                                                  >  1:1014413/1‑83 (MQ=255)                                                                    tGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGgcg                                                                                                                  <  2:1014413/83‑1 (MQ=255)                                                                          aaGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTc                                                    <  1:255203/139‑1 (MQ=255)                                                                          aaGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATCGAAGTTCCGGCACCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGccc                                                    <  2:569509/139‑3 (MQ=255)                                                                           ggTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                                                   >  2:916495/2‑139 (MQ=255)                                                                           aGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCTCCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                                                   <  2:936654/139‑1 (MQ=255)                                                                           aGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                                                   <  2:633046/139‑1 (MQ=255)                                                                             tGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAAGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCt                                                 >  1:169114/1‑139 (MQ=255)                                                                               ggCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGAt                                                                     >  1:157895/1‑117 (MQ=255)                                                                               ggCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGAt                                                                     <  2:157895/117‑1 (MQ=255)                                                                                 cGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGAtt                                             >  2:1166277/1‑139 (MQ=255)                                                                                       aGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTAt                                           <  2:43948/135‑1 (MQ=255)                                                                                       aGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTAt                                           >  1:43948/1‑135 (MQ=255)                                                                                            ccGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                                                   <  2:118168/122‑1 (MQ=255)                                                                                            ccGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCg                                                   >  1:118168/1‑122 (MQ=255)                                                                                              gCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAg                                                                                                        <  1:282869/67‑1 (MQ=255)                                                                                              gCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAg                                                                                                        >  2:282869/1‑67 (MQ=255)                                                                                                         aCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCa                     <  2:306830/139‑1 (MQ=255)                                                                                                            gATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGa                                                                                               >  2:244785/1‑62 (MQ=255)                                                                                                            gATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGa                                                                                               <  1:244785/62‑1 (MQ=255)                                                                                                            gATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGa                            >  2:941896/1‑129 (MQ=255)                                                                                                            gATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGa                            <  1:941896/129‑1 (MQ=255)                                                                                                               cACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGc                                                                                                      <  1:1339031/52‑1 (MQ=255)                                                                                                               cACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGc                                                                                                      >  2:1339031/1‑52 (MQ=255)                                                                                                                 ccGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGaa                                                                                              >  2:1420357/1‑58 (MQ=255)                                                                                                                 ccGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGaa                                                                                              <  1:1420357/58‑1 (MQ=255)                                                                                                                 ccGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTc                                                    >  2:355499/1‑100 (MQ=255)                                                                                                                 ccGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTc                                                    <  1:355499/100‑1 (MQ=255)                                                                                                                 ccGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGc             >  1:906346/1‑139 (MQ=255)                                                                                                                 ccGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGc             >  2:1064255/1‑139 (MQ=255)                                                                                                                 ccGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGA‑TATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCg            >  1:1441012/1‑139 (MQ=255)                                                                                                                       aaGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGCAGc       >  2:347976/1‑139 (MQ=255)                                                                                                                       aaGGCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGCAGc       >  1:479677/1‑139 (MQ=255)                                                                                                                         ggCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACt                                                                                        <  2:1009881/56‑1 (MQ=255)                                                                                                                         ggCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACt                                                                                        >  1:1009881/1‑56 (MQ=255)                                                                                                                         ggCGACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGCAGCTc     >  1:1207832/1‑139 (MQ=255)                                                                                                                            gACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGt                                                                            <  2:1074869/65‑1 (MQ=255)                                                                                                                            gACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGt                                                                            >  1:1074869/1‑65 (MQ=255)                                                                                                                            gACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCgg           >  2:333911/1‑130 (MQ=255)                                                                                                                            gACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCgg           <  1:333911/130‑1 (MQ=255)                                                                                                                            gACAAAGCCTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGCAGCTCCTg  >  2:334423/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                   GCGGCTGGCGTGAAAGAAGTTAACGTTCCGGATATCGGCGGTGACGAAGTTGAAGTGACTGAAGTGATGGTGAAAGTGGGCGACAAAGTTGCCGCTGAACAGTCACTGATCACCGTAGAAGGCGACAAAGCTTCTATGGAAGTTCCGGCGCCGTTTGCAGGCGTCGTGAAGGAACTGAAAGTCAACGTTGGCGATAAAGTGAAAACTGGCTCGCTGATTATGATCTTCGAAGTTGAAGGCGCAGCGCCTGCGGCAGCTCCTG  >  NC_000913/126301‑126562

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: ‑