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

Only 100 of 134 total aligned reads displayed.

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑3091351
                                              |                                              
gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              <  2:2689437/49‑1 (MQ=255)
 tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                             <  1:2304133/49‑1 (MQ=255)
   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           >  1:1535734/1‑49 (MQ=255)
   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           >  1:2166964/1‑49 (MQ=255)
   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           <  2:4605235/49‑1 (MQ=255)
   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           <  1:2419723/49‑1 (MQ=255)
    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  2:4628696/1‑49 (MQ=255)
    gACGGTAGCAACCAGGACTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  1:3475076/1‑49 (MQ=38)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                        <  2:1564320/50‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:3725191/1‑49 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         <  1:1117066/49‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  2:3527966/1‑49 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:246025/1‑49 (MQ=255)
      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                       <  1:2696037/50‑1 (MQ=255)
      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  2:1042648/1‑48 (MQ=255)
        gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                     >  2:483648/1‑50 (MQ=255)
         tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    >  2:2243516/1‑50 (MQ=255)
         tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    >  1:5436917/1‑50 (MQ=255)
          aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   >  1:112417/1‑50 (MQ=255)
           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  <  2:3532255/50‑1 (MQ=255)
            cAACCATGTCTTTGACGTCGAAATTACCAGCACCGGAAAAAAAAGCGTgg                                 <  2:5310574/50‑1 (MQ=25)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 <  2:5310615/50‑1 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  2:1841967/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:3160047/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  2:1558355/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  2:3106080/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  2:3520424/1‑49 (MQ=255)
              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:3189818/1‑48 (MQ=255)
              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:4052469/50‑1 (MQ=255)
              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:3181150/50‑1 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:4496239/50‑1 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  2:312370/50‑1 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  2:4328727/50‑1 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:1229113/50‑1 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  2:3159074/1‑50 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:225086/1‑50 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             <  2:2592012/50‑1 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:4046969/1‑50 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  2:5429360/1‑49 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            <  1:3838624/50‑1 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:3071560/1‑49 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  2:2191972/1‑49 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  1:4023139/50‑1 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  2:1167703/50‑1 (MQ=255)
                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:5250553/1‑48 (MQ=255)
                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa                          <  2:4850328/50‑1 (MQ=255)
                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         >  2:2110354/1‑50 (MQ=255)
                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         >  1:588520/1‑50 (MQ=255)
                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        >  2:5663107/1‑50 (MQ=255)
                      tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  2:1611438/1‑50 (MQ=255)
                       ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      >  2:2743682/1‑50 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  2:3789536/48‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     >  1:855035/1‑50 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     >  1:2781141/1‑50 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  2:965976/49‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  2:3253593/49‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  1:1759080/49‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  2:5455216/49‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   >  2:4471175/1‑50 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   >  1:3952161/1‑50 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   <  2:3228408/50‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   >  1:1634627/1‑50 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   >  2:638028/1‑50 (MQ=255)
                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  <  2:3140827/50‑1 (MQ=255)
                           cGCCGAAAGTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  <  2:3271626/50‑1 (MQ=38)
                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  >  2:2500014/1‑49 (MQ=255)
                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  >  1:3140541/1‑49 (MQ=255)
                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa                 <  2:3075147/50‑1 (MQ=255)
                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               >  1:3315698/1‑50 (MQ=255)
                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               <  2:4794228/50‑1 (MQ=255)
                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              <  2:930609/50‑1 (MQ=255)
                               gAAATTACCAGCACCGGCAAAAAAAGCGTGGAAGTGCAGGTGCTGGAAgg              >  2:2688594/1‑50 (MQ=38)
                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             <  1:298835/50‑1 (MQ=255)
                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5259066/1‑49 (MQ=255)
                                  aTTCCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:4105570/50‑1 (MQ=38)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5253131/1‑48 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:1561933/50‑1 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  2:370857/50‑1 (MQ=255)
                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  2:789784/50‑1 (MQ=255)
                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  2:976065/50‑1 (MQ=255)
                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  1:3837496/50‑1 (MQ=255)
                                    tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:3903967/1‑46 (MQ=255)
                                    tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:5469918/1‑46 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5471360/1‑44 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:28392/1‑44 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:3142690/1‑44 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:4946168/1‑44 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  2:2849063/49‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  2:2624398/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  2:1866003/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  2:555345/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  2:3167909/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       >  2:5581627/1‑49 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       >  1:5086944/1‑49 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     >  1:4766609/1‑50 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     <  2:4978904/50‑1 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     <  2:4716111/50‑1 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     >  2:4267719/1‑50 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     >  1:343236/1‑50 (MQ=255)
                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  <  1:4260236/50‑1 (MQ=255)
                                              |                                              
GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑3091351

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