Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	125,085	2 bp→CA	coding (2069‑2070/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	C	100.0%	141.9 / NA	48	I690T (ATC→ACA)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/24);  total (24/24)
*	NC_000913	125,086	0	C	A	100.0%	133.8 / NA	48	I690T (ATC→ACA)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (24/24);  total (24/24)

TCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTAT  >  NC_000913/124948‑125220                                                                                                                                          ||                                                                                                                                       tCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTa                                                                                                                                    >  1:295532/1‑143 (MQ=255) tCACAGCCACATTCAGACGCTGACTATCCCGAACCGTATCGCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCCGGAGCGTATGTACGGTGAAAACCAAGACAACGTTTACTCCTGCACAACTa                                                                                                                                    >  2:398893/1‑143 (MQ=255)                 tCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTa                                                                                                                    >  2:122022/1‑143 (MQ=255)                   gCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGACGTCGCTGTCCTCATGCATCACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTATACTAATACACAACTACGCTGAACCAAAACTAcc                                                                                                                  >  1:159575/1‑143 (MQ=255)                   cctggccaTCCCGGACTGCATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGCCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACGACGCTGAACGAAAACTAcc                                                                                                                  <  2:240231/136‑1 (MQ=255)                             ccGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCCTCATGCATGACGTTCTGGGGCGTATGTACGGGGAAAACCACGAGCACGCTTACTACTACACAACTACGCTGAACCAAAACTCCCCCATGCCCGc                                                                                                        >  2:170097/1‑143 (MQ=255)                              cGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTGGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCa                                                                                                       >  2:634628/1‑143 (MQ=255)                                gaCTGTATCTCTTACGACCCGGCTTACGCTTACGAAGGTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAAt                                                                                                     <  2:641557/142‑1 (MQ=255)                                    gTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCg                                                                                                 <  1:361671/143‑1 (MQ=255)                                     tATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGCATGTACGGTGAAAAACAAGAGCACGTTTACTACTACACAACTACGCTGAACGCAAACTACAAGATGCTACAAATGCCgg                                                                                                >  2:135929/1‑143 (MQ=255)                                       tctcTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTCCGCTGAACGAAAACTACCACATGCCAGCAAAGCCGGaa                                                                                              >  1:697837/1‑143 (MQ=255)                                        ctctTCCGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                                                             <  1:629429/143‑1 (MQ=255)                                         gcctACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTCCCACATGCCGGCAATGCCGGAAgg                                                                                            <  1:217962/140‑1 (MQ=255)                                                   cGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                  <  1:698248/143‑1 (MQ=255)                                                   cGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                  <  2:668874/143‑1 (MQ=255)                                                          cgcttacgAAGTTGCTGTCATCATGCATGACGGTTTGGATCATATGTACGGTGAAAAACAATAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATcc                                                                           <  2:390219/143‑1 (MQ=255)                                                           gcttacgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGATGCTGAGGAAGGTAGCCg                                                                          >  2:297975/1‑143 (MQ=255)                                                                cgAACTTGCTGTCATCATGGATGACGGTCTGGTGCGTATGTACGGTGAAAAACAATAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCTCATGTCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAg                                                                     <  2:414974/143‑1 (MQ=255)                                                                        atGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAGCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTAc                                                             <  2:798813/142‑1 (MQ=255)                                                                              ttatggatgACTGTATGTAGCGTATGTACAGTGAAAAAAAAGAGAACGAAAACGACAACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTc                                                       <  1:35618/141‑1 (MQ=255)                                                                                atgcatgACGGTCTGGAGCGTATGTCCGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGa                                                     <  1:859242/143‑1 (MQ=255)                                                                                  gcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGaaa                                                   <  2:444419/143‑1 (MQ=255)                                                                                  gcatgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAAATACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAAATCGaaa                                                   >  2:619035/1‑143 (MQ=255)                                                                                  gcatgACGGTCTGCAGCGTATGTACGGTGAAAAACAAGAGAACGTTGACTACTACACAACTACGCTGAATGAAAAATACGACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGaaa                                                   <  1:559538/143‑1 (MQ=255)                                                                                   catgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAAc                                                  <  1:249087/143‑1 (MQ=255)                                                                                     tgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTa                                                <  1:230245/143‑1 (MQ=255)                                                                                     tgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTa                                                <  2:3393/143‑1 (MQ=255)                                                                                        cGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTg                                             >  1:736493/1‑143 (MQ=255)                                                                                          gTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTGCGACATGCCGGCAATGCCGCAAGGTTCCGATGGAGGTATCCGTACAGGTATCTACAAACTCGAAACTATTGaa                                           >  1:454129/1‑143 (MQ=255)                                                                                             tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAgg                                         <  2:373313/142‑1 (MQ=255)                                                                                             tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAgg                                         >  1:373313/1‑142 (MQ=255)                                                                                               gAGCGTATGTACGGGGAAAAACAAGAGAACGTTTACTACTACACAGCTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAg                                      <  2:17349/143‑1 (MQ=255)                                                                                                 gccgacgTACGGGCAACAACAAGAGAACGTTTACTACTACACAACCACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCa                                    <  1:539895/137‑1 (MQ=255)                                                                                                   gTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaaa                                  <  1:719550/143‑1 (MQ=255)                                                                                                       gTACGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTa                              >  1:830340/1‑143 (MQ=255)                                                                                                         aCGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                            >  1:823256/1‑143 (MQ=255)                                                                                                          cGGTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAg                           <  1:122022/143‑1 (MQ=255)                                                                                                           ggTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACCAACTCGAAACTATTGTAGGTAGCAAAGGTAAAGt                          >  1:295182/1‑143 (MQ=255)                                                                                                           ggTGAAAAACAAGAGAACGTTAAATACTACACAACGACGATGAACGAAAACTACCAAATGCCGGCAATGACGGAAGGTGCTGAGCAAGGTAACAGAAAAGGTATCTACAAAATCGAAACTATTGAAGATAACAAAGGAAACGt                          >  2:79109/1‑143 (MQ=255)                                                                                                            gTGAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGtt                         <  1:360738/143‑1 (MQ=255)                                                                                                              gAAAAACAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTCAAGTTCa                       >  1:267721/1‑143 (MQ=255)                                                                                                                 aaaCAAGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgct                    >  2:880412/1‑143 (MQ=255)                                                                                                                  accAAGAGAACCGTTAGTGCTACCCACCTACGCTGGAGGAAAACTACCAGATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGCATCTACATACTCGATACGATTGAAGGTAGCAAAGGTAAAGTTCAgctg                   <  2:487128/141‑1 (MQ=255)                                                                                                                     aaGAGAACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTc                >  2:669968/1‑143 (MQ=255)                                                                                                                       gagaACGTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTACAGTTCAGCTGCTCgg              >  2:14715/1‑143 (MQ=255)                                                                                                                             gTTTACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTCTTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCgc        >  1:665912/1‑142 (MQ=255)                                                                                                                              tttACTACTACACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAACGGTAAAGTTCACCTGCTCGGCTACGGt       >  2:458722/1‑143 (MQ=255)                                                                                                                                   tactacACAACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTCAAGGTATCTACAAACTCGAAACTATTGACGGTAGCAAAGGTAAAGTTCAGCTCCGCGGCTCCGGTGCTac  >  2:661312/1‑142 (MQ=255)                                                                                                                                          ||                                                                                                                                       TCACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTAT  >  NC_000913/124948‑125220

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