Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,464,5320TG53.1% 25.4 / 29.8 33*198G (TGA→GGA) bdcRTetR/AcrR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (5/10);  new base G (17/0);  total (23/10)
Fisher's exact test for biased strand distribution p-value = 4.65e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATT  >  NZ_CP009273/4464450‑4464619
                                                                                  |                                                                                       
ggCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCgg                                                                                  <  1:170420/90‑1 (MQ=255)
    cGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCt                                                                              <  1:288797/90‑1 (MQ=255)
          gggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGAGGGATGCGGGGCCTGAccc                                                                        >  2:52503/1‑90 (MQ=255)
           ggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTTTTGAGGAGTGATGCGGGCCTTGACCCg                                                                       >  1:360437/1‑90 (MQ=255)
                 ctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGGGGGATGCGGGGCTTGACCCGAAAGGc                                                                 >  2:150857/1‑90 (MQ=255)
                      gAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAATATTTTTTTGGGGGGATGCGGGGCTTGACCCGAAAGGCgggga                                                            >  2:443958/1‑88 (MQ=255)
                        acaacTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAg                                                          >  1:4839/1‑90 (MQ=255)
                        acaacTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAAGCGGGgaaa                                                          >  2:61547/1‑89 (MQ=255)
                         caacTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTTTTGAGGAGGGATGCGGGCCTTTACCCCAAAGGCGGGaaaac                                                         >  1:389917/1‑88 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAATATTCTTGAGGAGAGAGGCGGGCCTTGACCCGAAAGGGGGGGAAgggc                                                       >  1:374653/1‑90 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTTACCCGAAAGGCGGGgaaaccc                                                       >  1:359921/1‑86 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                                                       >  2:57242/1‑90 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAAGGGGGGAAacgc                                                       >  2:211598/1‑90 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGTGGGGGGATGCGGGGCCTGGCCCGAAAGGCGGGGAAacgc                                                       >  2:321150/1‑90 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGGGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc                                                       >  2:91135/1‑90 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc                                                       >  2:204567/1‑90 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGt                                                                                         >  1:5290/1‑56 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGt                                                                                         <  2:5290/56‑1 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc                                                       >  1:207531/1‑90 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGAGGCGGGCCTTGAGCCGAAAGGCGGGAAAgcgc                                                       >  1:35060/1‑90 (MQ=255)
                           acTCTGTGCAACAGCTGCAATGGCGGGGGAAGACATAAAAAACATTCTTGAGGAAGGGTGTGGGTCTTAAGCCGAAAGGGGGGaaaatgg                                                       >  2:319197/1‑86 (MQ=255)
                             tctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTCTTGAGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAACGccc                                                     >  1:332166/1‑89 (MQ=255)
                                gtgCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGGGGGGAAGcgcccgcc                                                  >  1:448341/1‑90 (MQ=255)
                                    aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa                                              <  1:7880/90‑1 (MQ=255)
                                    aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa                                              <  1:241533/90‑1 (MQ=255)
                                     aCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAACGCCCGCCGATaa                                             >  2:197709/1‑90 (MQ=255)
                                          tGCAATGGCGGGGGAAGGGATAAAAAATATTCTTGAGGGGGGATGCGGGGCTTGACCCGAAAGGGGGGGAAGCGCCTGCCGGTAAGTtgg                                        >  2:328999/1‑89 (MQ=255)
                                          tGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa                                              >  1:180142/1‑84 (MQ=255)
                                          tGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa                                              <  2:180142/84‑1 (MQ=255)
                                             aagggCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGGGGGCCTTGCCCCGAAAGGGGGGAAAGCGCCTGCCGATAAGTTGTGAt                                     <  2:339725/87‑1 (MQ=255)
                                                          gCGATAAAAACTATTCTTGAGGAGTGAT‑‑GGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAc                      <  1:491506/90‑1 (MQ=255)
                                                            gATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAc                      <  1:81119/90‑1 (MQ=255)
                                                                              ggagTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGaa    <  1:149651/90‑1 (MQ=255)
                                                                               gagTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAAt   <  1:283538/90‑1 (MQ=255)
                                                                                agTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAAtt  >  1:336457/1‑90 (MQ=255)
                                                                                  |                                                                                       
GGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAATT  >  NZ_CP009273/4464450‑4464619

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