Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,164,0172.C50.0% 7.5 / 10.8 8coding (272/633 nt)comRTetR family copper‑responsive transcriptional repressor ComR
Reads supporting (aligned to +/- strand):  ref base . (2/2);  new base C (2/2);  total (4/4)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGATCATGAAGCAGCCTGCCGGGGTGTCTTTGCTGGTAAAGCAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTT‑‑TCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAA  >  NZ_CP009273/1163928‑1164089
                                                                                           |                                                                        
tgATCATGAAGCAGCCTGCCGGGGTGTCTTTGCTTGTAAAGCAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCttt                                                                            >  1:168530/1‑90 (MQ=255)
    cATGAAGCAGCCTGCCGGGGTGTCTTTGCTGGTAAAGCAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTT‑‑tctc                                                                      <  1:161773/90‑1 (MQ=255)
                        tgtCTTTGCTGGTAAAGCAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTT‑‑TCTCTTCACAAAACAGCTGTGCtt                                                  >  1:130705/1‑90 (MQ=255)
                                   gTAAAGCAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTTTCTCTCTTCACAAAACAGCTGTGCTTCATGCTTAg                                         >  2:304602/1‑90 (MQ=255)
                                   gTAAAGCAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTTTCTCTCTTCACAAAACAGCTGTGCTTCATGCTTAg                                         >  2:374363/1‑90 (MQ=255)
                                                                  cagccagcGCAGACTCCACGCTTT‑‑TCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCaa                                   <  1:107945/63‑1 (MQ=255)
                                                                  cagccagcGCAGACTCCACGCTTT‑‑TCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCaa                                   >  2:107945/1‑63 (MQ=255)
                                                                     ccagcGCAGACTCCACGCTTT‑‑TCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAATGATCGATATAGCGGTCGAGAACGGCAc     <  2:168530/90‑1 (MQ=255)
                                                                          gcAGACTCCACGCTTTTCTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGaa  <  1:304602/90‑1 (MQ=255)
                                                                          gcAGACTCCACGCTTTTCTCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGaa  <  1:374363/90‑1 (MQ=255)
                                                                                           |                                                                        
TGATCATGAAGCAGCCTGCCGGGGTGTCTTTGCTGGTAAAGCAGTTGGCGATGGCAGCAAAATAGTCAGCCAGCGCAGACTCCACGCTTT‑‑TCTCTTCACAAAACAGCTGTGCTTCATGCTTAGCGGCAAAACGATCGATATAGCGGTCGAGAACGGCACGAA  >  NZ_CP009273/1163928‑1164089

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