Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	126,807	T→C	R371R (CGT→CGC)	aceF →	pyruvate dehydrogenase, E2 subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	126,807	0	T	C	100.0%	115.5 / NA	33	R371R (CGT→CGC)	aceF	pyruvate dehydrogenase, E2 subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/9);  total (24/9)

TTCACGCGACTCCGCTGATCCGCCGTCTGGCACGCGAGTTTGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGTGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTG  >  NC_000913/126674‑126934                                                                                                                                      |                                                                                                                                ttCACGCGACTCCGCTGATCCGCCGTCTGGCACGCGAGTTTGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGaa                                                                                                                           <  1:1007731/140‑1 (MQ=255)   cACGCGACTCCGCTGATCCGCCGTCTGGCACGCGAGTTTGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAagc                                                                                                                         <  1:253179/140‑1 (MQ=255)                    ccgccgTCTGGCACGCGAGTTTGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTg                                                                                                        >  1:433365/1‑140 (MQ=255)                     cgccgTCTGGCACGCGAGTTTGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTgg                                                                                                       <  2:793048/140‑1 (MQ=255)                       ccgTCTGGCACGCGAGTTTGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGAc                                                                                                          >  1:1045360/1‑135 (MQ=255)                       ccgTCTGGCACGCGAGTTTGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGAc                                                                                                          <  2:1045360/135‑1 (MQ=255)                              gCACGCGAGTTTGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTAt                                                                                              >  1:184128/1‑140 (MQ=255)                                       tttGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCAt                                                                                     <  2:433365/140‑1 (MQ=255)                                              ttAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGGAACCCTGGCATGCTGccc                                                                              >  2:714432/1‑139 (MQ=255)                                                aaCCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTg                                                                            >  1:711542/1‑140 (MQ=255)                                                     tGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGa                                                                       >  1:268070/1‑140 (MQ=255)                                                      gCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGaa                                                                      >  1:671688/1‑140 (MQ=255)                                                      gCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGGGGTATCCCCGGCATGCTGCCGTGGCCGaa                                                                      >  1:149451/1‑140 (MQ=255)                                                          aaGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCCCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTg                                                                  <  2:803095/140‑1 (MQ=255)                                                             tGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGAc                                                               >  1:832455/1‑140 (MQ=255)                                                                    cACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTTGCCGAAGGTGGACTTCAGCa                                                        >  1:789137/1‑140 (MQ=255)                                                                      cTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGa                                                                       <  1:187939/123‑1 (MQ=255)                                                                      cTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGa                                                                       >  2:187939/1‑123 (MQ=255)                                                                      cTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAg                                                      >  1:171163/1‑140 (MQ=255)                                                                       tGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGGATCCCCGGCATGCTGCCCGGGCCGAAGGTGGACTTCAGCAAgg                                                     >  2:477729/1‑139 (MQ=255)                                                                                     gTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGaaa                                            >  1:410158/1‑135 (MQ=255)                                                                                     gTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGaaa                                            <  2:410158/135‑1 (MQ=255)                                                                                     gTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGGGGACTTCAGCAAGTTTGGTGAAATCgaa                                       >  2:915386/1‑140 (MQ=255)                                                                                            gcgcGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCCGGCATGCTGCCGCGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTgg                                >  1:888128/1‑140 (MQ=255)                                                                                             cgcgAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGa                               >  1:976767/1‑140 (MQ=255)                                                                                             cgcgAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGa                               >  2:886049/1‑140 (MQ=255)                                                                                             cgcgAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGGGAAATCGAAGAAGTGGa                               >  1:332006/1‑140 (MQ=255)                                                                                             cgcgAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCCGCCCTGGCCGAAGGGGGACTTCAGCAAGTTTGGTGAAAACGAAGAAggggg                               >  2:804687/1‑139 (MQ=255)                                                                                             cgcgAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCCGGCATGCTGCCGCGGGCGAAGGTGGGCTTCGGCAAGTTTGGGGAAATCGAAGAAGGGGa                               >  2:298502/1‑140 (MQ=255)                                                                                                              ttACGTGAAAGAAGCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGGATCCCTGGCATGCTGCCGTGGCCGAAGGGGGACTTCAGCAAGGTTGGTGAAAACGAAGAAGTGGGACTGGGCCGCATCCAGa              >  2:250478/1‑140 (MQ=255)                                                                                                                      aagaagCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATctct      >  1:999695/1‑140 (MQ=255)                                                                                                                      aagaagCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCATGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATctct      >  2:553795/1‑140 (MQ=255)                                                                                                                          agCTATCAAACGCGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTg  <  1:772515/140‑1 (MQ=255)                                                                                                                                      |                                                                                                                                TTCACGCGACTCCGCTGATCCGCCGTCTGGCACGCGAGTTTGGTGTTAACCTTGCGAAAGTGAAGGGCACTGGCCGTAAAGGTCGTATCCTGCGCGAAGACGTTCAGGCTTACGTGAAAGAAGCTATCAAACGTGCAGAAGCAGCTCCGGCAGCGACTGGCGGTGGTATCCCTGGCATGCTGCCGTGGCCGAAGGTGGACTTCAGCAAGTTTGGTGAAATCGAAGAAGTGGAACTGGGCCGCATCCAGAAAATCTCTGGTG  >  NC_000913/126674‑126934

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