Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 229,456 A→G H97R (CAC→CGC)  dkgB → methylglyoxal reductase DkgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913229,4560AG100.0% 99.6 / NA 29H97R (CAC→CGC) dkgBmethylglyoxal reductase DkgB
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/16);  total (13/16)

AAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTG  >  NC_000913/229380‑229532
                                                                            |                                                                            
aaaGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTgg                                                                          >  2:49466/1‑81 (MQ=255)
 aaGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCTATTATGTTGATCTGACGCTAATCCGCTGGc                                                                         >  1:49473/1‑81 (MQ=255)
   gACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCg                                                                       >  1:49469/1‑81 (MQ=255)
   gACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCg                                                                       <  1:49432/81‑1 (MQ=255)
    aCAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGt                                                                      <  1:49433/81‑1 (MQ=255)
       aaTTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCAc                                                                   <  1:49434/81‑1 (MQ=255)
         ttGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCa                                                                 <  1:49435/81‑1 (MQ=255)
           gATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCaaa                                                               >  2:49476/1‑81 (MQ=255)
                 aaGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGa                                                         <  1:49436/81‑1 (MQ=255)
                        aaaGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGtctc                                                   <  2:49438/80‑1 (MQ=255)
                        aaaGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTctct                                                  <  1:49437/81‑1 (MQ=255)
                        aaaGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTctct                                                  <  1:49439/81‑1 (MQ=255)
                           gagagCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGtt                                               <  2:49440/81‑1 (MQ=255)
                                ccTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGaa                                            >  2:49470/1‑79 (MQ=255)
                                              gTACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGCCACCAAACGATGACGTCTCTGTTGAAGAGTTTATGCAGgcgc                            >  1:49467/1‑81 (MQ=255)
                                               tACCGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCgct                           <  1:49441/81‑1 (MQ=255)
                                                 ccGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCgctgc                         >  2:49482/1‑81 (MQ=255)
                                                  cGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCgctgct                        >  1:49464/1‑81 (MQ=255)
                                                  cGATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCgctgct                        >  2:49471/1‑81 (MQ=255)
                                                   gATTATGTTGATCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGctgctg                       >  1:49483/1‑81 (MQ=255)
                                                             aTCTGACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCaaa              >  2:49472/1‑80 (MQ=255)
                                                                 gACGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAAc          >  1:49479/1‑80 (MQ=255)
                                                                  aCGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCaaaaa            <  2:49442/77‑1 (MQ=255)
                                                                  aCGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAAc          <  2:49443/79‑1 (MQ=255)
                                                                  aCGCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACa         <  1:49444/80‑1 (MQ=255)
                                                                    gCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAgg      <  2:49446/81‑1 (MQ=255)
                                                                    gCTAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAgg      <  1:49445/81‑1 (MQ=255)
                                                                      tAATCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGc    >  2:49480/1‑81 (MQ=255)
                                                                        aTCCGCTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTg  <  2:49447/81‑1 (MQ=255)
                                                                            |                                                                            
AAAGACAAATTGATCCCAAGTCTGAAAGAGAGCCTGCAAAAATTGCGTACCGATTATGTTGATCTGACGCTAATCCACTGGCCGTCACCAAACGATGAAGTCTCTGTTGAAGAGTTTATGCAGGCGCTGCTGGAAGCCAAAAAACAAGGGCTG  >  NC_000913/229380‑229532

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