Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,464,5180CA32.4% 49.2 / 10.5 35T193N (ACT→AAT) bdcRTetR/AcrR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (7/16);  new base A (11/0);  total (19/16)
Fisher's exact test for biased strand distribution p-value = 1.27e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAC  >  NZ_CP009273/4464430‑4464599
                                                                                        |                                                                                 
aTGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACt                                                                                  <  1:225887/90‑1 (MQ=255)
  gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGGGATAAAAAatat                                                                                >  2:288063/1‑90 (MQ=255)
  gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGGTAAAAAatat                                                                                >  1:297731/1‑90 (MQ=255)
  gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAt                                                                                <  2:173156/90‑1 (MQ=255)
  gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAt                                                                                <  1:254553/90‑1 (MQ=255)
  gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAt                                                                                <  1:276478/90‑1 (MQ=255)
           ttCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTgag                                                                         >  1:447787/1‑88 (MQ=255)
           ttCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTgag                                                                         <  2:447787/88‑1 (MQ=255)
            tCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAAAATTCTTtagggg                                                                      >  2:530892/1‑88 (MQ=255)
                 gAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATg                                                                 >  1:87781/1‑90 (MQ=255)
                          agagGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTg                                                        <  1:530892/90‑1 (MQ=255)
                           gagGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGa                                                       <  1:240313/90‑1 (MQ=255)
                                  aCTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTCTTTAGGAGGGATGCGGGGCTTGACCCGaaa                                                >  2:460641/1‑90 (MQ=255)
                                  aCTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGaaa                                                <  1:144179/90‑1 (MQ=255)
                                       ataGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAGgggg                                           >  1:346247/1‑90 (MQ=255)
                                          gAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGGGGGATGCGGGGCTTGACCCCAAAAGCGGGaa                                        >  2:302162/1‑90 (MQ=255)
                                             caacTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTAGGGGGGATGCGGGCCTTGACCCGAGAGGCGGGgaaag                                     >  2:453502/1‑88 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGGAGGGATAAAAAATATTTTTTAGGAGGGATGCGGGCCTTGACCCGAAAAGCGGGaaaaagc                                   >  2:290784/1‑86 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                                   >  1:585043/1‑90 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGTGATGCGGGCCTTGGCCCGAAAGGCGGGGAAacgc                                   >  2:22784/1‑90 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGGCTTGACCCGAAAGGCGGGGAAacgc                                   >  2:509429/1‑90 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGCAGTGACAAAAAGCATTTTTTAGGAGGGCGACAGGCCCTTCCCCGAAATGGGGGAAAacgc                                   >  1:511222/1‑90 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACt                                                                                  <  1:51807/43‑1 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACt                                                                                  >  2:51807/1‑43 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc                                   >  1:76354/1‑90 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc                                   >  1:192220/1‑90 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTTTTGGGGAGGGATGCGGGCCTTGACCCGAAAGGGGGGGAAgggc                                   >  1:438003/1‑90 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTTTTGAGGGGGGATGCGGGCCTTGACCCGAAAAGCGGGGAAacgc                                   >  2:212137/1‑90 (MQ=255)
                                               acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTCTTGTGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAacgc                                   >  2:375870/1‑90 (MQ=255)
                                                ctctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGcc                                  <  2:76354/90‑1 (MQ=255)
                                                    gtgCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTTGGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAAAGCCCGcc                              >  2:424220/1‑90 (MQ=255)
                                                        aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa                          <  1:22784/90‑1 (MQ=255)
                                                        aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa                          <  2:87781/90‑1 (MQ=255)
                                                                 aaTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGAt                 <  1:290784/90‑1 (MQ=255)
                                                                    ggCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAg              <  2:212922/90‑1 (MQ=255)
                                                                       gggggAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACa           <  1:218866/90‑1 (MQ=255)
                                                                                gATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAc  <  2:192220/90‑1 (MQ=255)
                                                                                        |                                                                                 
ATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAC  >  NZ_CP009273/4464430‑4464599

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