Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,906,0070TG63.6% 15.6 / 23.0 22H276P (CAC→CCC) htpXprotease HtpX
Reads supporting (aligned to +/- strand):  ref base T (0/8);  new base G (14/0);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 3.13e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.79e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAT  >  NZ_CP009273/1905919‑1906088
                                                                                        |                                                                                 
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                  <  2:692449/90‑1 (MQ=255)
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                  <  2:166488/90‑1 (MQ=255)
          cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGaaa                                                                        >  1:597674/1‑89 (MQ=255)
                    tATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGa                                                              >  2:839272/1‑90 (MQ=255)
                     aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg                                                             <  1:16087/90‑1 (MQ=255)
                        tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                  <  1:816010/66‑1 (MQ=255)
                        tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACCCACTTAGCGa                                                          >  1:310193/1‑90 (MQ=255)
                                  tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa                                                <  2:501394/90‑1 (MQ=255)
                                    cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGGTTTTTACTTAcc                                              >  1:904995/1‑90 (MQ=255)
                                      tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGATTTTGACTTACCGt                                            >  2:655138/1‑90 (MQ=255)
                                        aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGGTTTTGACTTACCGTTa                                          >  1:633568/1‑90 (MQ=255)
                                              tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag                                    <  2:75154/90‑1 (MQ=255)
                                              tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATAAACAAATCACTTAGCGATTTCTCCTTACCCTTAAAGCag                                    >  2:478078/1‑90 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTTATGAAAAACCCACAGAGGGATTTTGCCTTTCCGTTAATGCagagaggc                              >  2:490563/1‑88 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGcc                              >  2:100681/1‑90 (MQ=255)
                                                           gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc                       <  2:633568/90‑1 (MQ=255)
                                                           gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc                       >  2:859597/1‑90 (MQ=255)
                                                             gCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACACCCTCAGCGATTTTGACTTACCCTTAATGCAGAGAGGCATCATGCTg                     >  1:521744/1‑90 (MQ=255)
                                                             gCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTCAGCGGTTTTGCCTTTCCCTTTATGCAGAGAGCCATCATGCTg                     >  1:825051/1‑90 (MQ=255)
                                                                     tCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGCGCCCACATGCTGGTTGcttc             >  1:372808/1‑90 (MQ=255)
                                                                         ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc         <  2:324991/90‑1 (MQ=255)
                                                                         ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc         <  2:69073/90‑1 (MQ=255)
                                                                         ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc         <  1:409581/90‑1 (MQ=255)
                                                                         ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc         <  2:934275/90‑1 (MQ=255)
                                                                                gcggcgGGGGGGGCATGAACAACTCACTGAGCGGTTTTGACTTTCCCTTTATGCAGAGAGCCATCATGCTGGGTGCTTCTTTCGGTTcct  >  1:336495/1‑88 (MQ=255)
                                                                                        |                                                                                 
GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAT  >  NZ_CP009273/1905919‑1906088

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