Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273228,1750AG53.3% 9.4 / 18.8 15G189G (GGA→GGGyafDendonuclease/exonuclease/phosphatase family protein
Reads supporting (aligned to +/- strand):  ref base A (0/7);  new base G (8/0);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 1.55e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.72e-04
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTT  >  NZ_CP009273/228086‑228256
                                                                                         |                                                                                 
ttCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGga                                                                                   <  2:579310/90‑1 (MQ=255)
        gggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCaa                                                                           <  1:645326/90‑1 (MQ=255)
         ggCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGCGTTCaag                                                                          >  1:3207/1‑89 (MQ=255)
                  gTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGc                                                                 <  2:29590/90‑1 (MQ=255)
                  agcTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGc                                                                 <  1:530766/88‑1 (MQ=255)
                   tCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTTATGCCTGGGGcc                                                                >  2:661403/1‑90 (MQ=255)
                         gTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATTCCCGGGGGCGGAGaa                                                          >  2:88424/1‑90 (MQ=255)
                           aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATGTCTGGGGCCGGGGAAgg                                                        >  1:585/1‑90 (MQ=255)
                           aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGGTTTTATGCCTGGGGCCGGAGaaag                                                        >  1:672669/1‑88 (MQ=255)
                           aaGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGAGTTTAATTCCTGGGGGCGGGGAAgg                                                        >  2:360193/1‑90 (MQ=255)
                              cAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGATTTTAATGCCTGGGGACGTAGAAgggaa                                                     >  2:530766/1‑87 (MQ=255)
                                gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCTGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGaa                                                   <  2:515677/90‑1 (MQ=255)
                                gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGga                                                                                   >  2:205023/1‑58 (MQ=255)
                                gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGga                                                                                   <  1:205023/58‑1 (MQ=255)
                                gTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGGGGTTTTAATTCCTGGAGCCGGAGAAAGATGaa                                                   >  2:337973/1‑90 (MQ=255)
                                     ttCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGt                                              <  2:141832/90‑1 (MQ=255)
                                                                           gTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTg        <  2:349549/90‑1 (MQ=255)
                                                                                 aTGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGtttt  <  2:557151/90‑1 (MQ=255)
                                                                                         |                                                                                 
TTCAGTCTGGGCGTGGATGTCTATAGTAAGCAGTTACTTCCTATTGGCGATCAGATAGCTCACCACAGCGGCCCGGTCATTATGGCGGGAGATTTCAATGCCTGGAGCCGTAGAAGGATGAACGCGTTATATCGCTTTGCGCGGGAAATGTCGCTGCGCCAGGTGCGTTTT  >  NZ_CP009273/228086‑228256

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