Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273883,4700TG56.2% 8.8 / 14.9 16intergenic (+55/+120)rcdA/ybjLDNA‑binding transcriptional regulator RcdA/aspartate:alanine antiporter
Reads supporting (aligned to +/- strand):  ref base T (2/5);  new base G (9/0);  total (11/5)
Fisher's exact test for biased strand distribution p-value = 4.81e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGT  >  NZ_CP009273/883392‑883552
                                                                              |                                                                                  
aTGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTcc                                                                         >  2:74257/1‑90 (MQ=255)
         agGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGGTGGGGTCCACAAAAAgg                                                                >  1:417288/1‑90 (MQ=255)
           ggTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATgggg                                                              <  2:93990/90‑1 (MQ=255)
             ttGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGGGGTGGGCCCCAAAAAAAAGGGAc                                                            >  1:315117/1‑90 (MQ=255)
                              aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGTGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTg                                           >  2:264691/1‑90 (MQ=255)
                              aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTg                                           <  1:359204/90‑1 (MQ=255)
                              aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGGCATCAAAGAAAAGCAGTg                                           >  2:108148/1‑90 (MQ=255)
                              aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGGCATCAAAGAAAAGCAggg                                           >  2:228088/1‑88 (MQ=255)
                              aGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGGCATCAAAAAAAAGCAggg                                           >  2:490288/1‑88 (MQ=255)
                                 tttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAATGGGGGCATCAAAGAAAAACAGTggga                                        >  2:359204/1‑88 (MQ=255)
                                  ttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCAc                                       <  1:430205/90‑1 (MQ=255)
                                        aaCAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTa                                 <  2:458882/90‑1 (MQ=255)
                                                             gAACCGAAATGGCGGGGGTGATGGGGTCCCCAAAAAGGGGACATCAAAAAAAAACAGGGGCACTAAATAAAACTGATGCCCCGCGGaaaa            >  1:24411/1‑90 (MQ=255)
                                                             gAACCGAAATGGCGGGGGTGATGGGCTCCCCAAAAAAGGGACATCAAAAAAAAACAAGGGGACTAATTAAGACTGATGGCCTGCGGgaaa            >  1:353115/1‑90 (MQ=255)
                                                                  gAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGGCATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGGAAAGtttt       >  1:130645/1‑88 (MQ=255)
                                                                       ggCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGt  <  1:504380/90‑1 (MQ=255)
                                                                              |                                                                                  
ATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGCGGT  >  NZ_CP009273/883392‑883552

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