Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	125,814	C→T	D40D (GAC→GAT)	aceF →	pyruvate dehydrogenase, E2 subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	125,814	0	C	T	100.0%	129.0 / NA	37	D40D (GAC→GAT)	aceF	pyruvate dehydrogenase, E2 subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/21);  total (16/21)

AAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACC  >  NC_000913/125679‑125948                                                                                                                                        |                                                                                                                                       aaGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAg                                                                                                                                    <  2:451948/140‑1 (MQ=255)  agagGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGc                                                                                                                                   <  1:563695/140‑1 (MQ=255)  agagGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGc                                                                                                                                   <  2:484006/140‑1 (MQ=255)       tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATaaa                                                                                                                                     <  2:745980/133‑1 (MQ=255)       tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATaaa                                                                                                                                     >  1:745980/1‑133 (MQ=255)       tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGAGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTa                                                                                                                              >  1:370269/1‑140 (MQ=255)             aataatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGc                                                                                                                                   >  1:1174609/1‑129 (MQ=255)             aataatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGc                                                                                                                                   <  2:1174609/129‑1 (MQ=255)              ataatGGCTATCGAAATC‑AAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGtt                                                                                                                      <  2:1031805/140‑1 (MQ=255)              ataatGGCTATCGAAATC‑AAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGtt                                                                                                                      >  1:1031805/1‑140 (MQ=255)               taatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGtt                                                                                                                      >  2:389474/1‑140 (MQ=255)                   ggCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGt                                                                                                                  <  2:732002/140‑1 (MQ=255)                    gCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGtc                                                                                                                 >  2:604783/1‑140 (MQ=255)                       aTCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTcc                                                                                                              >  1:294871/1‑140 (MQ=255)                                aaaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg                                                                                                     <  1:389474/140‑1 (MQ=255)                                aaaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg                                                                                                     <  2:294871/140‑1 (MQ=255)                                aaaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg                                                                                                     <  1:159487/140‑1 (MQ=255)                                   gTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTAt                                                                                                  <  1:473899/140‑1 (MQ=255)                                           cATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAag                                                                                          <  1:461591/140‑1 (MQ=255)                                             tCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAagag                                                                                        >  1:773830/1‑140 (MQ=255)                                                 ggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCa                                                                                    <  1:604783/140‑1 (MQ=255)                                                      atgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCg                                                                                                <  2:520359/123‑1 (MQ=255)                                                      atgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCg                                                                                                >  1:520359/1‑123 (MQ=255)                                                      atgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGtc                                                                               >  1:219525/1‑140 (MQ=255)                                                      atgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGtc                                                                               >  1:429278/1‑140 (MQ=255)                                                                 aTCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGa                                                                    >  2:935826/1‑140 (MQ=255)                                                                  tCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGAt                                                                   <  2:219525/140‑1 (MQ=255)                                                                  tCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGAt                                                                   <  1:538476/140‑1 (MQ=255)                                                                      cGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGGATCGTTAAAAAGATCAAAATCTCTGTTGGCGATaaaa                                                               >  1:504911/1‑140 (MQ=255)                                                                                     aGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCAc                                                <  2:429278/140‑1 (MQ=255)                                                                                     aGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAAAGATCAAAGCCTCTGTTGGCGATAAAACCCAAACCCGCGCAc                                                >  2:170438/1‑140 (MQ=37)                                                                                          gCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGAtt                                           <  1:318457/140‑1 (MQ=255)                                                                                                                 tCGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACgg                    <  1:170438/140‑1 (MQ=255)                                                                                                                  cGCTGATCACCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGt                   <  2:376276/140‑1 (MQ=255)                                                                                                                          aCCGTAGAAGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGGATGCTTTAAAAGAGCAAAGTAGCTGATGGCGAGAAAACCCCCACCGGCGCACTGATTATGATTTTCGATTCCCCCGACCGGGCAGCAGa           >  2:892382/1‑140 (MQ=255)                                                                                                                                  aGGCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCAc   <  2:1175558/140‑1 (MQ=255)                                                                                                                                   ggCGATAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGGATCGTTAAAAAGATCAAAAGCTCTGTTGGCGATAAAAACCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCAcc  >  2:762985/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       AAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACGGTGCAGCAGACGCTGCACC  >  NC_000913/125679‑125948

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