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

GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTT  >  NZ_CP009273/1905919‑1906079
                                                                                        |                                                                        
gATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                         <  1:330348/90‑1 (MQ=255)
         gcgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGaa                                                                <  2:717715/90‑1 (MQ=255)
          cgcTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                         <  2:53757/80‑1 (MQ=255)
                     acagaCTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAg                                                    <  1:529670/86‑1 (MQ=255)
                        tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGa                                                 >  2:21822/1‑90 (MQ=255)
                        tACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGAACAACTCACTTAGCGa                                                 >  1:229671/1‑90 (MQ=255)
                                        aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCgggtg                                                                         <  1:173221/50‑1 (MQ=255)
                                        aGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGTTTAACAAATCACTGAGCGATTTCTACTTTCCCTTa                                 >  2:34025/1‑90 (MQ=255)
                                             ttCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATGAACAAACCACTGAGCGATTTTGACTTATCCCTAATGCa                            >  1:135068/1‑90 (MQ=255)
                                              tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCag                           <  2:532700/90‑1 (MQ=255)
                                              tCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCCCTGAACGATTTTGACTTCCCGGTAATGCag                           >  1:749392/1‑90 (MQ=255)
                                                aCCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCagag                         <  1:34025/90‑1 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGTGTTGAACCACTCACCGAGCGAGTTCGAATTACCCTTGATTCACGGAGcc                     >  2:438625/1‑90 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGGATTAACAACTCCCTGAGCGGTTTCGACTTACCCTTTATGCAGAGAGcc                     >  2:71848/1‑90 (MQ=255)
                                                    gTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGGGGGGCATGAAAAAATCACTGAGCGATTTTGACTTACCGTTAATGCAGAGAGcc                     >  1:532700/1‑90 (MQ=255)
                                                                 cAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTg        <  1:21822/90‑1 (MQ=255)
                                                                       gTTTATCCAGCGGCGGGGGGGGGATGAACAACTCACTGAGCGATTTTTACTTTCCCTTAATGGAGGGAGCCATCCTTTTGGGTGcttctt  >  2:42373/1‑90 (MQ=255)
                                                                                        |                                                                        
GATCAGGACGCGCTTTTTAGTATTTACTTCATATTACTTCAGGTATTCACCCGTACGCAGAGCTTCAATTCGTTTATCCAGCGGCGGGTGGGTCATGAACAACTCACTGAGCGATTTCGACTTACCGTTAATGCAGAGAGCCATCATGCTGGTTGCTTCTT  >  NZ_CP009273/1905919‑1906079

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