Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	125,085	Δ1 bp	coding (2069/2664 nt)	aceE →	pyruvate dehydrogenase E1 component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	125,085	0	T	.	100.0%	156.9 / NA	36	coding (2069/2664 nt)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (16/20);  total (16/20)

GTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGG  >  NC_000913/124947‑125214                                                                                                                                           |                                                                                                                                  gTCACAGCCACGTTCAGACGCTGACTATCCCGACCTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTa                                                                                                                              <  2:2323/143‑1 (MQ=255)   cacaGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTATCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACg                                                                                                                            <  1:61469/143‑1 (MQ=255)        ccACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGaa                                                                                                                       <  1:586254/143‑1 (MQ=255)              tCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCGTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGaaaa                                                                                                                 <  1:319720/143‑1 (MQ=255)                    gCTGACTATCCCGAACCGTACCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACca                                                                                                           <  1:250226/143‑1 (MQ=255)                         cTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCACGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGc                                                                                                      <  1:55717/143‑1 (MQ=255)                              ccGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTGCGGTGAAAAACAAGAGAACGTTTACTACTACA‑CAGTACGCTGAACGAAAACTACCACAAGCCAGCa                                                                                                 >  1:202347/1‑143 (MQ=255)                              ccGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCAGAACGAAAACTACCCCATGCCGGCa                                                                                                 >  2:467614/1‑143 (MQ=255)                                  aCTGTTTCTCTTACGATCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGTGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTCCGCTGAACGAAAACTACCACATGCCGGCAATGc                                                                                             <  2:50103/143‑1 (MQ=255)                                      tATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGa                                                                                         >  2:47974/1‑143 (MQ=255)                                       aTCTCTTACGGCCCGGCTTACTCTTACGTAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACTGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGaa                                                                                        <  2:58873/143‑1 (MQ=255)                                              aCGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTg                                                                                 >  1:346611/1‑143 (MQ=255)                                                            gcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGCCAGGAGCTGCGGAAGGTAGCCGt                                                                   >  2:495527/1‑143 (MQ=255)                                                              ttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGCAGGTATCCGtca                                                                 >  2:113526/1‑141 (MQ=255)                                                                  gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGt                                                             <  1:380961/143‑1 (MQ=255)                                                                    aggtggtGTCGTCATGCGTGGCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTAt                                                           <  2:177106/137‑1 (MQ=255)                                                                      gtGCGGTGATCATGCATGACGGGCTGGAGCGTACGGACGGTGAGAAACAAGAGAACGTTTCCTACTGCA‑CGCTCCGCGGAACGAAAACGACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCt                                                         <  2:618210/142‑1 (MQ=255)                                                                          tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACac                                                     >  1:124722/1‑142 (MQ=255)                                                                          tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaa                                                     >  2:559519/1‑143 (MQ=255)                                                                          tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaa                                                     >  1:152415/1‑143 (MQ=255)                                                                           gTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa                                                    >  1:431684/1‑143 (MQ=255)                                                                           gTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa                                                    <  1:205541/143‑1 (MQ=255)                                                                           gTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTCTACTCCTACA‑CACTACGCTGAACGAAAACTACCACATGACGACAAAACCGGAAAGTACAGAAAAAGGCAGACTGAAAGGTATCCACaaa                                                    >  1:523793/1‑143 (MQ=255)                                                                            tcatcaTGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAAc                                                   >  1:139032/1‑143 (MQ=255)                                                                               tcatGCCTGACGGTCCGGAGCGTAGGTACGGGGAAAGACACGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCCCATGCCGGCAATGCCGGAAGGTGCTGAGGAATGTATCCGTAAAGGTATCTACAAACTCg                                                <  1:59139/143‑1 (MQ=255)                                                                               tcatGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGACAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCg                                                <  2:433382/143‑1 (MQ=255)                                                                                catgcaggACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGa                                               <  1:507760/143‑1 (MQ=255)                                                                                   gcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCCGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGTTATCTACAAACTAGAAAc                                            >  1:277278/1‑143 (MQ=255)                                                                                          aatcaagAGCGTATGCACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGaa                                     <  1:353075/137‑1 (MQ=255)                                                                                               ggAGCGTATGTACGGTGAAAACCAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAg                                <  2:332453/143‑1 (MQ=255)                                                                                                   tgTATGTACAGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaaa                            <  1:463537/142‑1 (MQ=255)                                                                                                     tATGGACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAgg                          <  2:517719/143‑1 (MQ=255)                                                                                                              tGAAAAACAAGAGAACGTTTACTACTACA‑CTCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCa                 <  2:37643/143‑1 (MQ=255)                                                                                                                     cAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATCGAAGGTAGCCAAGGTAAAGTTCAGCTGATc          >  2:596597/1‑143 (MQ=255)                                                                                                                            aCGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGCAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCg   >  1:277691/1‑143 (MQ=255)                                                                                                                             cGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCgg  >  1:408712/1‑143 (MQ=255)                                                                                                                                           |                                                                                                                                  GTCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGG  >  NC_000913/124947‑125214

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