Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,906,0110TG58.8% 0.1 / 34.6 34T275P (ACC→CCC) htpXprotease HtpX
Reads supporting (aligned to +/- strand):  ref base T (2/12);  new base G (20/0);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 1.66e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT  >  NZ_CP009273/1905924‑1906092
                                                                                       |                                                                                 
ggACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCa                                                                                 <  2:311752/90‑1 (MQ=255)
     cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGa                                                                              <  2:674312/88‑1 (MQ=255)
     cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGaac                                                                            >  1:322985/1‑90 (MQ=255)
     cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGaaa                                                                            >  2:321852/1‑89 (MQ=255)
     cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaac                                                                            >  2:293132/1‑90 (MQ=255)
     cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGaaa                                                                            >  2:376475/1‑89 (MQ=255)
     cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGa                                                                              >  1:674312/1‑88 (MQ=255)
     cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCAAGaaa                                                                            >  2:658432/1‑89 (MQ=255)
                   tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGGACAAACCACTTAGGGa                                                              >  2:761789/1‑90 (MQ=255)
                   tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAACAAAACACCGAGCGa                                                              >  2:576384/1‑90 (MQ=255)
                               cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCACTGAGCGGTTTTGACTTTcc                                                  >  1:405546/1‑90 (MQ=255)
                               cTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTTGACTTTcc                                                  >  1:566226/1‑90 (MQ=255)
                                 tCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGt                                                >  1:114452/1‑90 (MQ=255)
                                 acAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGt                                                <  1:186179/89‑1 (MQ=255)
                                  cAGGTATTCACCCGTATGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACGGAGCGATTTCGACTTACCGtt                                               <  2:706372/90‑1 (MQ=255)
                                   aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGGCATGAACAACTCCCTGAGCGGTTTCTACTTACCGTTa                                              >  2:12067/1‑90 (MQ=255)
                                   aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCCTGAAAAACCCCCTGAGCGATTTTGACTTAcccgtg                                              >  2:711796/1‑86 (MQ=255)
                                   aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTa                                              >  1:554526/1‑90 (MQ=255)
                                       aTTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTTGACTTACCGTTAATGc                                          >  1:200578/1‑90 (MQ=255)
                                        ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACtt                                                     <  2:470989/78‑1 (MQ=255)
                                        ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACtt                                                     >  1:470989/1‑78 (MQ=255)
                                            cccGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagaga                                     <  2:200578/90‑1 (MQ=255)
                                               gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGCTGAAAAAATCACTTAGCGGTTTTGACTTTCCGTTAATTCagaggggc                                  >  2:601568/1‑86 (MQ=255)
                                               gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCAAGAACAAAACACCGAGCGGTTTCGACTTACCCCTAATGCagagctgt                                  >  2:186179/1‑86 (MQ=255)
                                                   gCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatca                              <  1:658432/90‑1 (MQ=255)
                                                   gCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCcatca                              <  1:321852/90‑1 (MQ=255)
                                                    cAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGGGAGCcatcat                             >  1:165338/1‑90 (MQ=255)
                                                      gagCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa                                             <  2:66715/72‑1 (MQ=255)
                                                      gagCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa                                             >  1:66715/1‑72 (MQ=255)
                                                          ttCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAACTCACTGAGCGATTTTGACTTTCCCTTAATGCAGAGAGCCATCATGCTggg                       >  2:298827/1‑89 (MQ=255)
                                                                  gTTTATCCAGCGGCGGGGGGGGCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGcttctt               >  2:52684/1‑90 (MQ=255)
                                                                    ttATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGc             <  1:376475/90‑1 (MQ=255)
                                                                               cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt  <  1:761789/90‑1 (MQ=255)
                                                                               cggGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCt  <  1:52684/90‑1 (MQ=255)
                                                                                     ggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa                                             <  2:242504/41‑1 (MQ=255)
                                                                                     ggtCATGAACAACTCACTGAGCGATTTCGACTTACCGTTaa                                             >  1:242504/1‑41 (MQ=255)
                                                                                       |                                                                                 
GGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTTGCGGTTCATAGCT  >  NZ_CP009273/1905924‑1906092

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