Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	125,804	T→C	V37A (GTA→GCA)	aceF →	pyruvate dehydrogenase, E2 subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	125,804	0	T	C	100.0%	155.6 / NA	43	V37A (GTA→GCA)	aceF	pyruvate dehydrogenase, E2 subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (18/25);  total (18/25)

GTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACG  >  NC_000913/125670‑125929                                                                                                                                       |                                                                                                                              gTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAg                                                                                                                          >  1:695633/1‑140 (MQ=255) gTCTGGCGTAAGAGGGAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAg                                                                                                                          >  2:714586/1‑140 (MQ=255)  tCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAgg                                                                                                                         >  1:308672/1‑140 (MQ=255)      gCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGAc                                                                                                                     <  1:1059182/140‑1 (MQ=255)        gTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACaa                                                                                                                   <  2:177197/140‑1 (MQ=255)         tAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACaaa                                                                                                                  <  1:245985/140‑1 (MQ=255)           agagGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGc                                                                                                                <  2:81549/140‑1 (MQ=255)              ggTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCctc                                                                                                             <  2:638835/140‑1 (MQ=255)              ggTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCctc                                                                                                             <  1:394273/140‑1 (MQ=255)              ggTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCctc                                                                                                             <  1:829883/140‑1 (MQ=255)                tAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTa                                                                                                           >  2:147821/1‑140 (MQ=255)                tAAAAGAACAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTa                                                                                                           >  1:469755/1‑140 (MQ=255)                 aaaaGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGc                                                                                                                        >  1:942226/1‑126 (MQ=255)                 aaaaGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGc                                                                                                                        <  2:942226/126‑1 (MQ=255)                        taatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACaaa                                                                                                                  >  2:479141/1‑125 (MQ=255)                        taatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACaaa                                                                                                                  <  1:479141/125‑1 (MQ=255)                        taatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGtt                                                                                                   >  2:374781/1‑140 (MQ=255)                         aatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTc                                                                                                  >  1:410739/1‑140 (MQ=255)                         aatGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTc                                                                                                  >  2:297936/1‑140 (MQ=255)                          atGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTcc                                                                                                 >  1:770507/1‑140 (MQ=255)                                aTCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTcc                                                                                           >  2:668655/1‑140 (MQ=255)                                aTCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTcc                                                                                           >  2:302659/1‑140 (MQ=255)                                 tCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCg                                                                                          <  2:308672/140‑1 (MQ=255)                                  cGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGc                                                                                         >  2:531364/1‑140 (MQ=255)                                    aaaTCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAg                                                                                       <  2:770507/140‑1 (MQ=255)                                    aaaTCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAg                                                                                       <  1:820979/140‑1 (MQ=255)                                         aaaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg                                                                                  <  1:714586/140‑1 (MQ=255)                                         aaaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg                                                                                  <  1:374781/140‑1 (MQ=255)                                         aaaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg                                                                                  <  1:668655/140‑1 (MQ=255)                                         aaaGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCggg                                                                                  <  1:725386/140‑1 (MQ=255)                                           aGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTa                                                                                <  2:481965/140‑1 (MQ=255)                                               ccGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGt                                                                            >  2:78605/1‑140 (MQ=255)                                                    cATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAag                                                                       <  2:410739/140‑1 (MQ=255)                                                       cGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGGATCGTTaaaaaaa                                                                    >  2:800687/1‑136 (MQ=255)                                                        ggggCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAGGAGAt                                                                   <  2:665237/140‑1 (MQ=255)                                                                tgaAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGtct                                                           <  1:705438/140‑1 (MQ=255)                                                                      tGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTg                                                     <  1:297936/140‑1 (MQ=255)                                                                            cACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATa                                               >  1:284894/1‑140 (MQ=255)                                                                                           cAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGgcg                                <  1:78605/140‑1 (MQ=255)                                                                                             aaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCa                              <  1:800687/140‑1 (MQ=255)                                                                                             aaGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGACTCTGTTGGCGATAAAAACCAGACCGgcgcc                              >  1:628525/1‑139 (MQ=255)                                                                                                   gCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAATCTCTGTTGGCGGTAAAAACCAAACCGGCGCACTGAtt                        >  2:134519/1‑140 (MQ=37)                                                                                                                aGCCGAACAGTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGAtt           <  1:134519/140‑1 (MQ=255)                                                                                                                         gTCGCTGATCACCGCAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACg  <  2:844628/140‑1 (MQ=255)                                                                                                                                       |                                                                                                                              GTCTGGCGTAAGAGGTAAAAGAATAATGGCTATCGAAATCAAAGTACCGGACATCGGGGCTGATGAAGTTGAAATCACCGAGATCCTGGTCAAAGTGGGCGACAAAGTTGAAGCCGAACAGTCGCTGATCACCGTAGAAGGCGACAAAGCCTCTATGGAAGTTCCGTCTCCGCAGGCGGGTATCGTTAAAGAGATCAAAGTCTCTGTTGGCGATAAAACCCAGACCGGCGCACTGATTATGATTTTCGATTCCGCCGACG  >  NC_000913/125670‑125929

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