Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,906,0110TG59.1% 10.4 / 24.3 22T275P (ACC→CCC) htpXprotease HtpX
Reads supporting (aligned to +/- strand):  ref base T (1/8);  new base G (13/0);  total (14/8)
Fisher's exact test for biased strand distribution p-value = 2.81e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.24e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAT  >  NZ_CP009273/1905929‑1906088
                                                                                  |                                                                             
cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGaaa                                                                        >  1:597674/1‑89 (MQ=255)
          tATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGa                                                              >  2:839272/1‑90 (MQ=255)
           aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg                                                             <  1:16087/90‑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:934275/90‑1 (MQ=255)
                                                               ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc         <  2:69073/90‑1 (MQ=255)
                                                               ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc         <  2:324991/90‑1 (MQ=255)
                                                               ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc         <  1:409581/90‑1 (MQ=255)
                                                                      gcggcgGGGGGGGCATGAACAACTCACTGAGCGGTTTTGACTTTCCCTTTATGCAGAGAGCCATCATGCTGGGTGCTTCTTTCGGTTcct  >  1:336495/1‑88 (MQ=255)
                                                                                  |                                                                             
CGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAT  >  NZ_CP009273/1905929‑1906088

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