Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273228,1910AG56.2% 6.3 / 17.2 16S195G (AGC→GGC) yafDendonuclease/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base A (1/6);  new base G (9/0);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 8.74e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.98e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA  >  NZ_CP009273/228104‑228261
                                                                                       |                                                                      
gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGc                                                                      <  2:527029/90‑1 (MQ=255)
     taGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAg                                                                 <  1:69289/90‑1 (MQ=255)
         aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTCTGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAgg                                                             >  2:27943/1‑90 (MQ=255)
         aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAAATCCTGGGGGCGGAGAAgg                                                             >  2:124341/1‑90 (MQ=255)
                 aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc                                                     <  1:247841/90‑1 (MQ=255)
                 aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc                                                     <  2:586968/90‑1 (MQ=255)
                      cTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATGGCTGGGGGCGGGGAAGGAGGAAagcggttt                                                >  2:98882/1‑86 (MQ=255)
                      cTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGGGGCGGAGAAGGGTGAACGCGGtat                                                >  1:244547/1‑90 (MQ=255)
                          tGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATTGCTGGGGCCGGGGAAGGGAGGAAACGTTatattg                                            >  1:424909/1‑88 (MQ=255)
                              gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTTATTCCTCGGGGCGCAGGAGGGGGAAAGCGGTTTATTGCttt                                        >  2:279416/1‑90 (MQ=255)
                                  aGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTgcgc                                    <  2:514037/90‑1 (MQ=255)
                                   gATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGAGTTTAATGCCCTGGGGCGTAGAAGGATGGACGCGGTATATTGCTTTtcggg                                   >  1:531892/1‑88 (MQ=255)
                                          caccacAGCGGCCCGGTCATTATGGCGGGAGGTTTCAATGCCTGGGGGCGTAGAAGGATGAAAGCGTTTTATCGCTTTTCGCGGGgaaag                            >  1:54062/1‑88 (MQ=255)
                                                 gCGGCCCGGTCATTATGGCGGGGGAGTTTAATGGCTGGGGGCGGAGGAGGAGGAAAGCGGTATATTGCTTTTTGCGGGGAAAGTTGCCgc                     >  2:199070/1‑90 (MQ=255)
                                                         gTCATTATGGCGGGAGGGTTTAATTCCTGGGGCCGGAGAAGGGTGGACGGGGTTTATCGCTTTTCGCGGGGAAAGTTGCTGCGCCCgggg             >  1:396595/1‑88 (MQ=255)
                                                                    gggAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga  <  1:401531/90‑1 (MQ=255)
                                                                                       |                                                                      
GTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA  >  NZ_CP009273/228104‑228261

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