Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273262,0600CA58.3% 1.7 / 16.1 12V69F (GTT→TTT) BW25113_RS25425hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base A (7/0);  total (7/5)
Fisher's exact test for biased strand distribution p-value = 1.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGT‑CAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTG  >  NZ_CP009273/261994‑262149
                                                                  |                                                                                          
gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGc                                                                                        <  2:265437/71‑1 (MQ=255)
gCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAAACTGc                                                                                        >  1:265437/1‑71 (MQ=255)
    atCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGT‑c                                                                <  1:313999/90‑1 (MQ=255)
              tagtagtaACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTCCAAGGGGGGGCCGGGGAAGTAGT‑CAAACcgtggc                                                      >  2:212052/1‑88 (MQ=21)
                     aaCTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGGAAAAACCCCAAGGGGGGGCCGGGGAAAAAGT‑CAAACCGTGACCcttctt                                               >  2:308714/1‑90 (MQ=255)
                       cTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTTCAAGGGGGGGGCGGGTAAAAAGA‑CAAACCGTGACCCttttttc                                             >  1:427988/1‑88 (MQ=255)
                       cTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGGAAGGAAAAACCGCAAGGGGGGGCCCGGAAAGGAGTCCAACCTG‑GGCCCttcttct                                             >  2:382551/1‑89 (MQ=255)
                           cATCCTCGCAGCCATCGGGCGACTGGCGGGGACGTAAAAACTGCAAGGGGGGGCCGGGTAAATAGT‑CAACCTGTGGCCCTTTTTTGCTGa                                         >  1:236820/1‑90 (MQ=255)
                                     gCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGT‑CAACCTGTGACCATTCTTCGCTGACATCATCCtg                               <  1:212052/90‑1 (MQ=255)
                                       cATCGGGCGACTGGCGGGGACGTAAAAACCGCAAAGGGGGGCCGGGGAAAAAGG‑CAACCTGTGGCCCTTTTTCGCTGGCATCCTCCtggg                             >  2:120873/1‑89 (MQ=255)
                                                 cTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGT‑CAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACcac                   <  1:382068/90‑1 (MQ=255)
                                                           cGTAAAACCTGCAAGGTGGGACCGGATAAGTAGT‑CAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGaaa         <  1:333977/90‑1 (MQ=255)
                                                                  ccTGCAAGGTGGGACCGGATAAGTAGT‑CAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTg  <  1:507399/90‑1 (MQ=255)
                                                                  |                                                                                          
GCATATCCACGAGGTAGTAGTAACTGCCATCCTCGCAGCCATCGGGCGACTGGCGGGTACGTAAAACCTGCAAGGTGGGACCGGATAAGTAGT‑CAACCTGTGACCATTCTTCGCTGACATCATCCTGGTGGCTGACCACATCACTGAAACGCGGTG  >  NZ_CP009273/261994‑262149

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