Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,464,5320TG57.1% 14.5 / 16.2 21*198G (TGA→GGA) bdcRTetR/AcrR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base G (12/0);  total (15/6)
Fisher's exact test for biased strand distribution p-value = 1.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-05
Rejected as consensus: Frequency below/above cutoff threshold.

CTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTA  >  NZ_CP009273/4464443‑4464605
                                                                                         |                                                                         
cTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGt                                                                           <  1:333796/90‑1 (MQ=255)
      aGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCg                                                                     <  1:70208/90‑1 (MQ=255)
                ggggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTTGGGGGGATGCGGGGCTTGAcc                                                           >  1:261981/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAacgc                                         >  2:122263/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                                         >  1:84329/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAAGCGGGGAAgcgc                                         >  2:231158/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTGGGAGGGATGCGGGCCTTTACCCGAAAGGCGGGGAAacgc                                         >  2:183498/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc                                         >  2:346504/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTTTTTAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc                                         >  2:165394/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc                                         >  2:114996/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTCTTTAGGGGGGATGCGGGCCTTGACCCGAAAGGGGGGGAAgagc                                         >  1:75130/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGAGATAAAAAAAATTTTTTAGGAATGAAGGGGGTGTTGACCCGAAAGGCGGGgaaaagc                                         >  1:330577/1‑86 (MQ=255)
                                    tctGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTCTTTAGGAGGGGTGCGGGGCCTGACCCGAAAGGCGGGGAAACGCCt                                       >  1:178445/1‑90 (MQ=255)
                                    tctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGaa                                              <  1:358339/83‑1 (MQ=255)
                                    tctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGaa                                              >  2:358339/1‑83 (MQ=255)
                                      tgtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGc                                     <  2:261981/90‑1 (MQ=255)
                                       gtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGcc                                    >  2:51432/1‑90 (MQ=255)
                                               gCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGGGGAGTGATGCGGCCTTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGtt                            <  2:330577/90‑1 (MQ=255)
                                               gCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGtt                            <  2:295346/90‑1 (MQ=255)
                                                        gCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAGGGCCTGCCGATAAGTTGTGATAAGa                   >  1:358445/1‑90 (MQ=255)
                                                            gggAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGGTAAGTTGTGATAAGACaat               >  2:95740/1‑90 (MQ=255)
                                                                         accTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTa  <  1:165394/88‑1 (MQ=255)
                                                                                         |                                                                         
CTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTA  >  NZ_CP009273/4464443‑4464605

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