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

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGAC  >  NC_000913/3091259‑3091349
                                              |                                            
gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                            >  1:2851662/1‑49 (MQ=255)
gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                            <  1:1972736/49‑1 (MQ=255)
 tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                           >  1:1154020/1‑49 (MQ=255)
   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                         >  2:1158830/1‑49 (MQ=255)
    gACGGTAGCAAACAGGTCTTTGACGCCGAGATTACCAGCGCCGGCaaaa                                        >  1:3534876/1‑49 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                      <  1:3378705/50‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                      <  2:1897/50‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                      <  1:1911363/50‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                      <  2:3549135/50‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       >  1:4470573/1‑49 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       <  2:1148559/49‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       <  1:2161218/49‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       >  2:2838546/1‑49 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                       <  1:1818021/49‑1 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                    >  2:1217759/1‑50 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                    >  1:694983/1‑50 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                    <  1:2730308/50‑1 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                    >  1:1779389/1‑50 (MQ=255)
        gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                   >  2:938899/1‑50 (MQ=255)
          aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                 >  2:2914209/1‑50 (MQ=255)
          aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                 >  2:883009/1‑50 (MQ=255)
          aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                 >  1:758928/1‑50 (MQ=255)
           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                <  1:2439246/50‑1 (MQ=255)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               <  2:1657939/50‑1 (MQ=255)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               <  1:1098212/50‑1 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  2:4331323/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  1:3137630/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  2:699746/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                               >  2:2438840/1‑49 (MQ=255)
              gccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                             <  2:655762/49‑1 (MQ=255)
              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                             <  2:1559191/50‑1 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                            <  1:1372988/50‑1 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                            >  2:1077376/1‑50 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                            <  1:4483218/49‑1 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           >  2:729303/1‑49 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                           >  1:3233245/1‑49 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                          <  2:2154743/50‑1 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                          >  2:4467803/1‑49 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                          >  1:3804093/1‑49 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                         <  2:2340651/50‑1 (MQ=255)
                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                          >  2:826797/1‑48 (MQ=255)
                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                          >  2:1221285/1‑48 (MQ=255)
                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa                        <  2:2861199/50‑1 (MQ=255)
                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                       <  1:4175383/50‑1 (MQ=255)
                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                       >  2:2636508/1‑50 (MQ=255)
                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                      >  1:2782407/1‑50 (MQ=255)
                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                      >  2:2027286/1‑50 (MQ=255)
                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                      >  1:3902012/1‑50 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                     <  1:1927831/48‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                     <  1:1641177/48‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   >  2:980622/1‑50 (MQ=255)
                         gACGCCGAAATTTCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   <  1:623690/49‑1 (MQ=38)
                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   <  1:2935999/49‑1 (MQ=255)
                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                   <  1:856817/49‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                  <  2:2151164/49‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 <  2:828896/50‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 <  2:2422593/50‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 <  2:2922299/50‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 <  2:2591761/50‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 <  2:4162706/50‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                 >  2:2884107/1‑50 (MQ=255)
                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                <  2:2605703/50‑1 (MQ=255)
                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                >  2:3859196/1‑49 (MQ=255)
                            gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa               <  2:3128498/50‑1 (MQ=255)
                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg             >  1:2304531/1‑50 (MQ=255)
                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg             <  2:3569220/50‑1 (MQ=255)
                              cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg             >  1:4110573/1‑50 (MQ=255)
                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            >  1:1957338/1‑50 (MQ=255)
                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            >  1:3861882/1‑50 (MQ=255)
                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg            >  2:1631892/1‑50 (MQ=255)
                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  2:305720/1‑50 (MQ=255)
                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:2804413/1‑50 (MQ=255)
                                aaaGTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCGGGAAGGc           <  2:4136237/50‑1 (MQ=38)
                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc          <  2:383334/50‑1 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  1:3089773/1‑48 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa         <  2:327927/50‑1 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa         <  2:1169759/50‑1 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa         <  1:644009/50‑1 (MQ=255)
                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg        <  2:9151/50‑1 (MQ=255)
                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg        <  2:4323378/50‑1 (MQ=255)
                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg        <  2:4469204/50‑1 (MQ=255)
                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGCGAAGGTGCTGGATGGc           >  2:3298243/1‑47 (MQ=37)
                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa       <  2:3911148/49‑1 (MQ=255)
                                     aCCAGCGCCGGCAAAAAAAGCGCGGAAGTGAAGGTGCTGGAAGGc           >  2:1474413/1‑45 (MQ=38)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc           >  2:3699149/1‑44 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt      <  2:931830/49‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt      <  1:1577144/49‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  1:3983321/50‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  2:3474011/50‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  1:1447461/50‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  2:268450/50‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     <  2:4148240/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg    <  1:1374369/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     >  2:4472732/1‑49 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc     >  2:2199066/1‑49 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga   >  1:2018965/1‑50 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga   >  2:149815/1‑50 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga   <  2:3362478/50‑1 (MQ=255)
                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac  <  1:2205144/50‑1 (MQ=255)
                                              |                                            
GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGAC  >  NC_000913/3091259‑3091349

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: ‑