Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,218,815 A→G T74A (ACG→GCG)  aceK → isocitrate dehydrogenase kinase/phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,218,8150AG100.0% 113.1 / NA 32T74A (ACG→GCG) aceKisocitrate dehydrogenase kinase/phosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/19);  total (13/19)

GCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTC  >  NC_000913/4218686‑4218947
                                                                                                                                 |                                                                                                                                    
gcgcAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                             <  1:592083/139‑1 (MQ=255)
   cagcagCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCAt                                                                                                                          >  1:85758/1‑139 (MQ=255)
        gcGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTta                                                                                                                     <  2:358308/139‑1 (MQ=255)
          gTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTCTtact                                                                                                                   <  2:123570/139‑1 (MQ=255)
             tCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACg                                                                                                                <  1:74555/139‑1 (MQ=255)
                    ggCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTaaa                                                                                                         <  1:102541/139‑1 (MQ=255)
                     gCCGACTGGCATGCTGTCCAGCAGTCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAag                                                                                                        <  1:467357/139‑1 (MQ=255)
                     gCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAag                                                                                                        <  1:358673/139‑1 (MQ=255)
                               aTGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                   >  2:435911/1‑118 (MQ=255)
                               aTGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTtact                                                                                                                   <  1:435911/118‑1 (MQ=255)
                                   tGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCgg                                                                                          <  2:138284/139‑1 (MQ=255)
                                               gATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGc                                                                                   <  1:324625/134‑1 (MQ=255)
                                               gATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGc                                                                                   >  2:324625/1‑134 (MQ=255)
                                                   aaaaaCCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTAcc                                                                          <  2:85758/139‑1 (MQ=255)
                                                          gTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCtt                                                                   <  2:826942/139‑1 (MQ=255)
                                                           tATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTc                                                                  >  2:772764/1‑139 (MQ=255)
                                                                          tcatcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGc                                                   <  1:772764/139‑1 (MQ=255)
                                                                           catcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtg                                                                                                              >  2:357745/1‑79 (MQ=255)
                                                                           catcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACgtg                                                                                                              <  1:357745/79‑1 (MQ=255)
                                                                             tcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                             <  1:585511/62‑1 (MQ=255)
                                                                             tcaCGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCgg                                                                                                                             >  2:585511/1‑62 (MQ=255)
                                                                                    ggtctggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACt                                         >  2:441445/1‑139 (MQ=255)
                                                                                        tggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCa                                                                                                                           <  1:241851/53‑1 (MQ=255)
                                                                                        tggtcGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCa                                                                                                                           >  2:241851/1‑53 (MQ=255)
                                                                                                     aCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                        <  2:740579/91‑1 (MQ=255)
                                                                                                     aCTGCGCTGCATTACTAACGGCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCcg                                                                        >  1:740579/1‑91 (MQ=255)
                                                                                                                        ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTa                                                                                                           >  2:533920/1‑37 (MQ=255)
                                                                                                                        ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCgg                                                                                          <  1:504837/54‑1 (MQ=255)
                                                                                                                        ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCgg                                                                                          >  2:504837/1‑54 (MQ=255)
                                                                                                                        ggCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTa     >  1:514185/1‑139 (MQ=255)
                                                                                                                        agCCAAAGTGCGGACGCGGCATTTTTACTACGTGTTa                                                                                                           <  1:533920/36‑1 (MQ=38)
                                                                                                                           cAAAGTGCGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCATGTACTGTCGGTTATTTGACCACCGCTCGCTTACTc  >  2:292778/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                    
GCGCAGCAGCGTTTCGAACAGGCCGACTGGCATGCTGTCCAGCAGGCGATGAAAAACCGTATCCATCTTTACGATCATCACGTTGGTCTGGTCGTGGAGCAACTGCGCTGCATTACTAACGGCCAAAGTACGGACGCGGCATTTTTACTACGTGTTAAAGAGCATTACACCCGGCTGTTGCCGGATTACCCGCGCTTCGAGATTGCGGAGAGCTTTTTTAACTCCGTGTACTGTCGGTTATTTGACCACCGCTCGCTTACTC  >  NC_000913/4218686‑4218947

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