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,3050AG87.6% 362.8 / 38.1 121S58G (AGC→GGC) rsmE16S rRNA m(3)U1498 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (7/8);  new base G (58/48);  total (65/56)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 122 total aligned reads displayed.

ATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091256‑3091351
                                                 |                                              
aTTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCa                                                <  1:84106/50‑1 (MQ=255)
 ttGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAg                                               <  2:3703190/50‑1 (MQ=255)
   gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              <  2:5374836/49‑1 (MQ=255)
   gTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                              >  1:3888300/1‑49 (MQ=255)
   gTTTGACGGAAGCAACCAGGTCTTTGACGCCGACATTACCAGCGCCGGc                                              >  2:3703326/1‑49 (MQ=38)
      tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           >  1:557445/1‑49 (MQ=255)
      tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                           <  1:2350718/49‑1 (MQ=255)
      tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCaaa                                           <  2:4398021/49‑1 (MQ=255)
        aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         <  1:3753338/49‑1 (MQ=255)
        aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:5382530/1‑49 (MQ=255)
        aCGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         <  2:1041653/49‑1 (MQ=255)
         cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                       <  1:4368400/50‑1 (MQ=255)
         cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaaaa                                       <  1:3152660/50‑1 (MQ=255)
         cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  2:4083509/1‑48 (MQ=255)
         cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  2:2399948/1‑48 (MQ=255)
         cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:3460044/1‑48 (MQ=255)
         cGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                         >  1:4582075/1‑48 (MQ=255)
         cGGTAGCAACCAGGTCTTTGACGCAGAAATTACCAGCGCCGGCaaaaa                                         >  2:1461778/1‑48 (MQ=38)
          ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      <  2:4578382/50‑1 (MQ=255)
          ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      <  2:3035670/50‑1 (MQ=255)
          ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAg                                      >  2:4929008/1‑50 (MQ=255)
           gTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                     <  2:3811783/50‑1 (MQ=255)
           gTAGCAACCAGGTATTTGCCGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                     <  1:3919910/50‑1 (MQ=255)
            tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCg                                    >  1:5011585/1‑50 (MQ=255)
             aGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   >  2:4945200/1‑50 (MQ=255)
              gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                   <  1:2578313/49‑1 (MQ=255)
              gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1398859/1‑50 (MQ=255)
              gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1595759/1‑50 (MQ=255)
               cAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTg                                  >  1:1548486/1‑49 (MQ=255)
                aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  2:2688049/1‑49 (MQ=255)
                aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  2:4788840/50‑1 (MQ=255)
                aaCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                                <  2:559854/50‑1 (MQ=255)
                 cccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                                 >  2:3614357/2‑48 (MQ=255)
                 aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  2:2783688/50‑1 (MQ=255)
                 aCCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTgg                                 >  1:2328986/1‑48 (MQ=255)
                 aCCAGGTCTTTGACGCCGAAAGTACCAGCGCCGGCAAAAAAAGCGTGGaa                               <  1:2060755/50‑1 (MQ=38)
                  ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  1:1674216/50‑1 (MQ=255)
                  ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              <  2:2036910/50‑1 (MQ=255)
                  ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAg                              >  1:2988985/1‑50 (MQ=255)
                  ccAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAg                              <  1:3615174/50‑1 (MQ=255)
                   cAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                             <  2:1675711/50‑1 (MQ=255)
                   cAGGCCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGCGGAAGt                             >  1:2446451/1‑50 (MQ=39)
                    aGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAGCGTGGAAGTg                            >  2:191516/1‑49 (MQ=38)
                     ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTg                            >  2:5249183/1‑49 (MQ=255)
                     ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  1:2202967/50‑1 (MQ=255)
                     ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGa                           <  2:4917500/50‑1 (MQ=255)
                       tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAGAGCGTGGAAGTGAAg                         >  2:2284692/1‑50 (MQ=39)
                       tCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAg                         >  1:4367631/1‑50 (MQ=255)
                        cTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        >  2:3642779/1‑50 (MQ=255)
                         tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:2346862/1‑50 (MQ=255)
                         tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:5362593/1‑50 (MQ=255)
                         tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  2:3877358/1‑50 (MQ=255)
                         tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  2:1215642/1‑50 (MQ=255)
                         tttGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       >  1:1163114/1‑50 (MQ=255)
                          ttGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTg                      >  2:1893706/1‑50 (MQ=255)
                           tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAgg                        <  1:1568639/47‑1 (MQ=255)
                           tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  1:594942/48‑1 (MQ=255)
                           tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGt                       <  2:5280449/48‑1 (MQ=255)
                           tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  2:3989031/49‑1 (MQ=255)
                           tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  1:2502245/49‑1 (MQ=255)
                           tGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     >  2:1274190/1‑50 (MQ=255)
                           tGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAgg                        <  1:1252942/47‑1 (MQ=255)
                           tGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTg                      <  2:4528036/49‑1 (MQ=255)
                            gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     <  1:1408848/49‑1 (MQ=255)
                            gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                     <  1:2921515/49‑1 (MQ=255)
                             aCGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTg                   >  1:5147252/1‑50 (MQ=255)
                              cGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  <  1:3962045/50‑1 (MQ=255)
                               gCCGAATTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  >  2:3178475/1‑49 (MQ=38)
                               gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  >  2:1197880/1‑49 (MQ=255)
                               gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTgg                  >  1:1988179/1‑49 (MQ=255)
                               gCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa                 <  1:5182451/50‑1 (MQ=255)
                               gCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCGGGa                 <  2:442028/50‑1 (MQ=255)
                                ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGaa                <  2:2273862/50‑1 (MQ=255)
                                 cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               >  2:3611718/1‑50 (MQ=255)
                                 cGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg               >  1:3196964/1‑50 (MQ=255)
                                  gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              <  2:4471970/50‑1 (MQ=255)
                                  gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              <  1:1907788/50‑1 (MQ=255)
                                  gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg              >  1:5225806/1‑50 (MQ=255)
                                   aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:17598/1‑50 (MQ=255)
                                   aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             <  1:2532334/50‑1 (MQ=255)
                                   aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:4750294/1‑50 (MQ=255)
                                   aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             <  1:1831861/50‑1 (MQ=255)
                                    aaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:3031480/1‑49 (MQ=255)
                                    aaTTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:4085270/1‑49 (MQ=255)
                                     aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:3290938/1‑48 (MQ=255)
                                     aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:1785931/1‑48 (MQ=255)
                                     aTTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa           <  1:1106948/50‑1 (MQ=255)
                                      ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:2619659/1‑47 (MQ=255)
                                      ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAg          <  1:2094877/50‑1 (MQ=255)
                                       tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:1896194/1‑46 (MQ=255)
                                        aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  2:5106026/49‑1 (MQ=255)
                                        aCCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGa         <  1:4209492/49‑1 (MQ=255)
                                         ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  1:5094205/1‑44 (MQ=255)
                                         ccAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc       <  1:3460902/50‑1 (MQ=255)
                                         ccAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc             >  2:2070820/1‑44 (MQ=255)
                                          cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  2:4702724/50‑1 (MQ=255)
                                          cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcg      <  2:4939352/50‑1 (MQ=255)
                                           aGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga     >  2:1086679/1‑50 (MQ=255)
                                            gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgac    <  2:2476790/50‑1 (MQ=255)
                                              gcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcgacga  >  1:2626634/1‑50 (MQ=255)
                                                 |                                              
ATTGTTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGACGA  >  NC_000913/3091256‑3091351

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