Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273883,4700TG52.9% 16.1 / 20.2 17intergenic (+55/+120)rcdA/ybjLDNA‑binding transcriptional regulator RcdA/aspartate:alanine antiporter
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base G (9/0);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.10e-04
Rejected as consensus: Frequency below/above cutoff threshold.

TTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGC  >  NZ_CP009273/883385‑883549
                                                                                     |                                                                               
tttGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGAt                                                                             <  1:526753/90‑1 (MQ=255)
                   gTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATgggg                                                           <  2:80214/89‑1 (MQ=255)
                   gTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATgggg                                                           >  1:80214/1‑89 (MQ=255)
                           gtagtaGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGGCATCAAAg                                                  >  1:224058/1‑90 (MQ=255)
                            tagtagATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCAAAAAAAGGGGGCATaaaaaa                                                 >  2:396251/1‑88 (MQ=255)
                               tagATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGaaaa                                              <  1:513967/90‑1 (MQ=255)
                                   tAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAAAGGGGGCATCAAAAAAAAACAg                                          >  2:124683/1‑90 (MQ=255)
                                   tAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAATGGGGAAATCAAAAAAAAGCAg                                          >  1:467254/1‑90 (MQ=255)
                                   tAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAGCAg                                          >  1:315611/1‑90 (MQ=255)
                                        tttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGGTCCACAAAAAAGGGACATCAAAAAAAAACAGTggga                                     >  2:634777/1‑88 (MQ=255)
                                        tttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGGCATCAAAAAAAAACAGTGGCa                                     >  2:214340/1‑90 (MQ=255)
                                         ttAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCAc                                    <  2:371143/90‑1 (MQ=255)
                                                aCAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCAAAAAAAAACAGGGGGACTAAttta                             >  1:688824/1‑88 (MQ=255)
                                                 cAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAg                            <  2:594496/90‑1 (MQ=255)
                                                                   tGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGaaa          <  2:315160/90‑1 (MQ=255)
                                                                          aaaTGGCGGGGGTGATGGGCTCCACAAAAAGGGGACATCAAAGAAAAACAGTGGCACTAATTAAAACTGATGCCCTGCGGGAAAGTTCTg   >  2:50706/1‑90 (MQ=255)
                                                                           aaTGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGc  <  1:50706/90‑1 (MQ=255)
                                                                                     |                                                                               
TTTGAGGATGGTTGAGAGGGTTGCAGGGTAGTAGATAAGTTTTAGATAACAAAAAACCCATCAACCTTGAACCGAAATGGCGGGGTTGATGGGCTCCACAAAATGGGGACATCAAAGAAAAGCAGTGGCACTAATTAAGACTGATGCCCTGCGGAAAAGTTCTGC  >  NZ_CP009273/883385‑883549

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