Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	125,085	T→A	I690N (ATC→AAC)	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	A	100.0%	160.3 / NA	52	I690N (ATC→AAC)	aceE	pyruvate dehydrogenase E1 component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (28/24);  total (28/24)

CACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCT  >  NC_000913/124949‑125223                                                                                                                                         |                                                                                                                                           cacaGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTAc                                                                                                                                      <  1:540170/143‑1 (MQ=255)   caGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGc                                                                                                                                    <  1:563609/143‑1 (MQ=255)       cacaTGCAGTCGCTGACTATCCCGCACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGaa                                                                                                                                <  1:685747/143‑1 (MQ=255)        acaTTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCCTGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGACCGTTTACTACTACAACACTACGCTGAAc                                                                                                                               >  2:244670/1‑143 (MQ=255)          aTTCAGTCGCTGACTATCCCGAACTGTCTCTCTTCCGACCCGGCTTACGCTTAAGAAGTTACTGTCAACATGCATGACAGACTAGAACGTACGAACGATGAAAAACAATAGAACGTACAATCCTACAACACTACGCAGAACGa                                                                                                                             >  2:263176/1‑143 (MQ=255)              aGTCGCAGACTATCCCAAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATTCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGCGAACCGTTACTACTACAACACTAAGCTGCACGaaaaa                                                                                                                         >  2:265332/1‑142 (MQ=255)                   cTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACca                                                                                                                    <  1:186040/143‑1 (MQ=255)                      aCTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACAt                                                                                                                 >  2:827043/1‑143 (MQ=255)                          tCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCg                                                                                                             >  2:226950/1‑143 (MQ=255)                                  tGTATCTCTTACGACCCTGCTTACGCTTAAACAGTTTCCGTCACCATGCATGACGGTCTGGAGCGTATGTACGCATAAAAACAAGAGAACGTTGACTACTACAACACGACGCTGAACGAAAACAACCACATGCCGGCAATGca                                                                                                     >  2:420768/1‑142 (MQ=255)                                   gTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGTTGAACGAAAACTACCACATGCCGGCAATGACg                                                                                                    >  1:500113/1‑143 (MQ=255)                                    tATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCgg                                                                                                   >  2:540110/1‑143 (MQ=255)                                      tctcTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGaa                                                                                                 <  1:740993/143‑1 (MQ=255)                                      tctcTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTGTACTACTACAACACTACGCTGAACGAAAACTACCACACGCCGGCAAAGCCggga                                                                                                 >  2:551584/1‑141 (MQ=255)                                       ctctTACGACCCGGCTTACGCTTACGTAGTTGTTGTCATCATGCATGACGGACTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAg                                                                                                <  1:641922/143‑1 (MQ=255)                                            aCGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCGGGCACTGCCGGAAGTTGCt                                                                                           >  1:696884/1‑143 (MQ=255)                                            aCGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCt                                                                                           <  2:230343/143‑1 (MQ=255)                                                 ccGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTAGTACAACAATACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGCAGGTGCAGACGa                                                                                      >  2:423469/1‑143 (MQ=255)                                                    gcttacgcttacGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAgg                                                                                   <  1:827043/143‑1 (MQ=255)                                                              acgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa                                                                         >  1:2416/1‑143 (MQ=255)                                                              acgAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGGGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa                                                                         <  2:88476/143‑1 (MQ=255)                                                              acgAAGTTGCTGTCATCACGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa                                                                         <  1:202090/143‑1 (MQ=255)                                                              acgAAGTGGCTGTCGTCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTaaa                                                                         <  2:437312/143‑1 (MQ=255)                                                                gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAgg                                                                       >  1:686772/1‑143 (MQ=255)                                                                gAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACACTGAACGAAAACTACCACATGCCGGCAATGCCGGAAAGTGCTGAGGAAGGTATCCGTAAAgg                                                                       >  1:435446/1‑143 (MQ=255)                                                                   gctGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACCACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTAt                                                                    <  2:89740/141‑1 (MQ=255)                                                                        tggcatcatGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACa                                                               <  1:823747/140‑1 (MQ=255)                                                                        tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                     <  1:414256/121‑1 (MQ=255)                                                                        tGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGaa                                                                                     >  2:414256/1‑121 (MQ=255)                                                                                    tgACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTa                                                   <  1:540110/143‑1 (MQ=255)                                                                                     gACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTCCCAAATGCCGGCAACGCCGGAAGGTGCTAAGGAAGGTATCCGTAAAGGTATCTACAAACTCCAAACTAt                                                  >  2:531197/1‑143 (MQ=255)                                                                                     gACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAt                                                  <  1:411498/143‑1 (MQ=255)                                                                                      aCGGTCTGGAGCTTATGTACGGTTAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAtt                                                 <  1:442174/143‑1 (MQ=255)                                                                                      aCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAtt                                                 >  2:720220/1‑143 (MQ=255)                                                                                      aCGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTAtt                                                 <  1:141751/143‑1 (MQ=255)                                                                                            tGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGt                                           >  1:585175/1‑143 (MQ=255)                                                                                             ggAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTa                                          <  1:508092/143‑1 (MQ=255)                                                                                                gCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCa                                       >  1:547471/1‑143 (MQ=255)                                                                                                gCGTATGTACGCTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCCGAACGAAAACTACCACATGCCGGCAATGCCCGAAGGTGCTGAGGAAGGTCTCCGTCCAGGCACCTCCAAACTCGAAACTATCGCAGgcagca                                       >  2:604852/1‑143 (MQ=255)                                                                                                  gTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCaaa                                     <  1:576746/143‑1 (MQ=255)                                                                                                     tGTACGGTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGACGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGt                                  >  1:625382/1‑143 (MQ=255)                                                                                                        aCGGGGAAACACAAGAGAACGTTTACTGCGACAACACGACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTaaa                               <  2:773451/143‑1 (MQ=255)                                                                                                           gTGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGAGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGtt                            >  1:409351/1‑143 (MQ=255)                                                                                                            tGAAAAACAAGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTc                           >  1:684672/1‑143 (MQ=255)                                                                                                               aaaaCAAGAGAACGTGTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgc                        <  1:609780/143‑1 (MQ=255)                                                                                                                    aaGAGAACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAgctgccc                   >  1:35623/1‑141 (MQ=255)                                                                                                                      gagaACGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCgg                 >  1:576760/1‑143 (MQ=255)                                                                                                                      gagaACGTTCACTCCTCCAACACTTCGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCCTAAAGGTATCTACAAACTCGAAACTATTGAATGTAGCAAAGGTAAAGCTCAGCTGCTCgg                 <  1:226950/143‑1 (MQ=255)                                                                                                                           cGTTTACTACTACAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCg            >  2:596635/1‑143 (MQ=255)                                                                                                                           cGTTTACTACTACAACACTAAGCGGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCg            >  2:755301/1‑143 (MQ=255)                                                                                                                                actactacAACACCACGCTGAACGAAAACGCCCACATGACGGCAATGCCGGAAGGTGCTGAGCAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCt       <  1:463754/143‑1 (MQ=255)                                                                                                                                     tacAACACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGACGGTAGCAAAGGTAAAGTTCCGCTGCTCGGCTCCGGTCCTATCCt  >  1:313700/1‑143 (MQ=255)                                                                                                                                         |                                                                                                                                           CACAGCCACATTCAGTCGCTGACTATCCCGAACTGTATCTCTTACGACCCGGCTTACGCTTACGAAGTTGCTGTCATCATGCATGACGGTCTGGAGCGTATGTACGGTGAAAAACAAGAGAACGTTTACTACTACATCACTACGCTGAACGAAAACTACCACATGCCGGCAATGCCGGAAGGTGCTGAGGAAGGTATCCGTAAAGGTATCTACAAACTCGAAACTATTGAAGGTAGCAAAGGTAAAGTTCAGCTGCTCGGCTCCGGTTCTATCCT  >  NC_000913/124949‑125223

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