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%	202.5 / NA	46	coding (2069/2664 nt)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base . (29/17);  total (29/17)

GCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATC  >  NC_000913/124953‑125221                                                                                                                                     |                                                                                                                                         gCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCGGa                                                                                                                               >  1:89661/1‑143 (MQ=255)  ccgcaTTCACTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGaa                                                                                                                              >  2:994539/4‑143 (MQ=255)  ccACATTCAGTCGCTGACTATCCCGAACTGGATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGaa                                                                                                                              <  2:667930/143‑1 (MQ=255)        tCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGaaaa                                                                                                                        >  1:88369/1‑143 (MQ=255)         cAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGATTACGATGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAAc                                                                                                                       <  1:706379/143‑1 (MQ=255)            tCGCAGACTCTCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTAc                                                                                                                    <  2:208311/143‑1 (MQ=255)               cTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGGAAAACAAGAGAACGTTAGCGACTACA‑CACTACGCGGCACGAAAACTACcac                                                                                                                 >  2:416391/1‑143 (MQ=255)                    tATCCCGAACTGTCTCTCTTACGACCCGGCTTCCGCTTCCGAAGTTGCTGTCATCATGCATGCCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTCCCACATGcc                                                                                                            >  1:222168/1‑143 (MQ=255)                      tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCAAATGCCgg                                                                                                          >  2:920350/1‑143 (MQ=255)                           gaCTGTATCTGGTGCGGCCCGGCTTACGCTGACGAAGTTGCTGTCATCATGCGTGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATg                                                                                                     <  2:417502/142‑1 (MQ=255)                                   ctctTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAACCTACCACATGCCGGCAATGCCGAAAgg                                                                                             >  1:697509/1‑143 (MQ=255)                                    tctTACGACCCGGCTTACGCTTACGACGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACCGTGAAAAACAAGAGAACGTTCACTACTACA‑CACCACGCTGCACGAACACTACCACACACCGGCCATGCCCGAAGGt                                                                                            >  2:16737/1‑143 (MQ=255)                                    tctTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTAACTACTACA‑CAATACGCTGAACGAAAACTACCACATGCCGGCAATGCAGGACGGt                                                                                            >  2:497453/1‑143 (MQ=255)                                     ctTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTg                                                                                           >  1:998053/1‑143 (MQ=255)                                      ttACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGc                                                                                          >  1:425649/1‑143 (MQ=255)                                             ccGGCTTACGCTGACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                   <  2:547039/143‑1 (MQ=255)                                             ccGGCTTACGCGTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                   <  2:1037322/143‑1 (MQ=255)                                              cGGCTTCCGCTTACGTACGTGCTGTCATCCTGCATGTCTGTCTGGAGCGTATGTACGGTGAAAAACAAGAGATCGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                                  <  2:911511/143‑1 (MQ=255)                                              cGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAg                                                                                  >  1:988938/1‑143 (MQ=255)                                                     cgcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCg                                                                           <  1:994539/143‑1 (MQ=255)                                                      gcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGt                                                                          <  1:1002202/143‑1 (MQ=255)                                                            gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGAGCGGAGGAAGATAGCCGTAAAGGt                                                                    >  2:885312/1‑143 (MQ=255)                                                                ttGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGCATCCGTAAAGGTATCt                                                                >  2:482908/1‑143 (MQ=255)                                                                     gtgagcatgcatgACGGTCTGGAGCGTATGTACGGTGAAACACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACaaa                                                           <  2:786555/138‑1 (MQ=255)                                                                       catcatGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCAGGAAGGTGCTGAGGAAGGTATCCGCAAAGGTATCGACAACCt                                                         >  2:92283/1‑143 (MQ=255)                                                                          catgcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGACGGTATCCGTAAAGGTATCTACAAACTCGa                                                      >  2:140735/1‑143 (MQ=255)                                                                           atgcatgACGGGCAGGAGCGGATGTAGGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGaa                                                     <  1:926308/143‑1 (MQ=255)                                                                           atacatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGCGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCCCATGCCGGCAATGCCGGAAGGTGCTGCGGAAGGTATCCGCAAACGTATCTACTAACTCGaa                                                     >  2:833127/4‑143 (MQ=255)                                                                                 gACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTagt                                               >  2:810252/1‑141 (MQ=255)                                                                                 gACGGCCTGGAGCGTATGTAGGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTGCGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAtt                                               <  2:202034/143‑1 (MQ=255)                                                                                   gtgggtGGGGCGTATGTACCGTTACAAACAAGAGAACGTGTACTACTACA‑CACGACGCTGAACGAAAACTACAACAGGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGa                                             <  1:174200/138‑1 (MQ=255)                                                                                                  gTACGGTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaa                              >  1:811529/1‑143 (MQ=255)                                                                                                   tACGGTGAAACCCGAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACGATTGAAGGTAGCAAAGGTaaa                             <  2:300623/143‑1 (MQ=255)                                                                                                   gaCGGTGAAAAACAAGAGAACGTTTACTACTACA‑CGCTGCGCTGGACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGGAAAGGTACCGACAAACTCGAAACTCTTTAAGGTAGCAAAGGTaaa                             <  2:176238/142‑1 (MQ=255)                                                                                                       gTGAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAACGGTACAGTTc                         >  1:69829/1‑143 (MQ=255)                                                                                                         gAAAAGCAAGAGAACGTGTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAg                       <  2:52433/143‑1 (MQ=255)                                                                                                         gAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACTGAGTGACGACAACGCCGGCAGTTGCGTATGAAGCTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAGGGTAAAGTTCAg                       >  2:219325/1‑143 (MQ=255)                                                                                                         gAAAAACAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATGGAGGGTAGAAAAGGTAAAGCTCAg                       >  2:243197/1‑143 (MQ=255)                                                                                                          aaaaaCAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgc                      >  1:813656/1‑143 (MQ=255)                                                                                                             aaCAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCAAATGCCGGCAATGCCGGAAGGTGCTGGGGAAGGTATCCGTAAAGGTATCTACAAACTAGAAAGGATTGAAGGTAGCGAGGGAAAAGTTAAgctgc                   >  1:408486/1‑143 (MQ=255)                                                                                                               cAAGAGAACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTc                 >  2:944324/1‑143 (MQ=255)                                                                                                                  gagaACGTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGc              >  1:329128/1‑143 (MQ=255)                                                                                                                     aaCGTTTACTACTACA‑CACTACGCTGAACGAAAGCTACCCCATGCCGACAATACCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTGCAAACACGAAACTATGGAAGGGAGCAACGGTGAAGTCAAGCAGCTAGGcgac           >  2:481845/1‑140 (MQ=255)                                                                                                                        gTTTACTACTACA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGt        >  2:445638/1‑143 (MQ=255)                                                                                                                         tttACTACTACA‑CACTACGCTGAACGAAAACTACCCCATGCCTGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGtt       <  1:87533/143‑1 (MQ=255)                                                                                                                              tactacA‑CACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATc  >  1:936160/1‑143 (MQ=255)                                                                                                                                     |                                                                                                                                         GCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATC  >  NC_000913/124953‑125221

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