Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,596,7190AG56.2% 0.4 / 13.0 16S267P (TCC→CCC) rsmC16S rRNA (guanine(1207)‑N(2))‑methyltransferase RsmC
Reads supporting (aligned to +/- strand):  ref base A (3/4);  new base G (9/0);  total (12/4)
Fisher's exact test for biased strand distribution p-value = 1.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.64e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACA  >  NZ_CP009273/4596632‑4596786
                                                                                       |                                                                   
cgcTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGgaga                                                                   <  1:205056/90‑1 (MQ=255)
           tGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGGGGATCatatta                                                        >  2:3715/1‑88 (MQ=255)
            gACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCATATCaa                                                       >  1:303546/1‑90 (MQ=255)
            gACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGAAGAACATATCaa                                                       >  2:172047/1‑90 (MQ=255)
               gcaccgcaccGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAAc                                                    <  1:237657/90‑1 (MQ=255)
                       ccgcGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGGTTATCATATCAAAACAACCTTTc                                            >  2:187674/1‑90 (MQ=255)
                                 gCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGATCATATCAAAACGACCCTTCACCTCGGgaa                                  >  2:225079/1‑90 (MQ=255)
                                          ccgcATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGc                         <  1:280183/90‑1 (MQ=255)
                                                      gggTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACtt             <  2:359666/89‑1 (MQ=255)
                                                       ggTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCATATCAAAACGACCTTTCACCTCGGAAAAAACGTTGCTGGGAAAGACTTc            >  1:340879/1‑90 (MQ=255)
                                                       ggTTTGCATCCCATCGTGGAACGGCGGGTTGGGGATGATCAAATTAAAAAGAACTTTTCCCCCGGGAAAAACGCTGGTCGCCAAAACTTc            >  2:282468/1‑90 (MQ=255)
                                                       ggTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc            >  1:29098/1‑90 (MQ=255)
                                                       ggTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTc            <  2:29098/90‑1 (MQ=255)
                                                       ggTTTGCATCCCATCGTGGAACGGCGGGGTGGGGATGATCATATCAAAAAGACCTTTTACCCCGGGAAAAACGTTTGTGGCAAAAAATTc            >  2:105546/1‑90 (MQ=255)
                                                            gCATCCCATCGTGGAACGGCGGGGTGGAGGTGGTTATATTAAAACGACCTTTCACCTCGGGAAAGACGTTGCTGGGGAAAAATTCCCCtt       >  1:313401/1‑90 (MQ=255)
                                                                 ccATCGTGGAACGGCGGGTTGGAGATGGTCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTTAaca  >  1:31266/1‑90 (MQ=255)
                                                                 ccATCGTGGAACGGCGGGGTGGGGATGATCATATCAAAACGACCTTTCCCCTCGGAAAAGACGTTGCTGGCAAAGACATCACCTTCAacc  >  1:359666/1‑89 (MQ=255)
                                                                                       |                                                                   
CGCTATTAAGATGACGCACCGCACCGCGAATCAGCGTTTGCGCCGCATCCAGGCTGGTTTGCATCCCATCGTGGAACGGCGGGTTGGAGATGATCATATCAAAACGACCTTTCACCTCGGAAAAGACGTTGCTGGCAAAGACTTCACCTTCAACA  >  NZ_CP009273/4596632‑4596786

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