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

TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091260‑3091351
                                             |                                              
tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                             >  1:3132706/1‑49 (MQ=255)
tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                             >  1:4518390/1‑49 (MQ=255)
tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                             >  1:4650385/1‑49 (MQ=255)
tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                             >  1:5910035/1‑49 (MQ=255)
  tGACGGTAGCAACCAGGTCTTTGACGCCTAAATTACCAGCGCCGGCaaa                                           >  1:5835883/1‑49 (MQ=39)
  tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           >  1:4081472/1‑49 (MQ=255)
   gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACGAGCGCCGGCaaaa                                          <  1:1756018/49‑1 (MQ=39)
    aCGGTAGCCCCCATGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                        <  1:2129770/50‑1 (MQ=37)
     cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:2904658/1‑48 (MQ=255)
     cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:1354295/1‑48 (MQ=255)
      ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:6067341/1‑50 (MQ=255)
      ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      <  1:131581/50‑1 (MQ=255)
      ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:3051137/1‑50 (MQ=255)
       gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                     >  1:6099633/1‑50 (MQ=255)
       gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                     >  1:6100531/1‑50 (MQ=255)
        tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    <  1:3410698/50‑1 (MQ=255)
         aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   >  1:566253/1‑50 (MQ=255)
          gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   <  1:1225883/49‑1 (MQ=255)
           cAACCAGGACTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:2405819/1‑49 (MQ=38)
            aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:413146/1‑49 (MQ=255)
            aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:915703/1‑49 (MQ=255)
            aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  1:2967769/50‑1 (MQ=255)
            aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  1:4608208/50‑1 (MQ=255)
             aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:4384197/1‑48 (MQ=255)
             aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:3790866/50‑1 (MQ=255)
             aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:2368209/50‑1 (MQ=255)
             aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:5818777/50‑1 (MQ=255)
             aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:4917651/50‑1 (MQ=255)
              ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:3683044/50‑1 (MQ=255)
              ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              >  1:4538753/1‑50 (MQ=255)
               cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:2805014/1‑50 (MQ=255)
               cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             <  1:6026877/50‑1 (MQ=255)
               cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             <  1:2477648/50‑1 (MQ=255)
               cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:2185302/1‑50 (MQ=255)
               cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:2462297/1‑50 (MQ=255)
                aGGTATTTGGCGACGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            <  1:3133152/50‑1 (MQ=37)
                 ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:1792525/1‑49 (MQ=255)
                 ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:3402068/1‑49 (MQ=255)
                 ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  1:331648/50‑1 (MQ=255)
                   tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         >  1:4081634/1‑50 (MQ=255)
                     tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:4458387/1‑50 (MQ=255)
                      ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      >  1:2202482/1‑50 (MQ=255)
                      ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      >  1:4607965/1‑50 (MQ=255)
                       tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        <  1:4156857/47‑1 (MQ=255)
                       tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  1:1451550/48‑1 (MQ=255)
                       tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  1:3802916/49‑1 (MQ=255)
                         aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  1:6317317/49‑1 (MQ=255)
                         aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  1:2437905/49‑1 (MQ=255)
                          cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  <  1:4753294/50‑1 (MQ=255)
                           gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGTAGTGAAGGTGCTgg                  >  1:2880916/1‑49 (MQ=38)
                           gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa                 <  1:4950644/50‑1 (MQ=255)
                             cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               <  1:6425461/50‑1 (MQ=255)
                              gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              <  1:135182/50‑1 (MQ=255)
                              gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              >  1:5132432/1‑50 (MQ=255)
                               aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             <  1:6412170/50‑1 (MQ=255)
                               aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5655659/1‑50 (MQ=255)
                                aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:3669640/50‑1 (MQ=255)
                                aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5569031/1‑49 (MQ=255)
                                 aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:651029/1‑48 (MQ=255)
                                 aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:6066466/1‑48 (MQ=255)
                                 aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:2541607/1‑48 (MQ=255)
                                  ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:3659168/1‑47 (MQ=255)
                                  ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  1:1438593/50‑1 (MQ=255)
                                  ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  1:2784662/50‑1 (MQ=255)
                                   tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:2284171/1‑46 (MQ=255)
                                   tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:1853982/1‑46 (MQ=255)
                                    aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  1:5727096/49‑1 (MQ=255)
                                    aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  1:1884409/49‑1 (MQ=255)
                                     ccAGCGCCGGCAAATAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  1:4440512/49‑1 (MQ=39)
                                     ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  1:4803658/49‑1 (MQ=255)
                                     ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:1455868/50‑1 (MQ=255)
                                     ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:2418728/50‑1 (MQ=255)
                                      cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  1:1703793/50‑1 (MQ=255)
                                      cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  1:45844/50‑1 (MQ=255)
                                      cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  1:2665334/50‑1 (MQ=255)
                                      cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       >  1:4761906/1‑49 (MQ=255)
                                        gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    <  1:1548414/50‑1 (MQ=255)
                                         cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacg   <  1:2600883/50‑1 (MQ=255)
                                         cgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacg   <  1:2715438/50‑1 (MQ=255)
                                          gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  >  1:1740809/1‑50 (MQ=255)
                                             |                                              
TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091260‑3091351

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