Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273228,1810CT60.6% 1.3 / 28.6 33F191F (TTC→TTT) yafDendonuclease/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base C (2/11);  new base T (20/0);  total (22/11)
Fisher's exact test for biased strand distribution p-value = 4.03e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.26e-05
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA  >  NZ_CP009273/228094‑228261
                                                                                       |                                                                                
gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCaa                                                                                <  1:660024/90‑1 (MQ=255)
     aggATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCt                                                                           <  1:384809/89‑1 (MQ=255)
      ggATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGATTTAATGCCTg                                                                          >  2:865620/1‑90 (MQ=255)
        aTGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGa                                                                        <  1:842077/90‑1 (MQ=255)
          gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCCCCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGc                                                                      <  1:916508/90‑1 (MQ=255)
          gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATTCCTGGGGc                                                                      >  1:133148/1‑90 (MQ=255)
          gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATGCCTGGGGc                                                                      >  1:790820/1‑90 (MQ=255)
          gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGc                                                                      <  1:58722/90‑1 (MQ=255)
                 gTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAAATCCCGGGGCCGGAgct                                                               >  2:384809/1‑88 (MQ=255)
                   aaacAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGGCTCGGGCCGGCGAAgg                                                             >  2:795408/4‑90 (MQ=255)
                   aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTTATTACTGGGGCCGGAGaaag                                                             >  2:457490/1‑88 (MQ=255)
                   aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATTCCTGGGGGCGTGGAAgg                                                             >  2:720591/1‑90 (MQ=255)
                   aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATGCCTCGAGGCGCGGAAgg                                                             >  1:138554/1‑90 (MQ=255)
                   aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATACCTGGGGCCGGGGAAgg                                                             >  1:218119/1‑90 (MQ=255)
                        gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGGCGCAGGAGGATGaa                                                        >  1:905243/1‑90 (MQ=255)
                        gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAATTTAATTCCTTGGGGCGTGGAAGGAGGaa                                                        >  2:832123/1‑90 (MQ=255)
                        gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGaa                                                        <  2:627566/90‑1 (MQ=255)
                        gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGaa                                                        <  1:941664/90‑1 (MQ=255)
                           aCTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTTATGCCTGGGGGCGCAGGAGGAGGaagaa                                                     >  1:915756/1‑87 (MQ=255)
                                cTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGGATTTAATTCCTTGGGGCGGAGAAGGAGGAACGCGGtat                                                >  2:505307/1‑90 (MQ=255)
                                        gATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGGGGCGGAGAAGGATGAAAGCGGTTTATTGCttt                                        >  2:25291/1‑90 (MQ=255)
                                                    caccacAGCGGCCCGGTCATTATGGCGGGGGGGTTTAATGCCTGGAGCCGGAGAAAGAAGAAAGCGGTTTATCGCTTTGCGCGGGgaaag                            >  2:740684/1‑88 (MQ=255)
                                                    caccacAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGCCGGAGAAGGGAGAAAGCGGTATATCGCTTTGCGCGGGgaaag                            >  2:759934/1‑88 (MQ=255)
                                                       cacaGCGGCCCGGTCATTATGGCGGGGGGTTTTAAAGCCTGGAGCCGGAGAAGGAGGAAAGCGTTATATTGCTTTTCGCGGGGAAAGTCg                         >  1:270271/1‑90 (MQ=255)
                                                         caGCGGCCCGGTCATTATGGCGGGGGATTTTAATTCCTTGGGCCGGGGAAGGGGGAAAGCGGTTTATTGCTTTTTGCGGGGAAAGTCGCt                       >  1:357048/1‑90 (MQ=255)
                                                         caGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCt                       <  2:930406/90‑1 (MQ=255)
                                                         caGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCt                       <  2:133148/90‑1 (MQ=255)
                                                           gCGGCCCGGTCATTATGGCGGGAGATTTTAATGCCTGGAGCCGTAGAAGGATGAACGCGGTATATCGCTTTGCGCGGGGAAAGTtgctgc                     >  2:170565/1‑90 (MQ=255)
                                                                ccGGTCATTATGGCGGGAGATTTCAATGCCTGGGGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGGAATGTCGCTGCGcccg                >  2:141876/1‑88 (MQ=255)
                                                                       ttATGGCGGGAGATTTTAATTCCTGGGGGCGCAGAAGGATGAACGCGGTTTATCGCTTTTCGCGGGGAAAGTGGCTGCGCCAGGGGCggt         >  2:434535/1‑88 (MQ=255)
                                                                         aTGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGtttt       <  2:138554/90‑1 (MQ=255)
                                                                             cGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGGGCGTTTTACCg   >  2:55833/1‑90 (MQ=255)
                                                                              gggAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCga  <  2:270271/90‑1 (MQ=255)
                                                                                       |                                                                                
GGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGA  >  NZ_CP009273/228094‑228261

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