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,3050AG95.7% 80.5 / ‑3.4 23S58G (AGC→GGC) rsmE16S rRNA m(3)U1498 methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (9/13);  total (10/13)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGA  >  NC_000913/3091260‑3091348
                                             |                                           
tttGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGc                                           >  1:1909705/1‑48 (MQ=255)
  tGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCaa                                         >  1:655120/1‑48 (MQ=255)
   gACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaa                                        >  1:1038775/1‑48 (MQ=255)
      ggTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCaaaaa                                      >  1:1160919/1‑47 (MQ=255)
        tAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGc                                  <  1:1984401/49‑1 (MQ=255)
          gCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGt                                <  1:694056/49‑1 (MQ=255)
             aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTgg                              >  1:1431230/1‑48 (MQ=255)
             aCCAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGa                             <  1:1117574/49‑1 (MQ=255)
              ccAGGTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGaa                            <  1:1814277/49‑1 (MQ=255)
                 ggTCTTTGACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGt                          >  1:1414743/1‑48 (MQ=255)
                        gACGCCGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGc                  >  1:194009/1‑49 (MQ=255)
                            ccGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGa              <  1:434914/49‑1 (MQ=255)
                             cGAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGaa             <  1:191214/49‑1 (MQ=255)
                              gAAATTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAg            <  1:1009605/49‑1 (MQ=255)
                               aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg           >  1:30594/1‑49 (MQ=255)
                               aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg           <  1:1079284/49‑1 (MQ=255)
                               aaaTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAgg           >  1:7581/1‑49 (MQ=255)
                                 aTTACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGcc         <  1:873502/49‑1 (MQ=255)
                                  ttACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCa        <  1:607344/49‑1 (MQ=255)
                                   tACCAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGc          >  1:1508637/1‑46 (MQ=255)
                                      cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc    <  1:2059363/49‑1 (MQ=255)
                                      cAGCGCCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATc    <  1:635768/49‑1 (MQ=255)
                                        gcgcCGGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATcga  <  1:953170/49‑1 (MQ=255)
                                             |                                           
TTTGACGGTAGCAACCAGGTCTTTGACGCCGAAATTACCAGCGCCAGCAAAAAAAGCGTGGAAGTGAAGGTGCTGGAAGGCCAGATCGA  >  NC_000913/3091260‑3091348

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