Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,464,5060AG61.1% 2.4 / 25.7 37E189G (GAA→GGA) bdcRTetR/AcrR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (8/6);  new base G (22/0);  total (31/6)
Fisher's exact test for biased strand distribution p-value = 1.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.13e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAG  >  NZ_CP009273/4464422‑4464587
                                                                                    |                                                                                  
gTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGa                                                                               >  1:59729/1‑90 (MQ=255)
gTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGa                                                                               >  1:367644/1‑90 (MQ=255)
     gTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATaaaa                                                                          >  2:474066/1‑90 (MQ=255)
          gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGGGGTAAAAAatat                                                                     >  1:499138/1‑90 (MQ=255)
          gTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAt                                                                     <  2:147222/90‑1 (MQ=255)
             tGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGGGCGATGAAAACAAttgt                                                                  >  1:560773/1‑88 (MQ=255)
             tGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAAAACAACTCTGTGCCAGAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCt                                                                  <  1:344390/90‑1 (MQ=255)
                   ttCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTgtgga                                                            >  1:487966/1‑90 (MQ=255)
                                      gggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGAccc                                         >  2:259024/1‑90 (MQ=255)
                                      gggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTTTTGAGGGGGGATGCGGGCCTTGAccc                                         >  2:250587/1‑90 (MQ=255)
                                            tctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTGGGGGGGGGGCGGGGCTTGGCCCGAAAgg                                  >  1:269420/1‑90 (MQ=255)
                                                  gAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGTGAAAAACTTTTCTTTAGGAGCGGTTCGGGCCGTGAACCAAAAAGCGGGaa                            >  2:259134/1‑90 (MQ=255)
                                                  gAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGGGGGATGCGGGCCTTGACCCGAAAGGCgggga                            >  2:531728/1‑88 (MQ=255)
                                                     caacTCTGTGCAACAGCTGCAATGGCGGGGGGAGGGATAAAAAATATTCTTTAGGAGGGATGCGGGCCTTGTCCCGAAAAGGGGGGAAgg                         >  1:507462/1‑89 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAATATTTTTGAGGGGGGAAGCGGGGCTTGACCCCAAAGGGGGGGAAacgc                       >  1:513703/1‑90 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCCAAAAGCGGGGAAgcgc                       >  1:41638/1‑90 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                       >  1:108869/1‑90 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGAGAGATGCGGGCCTTTACCCGAAAGGCGGGGAAacgc                       >  1:594081/1‑90 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGTGATGCGGGCCTTTACCCGAAAGGCGGGGAAacgc                       >  2:459735/1‑90 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc                       >  1:198618/1‑90 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                       >  1:128628/1‑90 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAACATTCTTTGGGAGGGATGCGGGCCGTGACCCGAAAGGCGGGGAAgggc                       >  1:574077/1‑90 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAt                                                                     <  1:11836/45‑1 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTAt                                                                     >  2:11836/1‑45 (MQ=255)
                                                       acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc                       >  2:149095/1‑90 (MQ=255)
                                                        ctctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGcc                      <  1:250587/90‑1 (MQ=255)
                                                         tctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGGGGGAAAACGCCt                     >  2:146314/1‑90 (MQ=255)
                                                         tctGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCt                     >  2:53911/1‑90 (MQ=255)
                                                           tgtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGc                   <  2:499138/90‑1 (MQ=255)
                                                            gtgCAACAGCTGCAATGGCGGGGGGAGGGATAAAAAAAATTTTTGAGGGGGGATGCGGGCCTTGGCCCGAAAGGGGGGGAAAGGGccgcc                   >  1:71626/1‑87 (MQ=255)
                                                            gtgCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTTTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGGGGGGAAGCGCCTGcc                  >  2:269501/1‑90 (MQ=255)
                                                            gtgCAACAGCTGCAATGGCGGGGGCAGGGATAAAAACTAATCTTGAGGGGCGATGCGGGGCTTGACCCGAAAGGCGGGAAAGcgcccgcc                  >  2:227659/1‑90 (MQ=255)
                                                            gtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGcc                  >  2:31810/1‑90 (MQ=255)
                                                                    gCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGGCTTGGCCCGAAAGGGGGGGAAAAGCCTGCCGGTAAGtt          >  1:289548/1‑90 (MQ=255)
                                                                    gCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAggt          >  1:171456/1‑88 (MQ=255)
                                                                      tGCAATGGCGGGGGAAGCGATAAAAACTATTCTTTAGGGGTGATGCGGGCCTTTACCCGAAAGGCGGGAAAGCGCCCGCCCATAAATtgt        >  2:268151/1‑90 (MQ=255)
                                                                            ggCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGCCTTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAg  <  2:560773/90‑1 (MQ=255)
                                                                                    |                                                                                  
GTACCGTGATGTCTGGGCTTTCTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAG  >  NZ_CP009273/4464422‑4464587

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: