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

GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑3091351
                                              |                                              
gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              <  1:2200515/49‑1 (MQ=255)
gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              >  1:6422611/1‑49 (MQ=255)
gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              <  1:635811/49‑1 (MQ=255)
gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              >  1:4260631/1‑49 (MQ=255)
 tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCa                                             <  1:1168607/49‑1 (MQ=255)
   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           >  1:6762413/1‑49 (MQ=255)
   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           <  1:370054/49‑1 (MQ=255)
   tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           <  1:3829351/49‑1 (MQ=255)
    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  1:3868231/1‑49 (MQ=255)
    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  1:6419043/1‑49 (MQ=255)
    gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaa                                          >  1:4458079/1‑49 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                        <  1:6310947/50‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaa                                        <  1:4681893/50‑1 (MQ=255)
     aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:4672680/1‑49 (MQ=255)
      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                       <  1:62244/50‑1 (MQ=255)
      cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:3862409/1‑48 (MQ=255)
       ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  1:6858285/1‑50 (MQ=255)
        gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                     >  1:3988919/1‑50 (MQ=255)
         tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCg                                    >  1:3405527/1‑50 (MQ=255)
           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:4953526/1‑50 (MQ=255)
           gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1628641/1‑50 (MQ=255)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 <  1:1966510/50‑1 (MQ=255)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 <  1:2727425/50‑1 (MQ=255)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 <  1:31463/50‑1 (MQ=255)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1830871/1‑49 (MQ=255)
            cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1526981/1‑49 (MQ=255)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCTGGTAAAAAAAGCGTgg                                 >  1:7615760/1‑49 (MQ=38)
             aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  1:1707840/50‑1 (MQ=255)
              aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  1:3330142/1‑48 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:7398494/50‑1 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:4159278/50‑1 (MQ=255)
               ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:7535365/50‑1 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:5679580/1‑50 (MQ=255)
                cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:990727/1‑50 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGTGCCGGCAAAAAAAGCGTGGAAGTg                            <  1:2833583/50‑1 (MQ=38)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             >  1:753722/1‑49 (MQ=255)
                 aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            <  1:1640678/50‑1 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:7151854/1‑49 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  1:6857812/1‑49 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  1:6162366/50‑1 (MQ=255)
                  ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  1:5539121/50‑1 (MQ=255)
                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa                          <  1:3787559/50‑1 (MQ=255)
                   gTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGaa                          <  1:4526043/50‑1 (MQ=255)
                    tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         <  1:1405546/50‑1 (MQ=255)
                     cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        >  1:2891235/1‑50 (MQ=255)
                      tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:6897122/1‑50 (MQ=255)
                      tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:6247271/1‑50 (MQ=255)
                       ttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      >  1:5417012/1‑50 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  1:130244/48‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  1:1565523/49‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  1:463608/49‑1 (MQ=255)
                        tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  1:4099299/49‑1 (MQ=255)
                         gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    >  1:1359596/1‑50 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCt                    <  1:4180824/49‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   <  1:5670705/50‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   <  1:6707372/50‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   <  1:6864865/50‑1 (MQ=255)
                          aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   >  1:42272/1‑50 (MQ=255)
                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  <  1:5791953/50‑1 (MQ=255)
                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   >  1:6533668/1‑49 (MQ=255)
                           cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   >  1:2536591/1‑49 (MQ=255)
                             ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  >  1:5860362/1‑48 (MQ=255)
                             ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGaa                <  1:5987296/50‑1 (MQ=255)
                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              >  1:7191686/1‑50 (MQ=255)
                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              >  1:1333654/1‑50 (MQ=255)
                               gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              >  1:242175/1‑50 (MQ=255)
                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:4294936/1‑50 (MQ=255)
                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:2887647/1‑50 (MQ=255)
                                aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             <  1:1723258/50‑1 (MQ=255)
                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:7504607/50‑1 (MQ=255)
                                 aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc            <  1:2503969/50‑1 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:423460/1‑48 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:4859075/50‑1 (MQ=255)
                                  aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:4301280/50‑1 (MQ=255)
                                   ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  1:6564583/50‑1 (MQ=255)
                                    tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:1680748/1‑46 (MQ=255)
                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  1:6310088/49‑1 (MQ=255)
                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  1:1438479/49‑1 (MQ=255)
                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  1:1397707/49‑1 (MQ=255)
                                     aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  1:108440/49‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:992375/1‑44 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGAt        <  1:4899623/49‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:3971184/50‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:3531413/50‑1 (MQ=255)
                                      ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:1523078/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  1:4269762/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  1:820238/50‑1 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       >  1:6446934/1‑49 (MQ=255)
                                       cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       >  1:7055215/1‑49 (MQ=255)
                                        aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     <  1:1671737/50‑1 (MQ=255)
                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    <  1:1846329/50‑1 (MQ=255)
                                         gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    <  1:1111306/50‑1 (MQ=255)
                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  >  1:5384934/1‑50 (MQ=255)
                                           gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  <  1:837673/50‑1 (MQ=255)
                                              |                                              
GTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091259‑3091351

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