Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,597,860 C→T T36T (ACC→ACT rsmD → 16S rRNA (guanine(966)‑N(2))‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,597,8600CT82.3% 39.4 / 1.8 17T36T (ACC→ACTrsmD16S rRNA (guanine(966)‑N(2))‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (7/7);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

AATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACCGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGATTGCTTCG  >  NZ_CP009273/3597785‑3597933
                                                                           |                                                                         
aaTCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGa                                                             <  2:95126/90‑1 (MQ=255)
                           aGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGt                                  >  1:389659/1‑90 (MQ=255)
                                     cTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGAc                        <  1:213655/90‑1 (MQ=255)
                                      tCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACg                       >  1:129216/1‑90 (MQ=255)
                                       cccGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGc                      <  1:2910/90‑1 (MQ=255)
                                       cccGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGc                      >  2:59128/1‑90 (MQ=255)
                                          ggTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCa                   <  2:55618/90‑1 (MQ=255)
                                          ggTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCa                   <  2:57904/90‑1 (MQ=255)
                                          ggTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCa                   <  2:94335/90‑1 (MQ=255)
                                             tCCTGATAGCCCAGGTCTGCGCCCCACCACCGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATg                >  2:184965/1‑90 (MQ=255)
                                                    aGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGAt         <  2:139563/90‑1 (MQ=255)
                                                     gCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGAtt        >  1:246255/1‑90 (MQ=255)
                                                     gCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGAtt        >  1:309746/1‑90 (MQ=255)
                                                     gCCCAGGTCTGCGCCCCACCACTGACCGCGCACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGAtt        >  1:229247/1‑90 (MQ=255)
                                                           gTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGATTGCTTCg  >  2:119579/1‑90 (MQ=255)
                                                                    ccaccaCCGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTc                                      >  1:223884/1‑45 (MQ=255)
                                                                    ccaccaCCGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTc                                      <  2:223884/45‑1 (MQ=255)
                                                                           |                                                                         
AATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACCGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGATTGCTTCG  >  NZ_CP009273/3597785‑3597933

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