Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,091,305 A→G S58G (AGC→GGC)  rsmE → 16S rRNA m(3)U1498 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,091,3050AG100.0% 101.3 / NA 27S58G (AGC→GGC) rsmE16S rRNA m(3)U1498 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (13/14);  total (13/14)

TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATC  >  NC_000913/3091260‑3091346
                                             |                                         
tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                         <  1:1907372/48‑1 (MQ=255)
  tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaa                                       >  1:693228/1‑48 (MQ=255)
     cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                   <  1:697517/49‑1 (MQ=255)
     cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                   <  1:1315236/49‑1 (MQ=255)
     cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                   <  1:471636/49‑1 (MQ=255)
      ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                  <  1:261578/49‑1 (MQ=255)
             aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                            >  1:1709413/1‑48 (MQ=255)
               cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                         >  1:262144/1‑49 (MQ=255)
                aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                        <  1:717318/49‑1 (MQ=255)
                aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                        >  1:1347792/1‑49 (MQ=255)
                    cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                    >  1:856512/1‑49 (MQ=255)
                     tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                   >  1:318574/1‑49 (MQ=255)
                       tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                  <  1:1678987/48‑1 (MQ=255)
                       tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                  <  1:1075186/48‑1 (MQ=255)
                       tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                 >  1:1940982/1‑49 (MQ=255)
                        gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                >  1:193145/1‑49 (MQ=255)
                          cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg              <  1:1760510/49‑1 (MQ=255)
                            ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg             >  1:366059/1‑48 (MQ=255)
                               aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg         <  1:2040657/49‑1 (MQ=255)
                                aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc        >  1:289014/1‑49 (MQ=255)
                                aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc        >  1:1302103/1‑49 (MQ=255)
                                 aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc       <  1:502513/49‑1 (MQ=255)
                                 aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc       <  1:690697/49‑1 (MQ=255)
                                  ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc        >  1:1266343/1‑47 (MQ=255)
                                   tACCAGCGCCGGCAAAAAAAGCGTGGAAGGGAAGGTGCTGGAAGGc        >  1:1919841/1‑46 (MQ=38)
                                    aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg     <  1:663594/48‑1 (MQ=255)
                                      cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc  <  1:713931/49‑1 (MQ=255)
                                             |                                         
TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATC  >  NC_000913/3091260‑3091346

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