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% 185.5 / NA 46S58G (AGC→GGC) rsmE16S rRNA m(3)U1498 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (24/22);  total (24/22)

ACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTG  >  NC_000913/3091273‑3091334
                                |                             
aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                             >  1:1394717/1‑35 (MQ=255)
aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                             <  1:4647603/35‑1 (MQ=255)
 ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                            <  1:425921/35‑1 (MQ=255)
 ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                            >  1:3460402/1‑35 (MQ=38)
    ggTCTTTGACGCCGAAATTACCATCGCCGGCaaaa                         >  1:1792007/1‑35 (MQ=25)
    ggTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                         <  1:1959716/35‑1 (MQ=255)
    ggTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                         >  1:2109485/1‑35 (MQ=255)
     gTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                        <  1:5876381/35‑1 (MQ=38)
     gTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                        <  1:5783615/35‑1 (MQ=255)
     gTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                        <  1:6419198/35‑1 (MQ=255)
     gTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                        <  1:6707296/35‑1 (MQ=38)
      tCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                       <  1:5351723/35‑1 (MQ=37)
        tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                     >  1:7520668/1‑35 (MQ=38)
        tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                     <  1:460624/35‑1 (MQ=255)
        tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                     >  1:2848540/1‑35 (MQ=38)
        tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                     <  1:836163/35‑1 (MQ=255)
        tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                     >  1:1317516/1‑35 (MQ=38)
         ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                    >  1:3123517/1‑35 (MQ=255)
         ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                    >  1:4449326/1‑35 (MQ=255)
         ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                    >  1:1661461/1‑35 (MQ=38)
         ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                    >  1:1661418/1‑35 (MQ=38)
           gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCgg                  >  1:2999618/1‑34 (MQ=37)
            aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                 >  1:1272139/1‑35 (MQ=38)
             cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                <  1:634864/35‑1 (MQ=38)
              gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa               <  1:2713605/35‑1 (MQ=255)
                  aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg           <  1:7713036/35‑1 (MQ=255)
                  aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg           <  1:6046952/35‑1 (MQ=37)
                   aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa          <  1:130046/35‑1 (MQ=255)
                   aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa          >  1:7442342/1‑35 (MQ=255)
                   aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa          >  1:6840575/1‑35 (MQ=37)
                   aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa          >  1:1946348/1‑35 (MQ=37)
                    aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa         <  1:6158955/35‑1 (MQ=255)
                    aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa         >  1:6558760/1‑35 (MQ=255)
                     ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg        >  1:6336212/1‑35 (MQ=37)
                     ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg        >  1:5402623/1‑35 (MQ=255)
                     ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg        <  1:2025213/35‑1 (MQ=37)
                     ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg        >  1:3186904/1‑35 (MQ=255)
                      tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg       >  1:3872775/1‑35 (MQ=37)
                       aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt      >  1:913815/1‑35 (MQ=37)
                       aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt      >  1:5457642/1‑35 (MQ=255)
                        ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg     >  1:4494406/1‑35 (MQ=37)
                        ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg     <  1:6925993/35‑1 (MQ=37)
                           gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg  <  1:7218348/35‑1 (MQ=255)
                           gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg  <  1:3834875/35‑1 (MQ=37)
                           gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg  <  1:5837726/35‑1 (MQ=37)
                           gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg  >  1:5459583/1‑35 (MQ=37)
                           gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg  <  1:5814841/35‑1 (MQ=37)
                                |                             
ACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTG  >  NC_000913/3091273‑3091334

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