Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,906,0070TG57.1% 19.2 / 24.5 21H276P (CAC→CCC) htpXprotease HtpX
Reads supporting (aligned to +/- strand):  ref base T (2/7);  new base G (12/0);  total (14/7)
Fisher's exact test for biased strand distribution p-value = 3.10e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.95e-04
Rejected as consensus: Frequency below/above cutoff threshold.

GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG  >  NZ_CP009273/1905919‑1906093
                                                                                        |                                                                                      
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  2:102504/90‑1 (MQ=255)
          cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaaa                                                                             >  2:255588/1‑89 (MQ=255)
                    tATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGa                                                                   >  2:141093/1‑90 (MQ=255)
                        tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGa                                                               >  2:293342/1‑90 (MQ=255)
                             cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       >  2:307858/1‑61 (MQ=255)
                             cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  1:307858/61‑1 (MQ=255)
                             cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCg                                                          >  2:119967/1‑90 (MQ=255)
                               tatTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGTATGAAAAAAACACTGAGCGAgtttgtc                                                        >  1:192167/1‑86 (MQ=255)
                                 ttACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACtt                                                      <  2:844030/90‑1 (MQ=255)
                                  tACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTa                                                     <  2:950139/90‑1 (MQ=255)
                                    cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATg                                                                                <  2:143576/61‑1 (MQ=255)
                                    cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATg                                                                                >  1:143576/1‑61 (MQ=255)
                                        aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGATTTTGACTTACCGTTa                                               >  2:767096/1‑90 (MQ=255)
                                         ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa                                              >  2:227009/1‑90 (MQ=255)
                                              tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGGTTTCGACTTTCCGTTAATGCag                                         >  2:578486/1‑90 (MQ=255)
                                              tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACCGAGCGGTTTTGGCTTACCGTTAATGGag                                         >  1:675822/1‑90 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  1:129568/38‑1 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       >  2:129568/1‑38 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCAGAAAAACACACTGCGCGATTTTGACTTACCGCTAAAGCagggaggc                                   >  1:801052/1‑88 (MQ=255)
                                                         cAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatcat                              >  2:228053/1‑90 (MQ=255)
                                                           gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGGTTTCGACTTACCGTTTATGCAGAGAGCCATCATGc                            >  1:134547/1‑90 (MQ=255)
                                                                        tttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAAAGCAGAGAGCCATCATGCTGGGTGcttcttt               >  2:311080/1‑89 (MQ=255)
                                                                         ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc              <  1:227009/90‑1 (MQ=255)
                                                                         ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc              <  1:311080/90‑1 (MQ=255)
                                                                          tATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCg             <  1:676425/90‑1 (MQ=255)
                                                                                    cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt   <  1:578486/90‑1 (MQ=255)
                                                                                     gggtggggCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTTCGGGTCATAGCTg  >  1:24853/1‑90 (MQ=255)
                                                                                        |                                                                                      
GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG  >  NZ_CP009273/1905919‑1906093

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