Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,906,0110TG62.5% 1.0 / 10.8 16T275P (ACC→CCC) htpXprotease HtpX
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base G (10/0);  total (12/4)
Fisher's exact test for biased strand distribution p-value = 8.24e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.57e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG  >  NZ_CP009273/1905929‑1906093
                                                                                  |                                                                                  
cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaac                                                                             >  1:166898/1‑90 (MQ=255)
 gcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaaca                                                                            >  1:45183/1‑90 (MQ=255)
           aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg                                                                  <  1:272588/90‑1 (MQ=255)
              tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCgg                                                               >  2:10796/1‑89 (MQ=255)
                          cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAAACACTGAGCGCTTTTGTCTTTcc                                                   >  1:277621/1‑90 (MQ=255)
                            tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCCCTGAGCGATTTCTCCTTTCCGt                                                 >  2:330851/1‑90 (MQ=255)
                               ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTAcc                                                   >  2:111661/1‑85 (MQ=255)
                               ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTAcc                                                   <  1:111661/85‑1 (MQ=255)
                               ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGGGCGATTTTGACTTACCGTTaa                                              >  2:258431/1‑90 (MQ=255)
                                    tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGGTTTTGACTTTCCGTTTATGCag                                         >  1:240302/1‑90 (MQ=255)
                                         cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCCCTGAGCGATTTTGCCTTACCGTTAATGCAGAGAGc                                    >  1:578315/1‑90 (MQ=255)
                                          gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAACCACTGAGCGGTTTTTACTTTCCCGTTATGCagaggggc                                   >  1:589708/1‑86 (MQ=255)
                                                      tCAATTCGTTTATCCAGCGGCGGGGGGGGGATTAAAAACCCACCCAGCGATTTTTACTTTCCCTTAAAGCAGAGGGCCCTCATGCTGGtt                       >  1:339097/1‑90 (MQ=255)
                                                             gTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGcttctt                >  2:251805/1‑90 (MQ=255)
                                                               ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc              <  1:445096/90‑1 (MQ=255)
                                                               agaTCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc              <  2:159160/88‑1 (MQ=255)
                                                                           gggggggtCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTg  >  2:220329/1‑90 (MQ=255)
                                                                           ggggggggCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGTTCATAGcct  >  1:159160/1‑88 (MQ=255)
                                                                                  |                                                                                  
CGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG  >  NZ_CP009273/1905929‑1906093

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