Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273228,1910AG54.5% 13.1 / 26.5 22S195G (AGC→GGC) yafDendonuclease/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base A (1/9);  new base G (12/0);  total (13/9)
Fisher's exact test for biased strand distribution p-value = 2.01e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.45e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGA‑TTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA  >  NZ_CP009273/228104‑228261
                                                                                        |                                                                      
gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGA‑TTTCAATGCCTGGAGc                                                                      <  2:712081/90‑1 (MQ=255)
      aGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGG‑GTTTAATAGCTTGAGGCGTAGa                                                                >  1:196624/1‑90 (MQ=255)
      aGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGA‑TTTCAATGCCTGGAGCCGTAGa                                                                <  2:538496/90‑1 (MQ=255)
         aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGG‑TTTTAATGCCTGGGGGCGTAGAAgg                                                             >  2:269139/1‑90 (MQ=255)
         aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGA‑TTTTAATGCCTGGGGGCGGAGAAgg                                                             >  2:180070/1‑90 (MQ=255)
         aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGA‑GTTTAATATCTTGGGCCGGCGAAgg                                                             >  1:129518/1‑90 (MQ=255)
              gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGA‑TTTTAATGCCTGGGGCCGTAGAAGGATGaa                                                        >  2:219773/1‑90 (MQ=255)
              gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGA‑TTTTAAAGCCTGGGGCCCTAGAAGGAGGaa                                                        >  2:734052/1‑90 (MQ=255)
                 aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGA‑TTTCAATGCCTGGAGCCGTAGAAGGATGAAcgc                                                     <  1:607860/90‑1 (MQ=255)
                        aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGA‑TTTCAATGCCTGGa                                                                        <  1:404224/64‑1 (MQ=255)
                        aTTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGA‑TTTCAATGCCTGGa                                                                        >  2:404224/1‑64 (MQ=255)
                                         tCACCACAGCGGCCCGGTCATTATGGCGGGGGA‑GTTTAATTCCCGGGGCCGGGGAAGGGGGAAGCCGTTATTTTGCCTTTCGCGGGgaaa                             >  1:390524/1‑89 (MQ=255)
                                                 gCGGCCCGGTCATTATGGCGGGGTA‑TTTTAAAGGCCGGGGCCGGGGAAAGGTGAAAGCGGTATTTTGCGTTTCGGGGGGGATGGtgctgc                     >  1:350753/1‑90 (MQ=255)
                                                 gCGGCCCGGTCATTATGGCGGGGGG‑GTTTAATGCCTGGGGGCGGAGAAGGGGGGAAGCGGTTTATTGTTTTTCGCGGGGGAAGGCGCTgc                     >  1:634427/1‑90 (MQ=255)
                                                 gCGGCCCGGTCATTATGGCGGGAGG‑TTTTAATGCCTGGGGGCGGAGAAGGAAGAAAGCGGTATATCGCTTTTCGCGGGGAAAGGCGGTgc                     >  2:657699/1‑90 (MQ=255)
                                                 gCGGCCCGGTCATTATGGCGGGAGAGTTTCAA‑GGCTCGGGCCGCAGGAGGGTGGAGGAGTTATATTGCTTTGCGCGGGGAAGGgcgcggg                     >  1:684375/1‑87 (MQ=255)
                                                 gCGGCCCGGTCATTATGGCGGGAGA‑TTTTAATGCCTGGGGCCGTAGAAGGGTGAAAGCGTTATATCGCTTTGCGCGGGGAATGTCGCTgc                     >  2:213991/1‑90 (MQ=255)
                                                   ggCCCGGTCATTATGGCGGGAGA‑TTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTgcgc                   <  1:639702/90‑1 (MQ=255)
                                                    gCCCGGTCATTATGGCGGGAGA‑TTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGcc                  <  1:250982/90‑1 (MQ=255)
                                                      ccGGTCATTATGGCGGGAGA‑TTTTAATGCCTGGGGCCGGAGAAGGGTGAAAGCGGTATATCGCTTTTCGCGGGGAAAGTCGCTGCGcccg                >  1:98717/1‑88 (MQ=255)
                                                               aTGGCGGGAGA‑TTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGtttt       <  1:573916/90‑1 (MQ=255)
                                                               aTGGCGGGAGA‑TTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGtttt       <  1:358067/90‑1 (MQ=255)
                                                                    gggAGA‑TTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga  <  1:96049/90‑1 (MQ=255)
                                                                    gggAGA‑TTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga  <  2:443400/90‑1 (MQ=255)
                                                                                        |                                                                      
GTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGA‑TTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA  >  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 < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑