Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,906,0070TG63.3% 11.4 / 19.5 30H276P (CAC→CCC) htpXprotease HtpX
Reads supporting (aligned to +/- strand):  ref base T (5/6);  new base G (19/0);  total (24/6)
Fisher's exact test for biased strand distribution p-value = 7.78e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.58e-02
Rejected as consensus: Frequency below/above cutoff threshold.

GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTG  >  NZ_CP009273/1905919‑1906093
                                                                                        |                                                                                      
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  1:1152084/90‑1 (MQ=255)
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  1:1694281/90‑1 (MQ=255)
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  2:41616/90‑1 (MQ=255)
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  2:492004/90‑1 (MQ=255)
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  1:1909374/90‑1 (MQ=255)
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  2:491167/90‑1 (MQ=255)
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  2:449978/90‑1 (MQ=255)
   cAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtgggt                                                                                    <  2:787006/90‑1 (MQ=255)
                     aTTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg                                                                  <  1:881874/90‑1 (MQ=255)
                      tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGc                                                                 >  1:308782/1‑90 (MQ=255)
                      tttACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGc                                                                 >  1:216942/1‑90 (MQ=255)
                             cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa                                                                                  <  2:1027348/66‑1 (MQ=255)
                             cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCa                                                                                  >  1:1027348/1‑66 (MQ=255)
                             cATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGGttttg                                                          >  1:1746339/1‑88 (MQ=255)
                                      tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATTAACAACTCACTGGGGGATTTTGTCATACCGt                                                 >  1:410056/1‑90 (MQ=255)
                                      tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAAAAACTCACTGAGCGATTTTTACTTACCgg                                                 >  2:1603407/1‑89 (MQ=255)
                                      tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACACACTGAGCGATTTTGACTTACCGt                                                 >  2:1677240/1‑90 (MQ=255)
                                      tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAAATCACTGAGCGATTTTGACTTTCCgg                                                 >  2:1280105/1‑89 (MQ=255)
                                        aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATTGACAAACCACCGAGCGATTTTGTCTTACCCTTa                                               >  2:1285680/1‑90 (MQ=255)
                                         ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa                                                                                  <  1:568705/54‑1 (MQ=255)
                                         ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa                                              >  1:379970/1‑90 (MQ=255)
                                         ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACACACTGAGCGATTTCGACTTACCGTTaa                                              >  1:1756102/1‑90 (MQ=255)
                                         ggTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCa                                                                                  >  2:568705/1‑54 (MQ=255)
                                              tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag                                         <  1:341129/90‑1 (MQ=255)
                                                   cGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACCCACTGAGCGGTTTTGACTTACCGGTTAAGGAGAGAGc                                    >  2:213966/1‑90 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                                       <  1:1123378/38‑1 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGcc                                   >  1:1800021/1‑90 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAAAAAATCACTGAGCGATTTTGACTTTCCGTTAATGCagagaggc                                   >  1:1664356/1‑88 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGGTTTCGACTTACCGTTAATGCAGAGAGcc                                   >  1:738349/1‑90 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGcc                                   >  1:626062/1‑90 (MQ=255)
                                                           gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGc                            >  1:673186/1‑90 (MQ=255)
                                                           gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCCTCATGc                            >  1:240072/1‑90 (MQ=255)
                                                                       gTTTATCCAGCGGCGGGGGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGcttctt                >  1:54425/1‑90 (MQ=255)
                                                                       gTTTATCCAGCGGCGGGGGGGAGATGAACAACCCACTGAGCGGTTTCGAATTACCGTTTATGCAAAGAGCCATTCTTCTTGGTGcttctt                >  1:1278500/1‑90 (MQ=255)
                                                                         ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc              <  2:1746339/90‑1 (MQ=255)
                                                                           aTCCAGCGGCGGGTGGGTCATGAACAACTCCCTGAGCGATTTTGCCTTACCGTTAATGCAGAGGGCCATCAAGCTGGTTTCTTCTTTCgg            >  1:29148/1‑90 (MQ=255)
                                                                             ccAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGtt          <  2:54425/90‑1 (MQ=255)
                                                                                gcggcgGGGGGGGCATGAACAACTCACTTAGCGATTTTGACTTTCCCGTAATGCAGAGAGCCACCATGCTGGTTGCTTTTTTCGGGTcct       >  2:1660139/1‑88 (MQ=255)
                                                                                gcggcgGGGGGGGCATGAACAACTCACTGAGCGGTTTTGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCAt       >  2:1035207/1‑90 (MQ=255)
                                                                                     gggtggggCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCTg  >  2:746241/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 < 8 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: