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

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑3091351
                                              |                                              
gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              >  1:2972495/1‑49 (MQ=255)
 tttGACGGTAGCAACCCGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                             >  1:4947034/1‑49 (MQ=38)
   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           <  1:2911185/49‑1 (MQ=255)
    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          <  1:4051454/49‑1 (MQ=255)
    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  1:5750009/1‑49 (MQ=255)
    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          <  1:2279735/49‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                        <  1:1777350/50‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         <  1:3645375/49‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         <  1:1135680/49‑1 (MQ=255)
      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                       <  1:5015189/50‑1 (MQ=255)
      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                       <  1:410113/50‑1 (MQ=255)
      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:1006814/1‑48 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      <  1:5668300/50‑1 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      <  1:33483/50‑1 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:4293557/1‑50 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      <  1:6642212/50‑1 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:6197956/1‑50 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      <  1:1066201/50‑1 (MQ=255)
        gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                     >  1:5418951/1‑50 (MQ=255)
        gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                     <  1:3930147/50‑1 (MQ=255)
         tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    <  1:7250337/50‑1 (MQ=255)
          aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   >  1:2917913/1‑50 (MQ=255)
           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   <  1:7458864/49‑1 (MQ=255)
           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   <  1:161540/49‑1 (MQ=255)
           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  <  1:6540931/50‑1 (MQ=255)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 <  1:6628291/50‑1 (MQ=255)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:2731540/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:3592228/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:3599315/1‑49 (MQ=255)
              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:4292211/1‑48 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              >  1:1067343/1‑50 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:5471254/50‑1 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:3059011/50‑1 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:605636/1‑50 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:5368297/1‑50 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:6573140/1‑50 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:7210291/1‑49 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:3920435/1‑49 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:315469/1‑49 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            <  1:2108542/50‑1 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  1:5935721/50‑1 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  1:3608256/50‑1 (MQ=255)
                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa                          <  1:7480227/50‑1 (MQ=255)
                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         <  1:2433196/50‑1 (MQ=255)
                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         <  1:3446181/50‑1 (MQ=255)
                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         >  1:3484867/1‑50 (MQ=255)
                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        >  1:7121193/1‑50 (MQ=255)
                      tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:1332512/1‑50 (MQ=255)
                       ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      >  1:7055436/1‑50 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        <  1:4533322/47‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        <  1:4936345/47‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  1:5330781/49‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  1:1053210/49‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  1:4466668/49‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     >  1:7157518/1‑50 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     >  1:7468249/1‑50 (MQ=255)
                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     <  1:633380/49‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  1:549204/49‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   <  1:7224009/50‑1 (MQ=255)
                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  <  1:808010/50‑1 (MQ=255)
                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  >  1:7438064/1‑49 (MQ=255)
                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa                 <  1:6241004/50‑1 (MQ=255)
                             ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  >  1:5989012/1‑48 (MQ=255)
                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              >  1:4421150/1‑50 (MQ=255)
                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              <  1:6779619/50‑1 (MQ=255)
                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5487669/1‑50 (MQ=255)
                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:1803046/1‑50 (MQ=255)
                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:4843502/50‑1 (MQ=255)
                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:7322901/50‑1 (MQ=255)
                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:1316215/1‑49 (MQ=255)
                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5239973/1‑49 (MQ=255)
                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:7371235/1‑49 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:4894670/1‑48 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:6003434/1‑48 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:6665167/50‑1 (MQ=255)
                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:3836545/1‑47 (MQ=255)
                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:7176762/1‑45 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:4668288/1‑44 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:4055875/1‑44 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  1:4577230/49‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:6184952/50‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:303445/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  1:474247/50‑1 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     >  1:4149898/1‑50 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     >  1:3513419/1‑50 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     >  1:2637367/1‑50 (MQ=255)
                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    >  1:5026077/1‑50 (MQ=255)
                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    <  1:6838956/50‑1 (MQ=255)
                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  >  1:5100145/1‑50 (MQ=255)
                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  <  1:4883260/50‑1 (MQ=255)
                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  <  1:2156763/50‑1 (MQ=255)
                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  >  1:5551883/1‑50 (MQ=255)
                                              |                                              
GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑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: ‑