Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP00927315,4100TG53.3% 9.8 / 14.9 15intergenic (+112/‑35)dnaJ/BW25113_RS00080molecular chaperone DnaJ/IS4‑like element IS421 family transposase
Reads supporting (aligned to +/- strand):  ref base T (1/6);  new base G (8/0);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 1.40e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.15e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC  >  NZ_CP009273/15321‑15473
                                                                                         |                                                               
ggCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt                                                                 <  1:294636/90‑1 (MQ=255)
   aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtgt                                                              <  2:250311/90‑1 (MQ=255)
   aGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtgt                                                              <  2:56738/90‑1 (MQ=255)
    ggCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtgtg                                                             <  1:229426/90‑1 (MQ=255)
           gTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                                      <  2:107891/90‑1 (MQ=255)
                   aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCTAACTTAAGGGGTGTGGGATTACGCCTGgta                                              >  2:140892/1‑90 (MQ=37)
                   aGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                                              <  1:189145/90‑1 (MQ=255)
                    gggTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGGGACTCCCCTAAACGCTAAATTTAAGGGTGTGGGGTTTTGCCTGatat                                             >  1:107891/1‑90 (MQ=255)
                      gTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGGTTACTCCCCTAAACGCTAACTTTAAGGGTGTGGGGTTTCGCCTGGTATGa                                           >  2:55419/1‑90 (MQ=255)
                                  tatGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCCAAACGCCAAATTAAAGGGTGGGGGATTTAGCCTGATATGATTTTAAGTGGCg                               >  1:71937/1‑90 (MQ=255)
                                    tGCGAGCACCTGTAAAGTGGCGGGGGTCACTCCCCCAAACGCGAACTTTAGGGTTGTGGGATTACGCCCGGTATGATTTTAAGTGGCGAt                             >  2:336160/1‑90 (MQ=255)
                                       gAGCACCTGTAAAGTGGCGGGGGTCACTCCCCTAAACGCTAACTTAAAGGGTGTGGGATTACGCCCGATATGATTTTAAGTGCCGGTGaa                          >  2:265867/1‑90 (MQ=255)
                                       gAGCACCTGTAAAGTGGCGGGGGTCACTCCCCCAAGCGCTAACTTAAGGGGTGGGGGGTTTTGCCTGATATGATTTTACGTGCCGCTGaa                          >  1:44520/1‑90 (MQ=18)
                                       gAGCACCTGTAAAGTGGCGGGGGGCACTCCCCTAAGCGCTAAATTAAAGGGTGTGGGGTTAAGCCTGGTATTATTTTAAGTGCCGATGaa                          >  1:262746/1‑90 (MQ=255)
                                                      ggCGGGGGTCACTCCCCCAAACGCTAACTTAAGGGGTGGGGGGTTTCGCCCGATATGATTTTACATGCCGGGGAATTAATCTCCCCATaa           >  1:143839/1‑90 (MQ=9)
                                                               cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGc  <  1:336160/90‑1 (MQ=33)
                                                                                         |                                                               
GGCAGGCCTGGGTAAAAATAGGGTGCGTTGAAGATATGCGAGCACCTGTAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGC  >  NZ_CP009273/15321‑15473

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