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% 37.2 / ‑0.1 17T36T (ACC→ACTrsmD16S rRNA (guanine(966)‑N(2))‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (10/4);  total (12/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACCGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGATTGCTTCGCCGG  >  NZ_CP009273/3597780‑3597937
                                                                                |                                                                             
ggCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTa                                                                      <  1:348851/90‑1 (MQ=255)
        ccGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAAc                                                              <  2:901140/90‑1 (MQ=255)
                  ggCGGGCAGTGGCGAGGCCGGAAACTCCCGGTTCCTGATAACCCCGGGCTGCGCCCCACCCCTGACCCCGTACGCGAAACGTTGTTTAAc                                                    >  2:808734/1‑90 (MQ=255)
                     gggCAGGGGCGAGGCCCTAAACTCCCGGGTCCTTAAAGCCCACGGCTTCGCCCCACCACCGACCGCGGACGCGGAACGCTTTTTTAatgg                                                 >  2:331545/1‑90 (MQ=255)
                       gCAGTGGCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAActggct                                               <  1:78183/90‑1 (MQ=255)
                          gTGGCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACtggctggc                                            >  1:834068/1‑90 (MQ=255)
                            ggCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTc                                          >  2:607096/1‑90 (MQ=255)
                                 ggCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACCGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTc                                     >  1:794972/1‑90 (MQ=255)
                                      tAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCAttgt                                >  1:610115/1‑90 (MQ=255)
                                             ccGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGcc                         >  1:177679/1‑90 (MQ=255)
                                               ggTTCCTGATAGCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCa                       <  1:448103/90‑1 (MQ=255)
                                                  tCCTGATAGCCCAGGTCTGCGCCCCACCACCGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATg                    <  2:794972/90‑1 (MQ=255)
                                                          gCCCAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCCCCGGTCATTGTTGACGCCCAATGTCTGGAtt            >  2:41981/1‑90 (MQ=255)
                                                             cAGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGATTGCt         >  2:684980/1‑90 (MQ=255)
                                                              aGGTCTGCGCCCCACCACTGACCGCGTACGCGAAACGTTGTTCAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGATTGCtt        >  2:790117/1‑90 (MQ=255)
                                                                    gcgcCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGATTGCTTCGCCgg  >  2:115307/1‑90 (MQ=255)
                                                                    gcgcCCCACCACTGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGATTGCTTCGCCgg  >  1:249820/1‑90 (MQ=255)
                                                                                |                                                                             
GGCCAAATCCGCATTATTGGCGGGCAGTGGCGAGGCCGTAAACTCCCGGTTCCTGATAGCCCAGGTCTGCGCCCCACCACCGACCGCGTACGCGAAACGTTGTTTAACTGGCTGGCTCCGGTCATTGTTGACGCCCAATGTCTGGATTGCTTCGCCGG  >  NZ_CP009273/3597780‑3597937

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