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

CTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAA  >  NZ_CP009273/4464443‑4464617
                                                                                         |                                                                                     
cTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGt                                                                                       <  1:688387/90‑1 (MQ=255)
     aaGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGc                                                                                  <  2:581043/90‑1 (MQ=255)
      aGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCg                                                                                 <  2:826827/90‑1 (MQ=255)
      aGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCg                                                                                 <  2:105830/90‑1 (MQ=255)
                ggggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGGGATAAAAAATATTTTTGAGGAGTGGTGGGGGGCTTGAcc                                                                       >  1:433390/1‑90 (MQ=255)
                 gggCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATTTTCTTGAGGGGAGATGCGGGGCTTGAccc                                                                      >  1:739710/1‑90 (MQ=255)
                   gCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGAGGGATGCGGGGCTTGACCCGa                                                                    >  2:367895/1‑90 (MQ=255)
                     aCTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGGTGCGGGCCTTGACCCGaaa                                                                  >  1:114716/1‑90 (MQ=255)
                      ctctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAAAATTTTTTAGGAGGGAGGCGGGGCTTTGCCCGGAAg                                                                 >  2:940460/1‑90 (MQ=255)
                      ctctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTTTTGGGGGGGGGTGCGGGGCTTGGGCCGGAAg                                                                 >  2:866439/1‑90 (MQ=255)
                       tctATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAgg                                                                >  1:5474/1‑90 (MQ=255)
                          ataGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTGGGGGTGATGCGGGGCTTGACCCGAAAGGCgg                                                             >  2:673041/1‑90 (MQ=255)
                           taGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGGGGGATGCGGGGCCTGACCCGAAAAGCggg                                                            >  2:112273/1‑90 (MQ=255)
                                caacTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTAGGGGGGATGCGGGGCTTGACCCGAAAGGCGGGGAAga                                                       >  2:506535/1‑89 (MQ=255)
                                caacTCTGTGCAACAGCTGCAATGGCGGGGGGAACGATAAAAATATTCTTTGGGAGAGATGCGGGCCCTGACCCGAAAGGGGGGGAAacg                                                      >  2:878104/1‑90 (MQ=255)
                                caacTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAgc                                                       >  1:932929/1‑90 (MQ=255)
                                caacTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgc                                                       >  1:338899/1‑90 (MQ=255)
                                caacTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgc                                                       >  1:330633/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAGTGGCGGGGGGAGCGATAAAAACTATTTTTGAGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc                                                     >  1:346857/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAATATTCTTTAGGGGGGATGGGGGCCTTTGCCCGAAAGGGGGGAAAgcgc                                                     >  1:478634/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTTTTTTGGGGGGATGGGGGCCTTGACCCGAAAGGGGGGGAAacgc                                                     >  2:596280/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTTTTGAGGGGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc                                                     >  1:733172/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                                                     >  2:337584/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                                                     >  2:821257/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTTAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                                                     >  2:386825/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGTGGAGGGATGCGGGGCCTGGCCCGAAAGGCGGGAAAacgc                                                     >  2:498765/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGGCTTGACCCGAAAGGCGGGGAAacgc                                                     >  2:826993/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                                                     >  2:241455/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTTAGGGGGGATGCGGGCCTTGACCCGAAAAGCGGGGAAgcgc                                                     >  1:885997/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc                                                     >  1:787916/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTCTTTAGGAGGGATGCGGGGCTTGACCCCAAAGGCGGGGAAacgc                                                     >  2:238521/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGGAGCGAAAAAAACTATTTTTTAGGAGGGATGCGGGCCTTGACCCGAAAAGCGGGaaaaagc                                                     >  2:640176/1‑86 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGGAAacgc                                                     >  1:105964/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc                                                     >  1:662953/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc                                                     >  1:171789/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGAAAgcgc                                                     >  1:374515/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGCGGGGAAgcgc                                                     >  1:785691/1‑90 (MQ=255)
                                  acTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCTTGACCCGAAAAGCGGGGAAacgc                                                     >  2:191672/1‑90 (MQ=255)
                                   ctctGTGCAACAGCTGCAATGGGGGGGGAAGCGAAAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGcc                                                    <  2:895041/90‑1 (MQ=255)
                                   ctctGTGCAACAGCTGAAAGGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGcc                                                    <  2:76630/90‑1 (MQ=255)
                                    tctGTGCAACAGCTGCAATGGCGGGGGGAGGGATAAAAAAAATTTTTTAGGGGGGATGCGGGGCTTGACCCGAAAGGCGGGGAAAAGCCt                                                   >  2:901629/1‑90 (MQ=255)
                                    tctGTGCAACAGCTGCAATGGCGGGGGGAGCGATAAAAAAAATTTTTTAGGAGGGATGCGGGCCCTGAGCCGGAAGGCGGGGAAACGccc                                                   >  2:584614/1‑89 (MQ=255)
                                       gtgCAACAGCTGCAATGGCGGGGGGAGCGGTAAAAAATATTTTTGAGGGGGGATGCGGGCCTTGACCCGAAAGACGGGAAAACGCCTGcc                                                >  2:156115/1‑90 (MQ=255)
                                       gtgCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGcc                                                >  1:59391/1‑90 (MQ=255)
                                          cAACAGCTGCAATGGCGGGGGGAGCGATAAAAACTATTCTTGAGGAGGGATGCGGGCCTTGACCCCAAAGGCGGGAAAACGCCCGCCGAt                                             >  1:951100/1‑90 (MQ=255)
                                           aaCAGGTGCAATGGCGGGGGAAGCGATAAAAACTTTTCTTGGGGAGTGAGGGGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATa                                            <  2:721886/90‑1 (MQ=255)
                                           aaCAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGACGGGAAAGCGCCTGCCGATa                                            <  1:156115/90‑1 (MQ=255)
                                              aGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGt                                         <  2:338899/90‑1 (MQ=255)
                                               gCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGtt                                        <  1:7609/90‑1 (MQ=255)
                                               gCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGtt                                        <  2:82231/90‑1 (MQ=255)
                                                 tGCAATGGCGGGGGGAGGGGTAAAAACTATTTTTGAGGAGGGGTGCGGGGCCTGACCCGAAAGGCGGGGGAGCGCCCGCCGATAAggttt                                      >  2:672073/1‑86 (MQ=255)
                                                 tGCAATGGCGGGGGGAGCGATAAAAAATATTTTTGAGGAGGGATGCGGGCCTTGACCCGAAAGGGGGGAAAAGGCCTGCCGATAAGTtgg                                      >  2:516367/1‑89 (MQ=255)
                                                 tGCAATGGCGGGGGGAGCGATAAAAAATATTCTTGAGGAGGGATGCGGGCCCTGACCCGAAAGGCGGGAAAGCGCCTGCCGCTAAGTtgt                                      >  2:918160/1‑90 (MQ=255)
                                                     aTGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATa                                  <  1:669892/90‑1 (MQ=255)
                                                                 gCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTc                      >  2:535648/1‑90 (MQ=255)
                                                                   gATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAc                    <  1:579092/90‑1 (MQ=255)
                                                                   gATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCAc                    <  1:901629/90‑1 (MQ=255)
                                                                         aaCTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTa              <  1:830253/90‑1 (MQ=255)
                                                                         aaCTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTa              <  1:506535/90‑1 (MQ=255)
                                                                            tgttCTTGGGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAgg           <  2:341294/88‑1 (MQ=255)
                                                                                  tGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCg     >  2:791996/1‑90 (MQ=255)
                                                                                   gaggagTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCgg    >  2:730930/1‑90 (MQ=255)
                                                                                     ggagTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGaa  <  2:662953/90‑1 (MQ=255)
                                                                                         |                                                                                     
CTGCGAAGGCACGAGAGGGGCACTCTATAGAACAACTCTGTGCAACAGCTGCAATGGCGGGGGAAGCGATAAAAACTATTCTTGAGGAGTGATGCGGGCCTTGACCCGAAAGGCGGGAAAGCGCCTGCCGATAAGTTGTGATAAGACAATAATTCACGCGTTAAGGCTAGCGGAA  >  NZ_CP009273/4464443‑4464617

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: