Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,781,3850TG57.1% 8.1 / 21.2 21G396G (GGT→GGGenvCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base T (5/4);  new base G (12/0);  total (17/4)
Fisher's exact test for biased strand distribution p-value = 2.11e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.15e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATT  >  NZ_CP009273/3781298‑3781471
                                                                                       |                                                                                      
cTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCg                                                                                      >  2:258037/1‑90 (MQ=255)
   taATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGcc                                                                                   <  2:323742/90‑1 (MQ=255)
         gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCCGGGGCGGccctccct                                                                             >  1:491103/1‑90 (MQ=255)
         gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCAcc                                                                             >  2:385299/1‑89 (MQ=255)
               aCTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTAttt                                                                       <  2:488318/90‑1 (MQ=255)
                   gTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGaa                                                                   >  2:446754/1‑90 (MQ=255)
                       gCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTc                                                               >  2:352632/1‑90 (MQ=255)
                          ttGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCCCTCTATTTCGAAATTcgcc                                                            >  2:399323/1‑90 (MQ=255)
                            ggctCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCCTCACTCTATTTTGAAATTcgccgc                                                          >  1:435537/4‑90 (MQ=255)
                            ggTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCCTCCCTCCATTTCGAAATTcgcccc                                                          >  1:407093/1‑88 (MQ=255)
                                       tCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAggg                                               >  2:194131/1‑89 (MQ=255)
                                       tCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCAGGGGCGGCCTTTACTCGCTTTTGGAATTTGCCGCCAGGGGCAggg                                               >  1:132049/1‑89 (MQ=255)
                                           cagggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc                                           <  2:519521/88‑1 (MQ=255)
                                                        ttGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGTCAATCCACAGCCGt                              >  2:201566/1‑90 (MQ=255)
                                                                         agTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGGTGGGAAGATAAGtt             >  2:309378/1‑90 (MQ=255)
                                                                         agTGGCGGTCAGGGGCGGCCTTCCCTCTATTTTGAAATTCGCCGCCAGGGGCAGGCGGGCAATCCACAGCCGGGGTTGGGGAGATAAGtt             >  2:478113/1‑90 (MQ=255)
                                                                         agTGGCGGTCAGGGGCGGCCTTCCCTCTATTTCGAAATTTGCCGCCCGGGGCAGGGGGGCAATCCCCCGCCGTGGGTGGGGAGATAAggt             >  1:591164/1‑88 (MQ=255)
                                                                            ggCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTg          >  2:53128/1‑90 (MQ=255)
                                                                                    gggTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCAtt  >  2:151481/1‑90 (MQ=255)
                                                                                      gTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc                                           >  1:529247/1‑47 (MQ=255)
                                                                                      gTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc                                           <  2:529247/47‑1 (MQ=255)
                                                                                       |                                                                                      
CTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATT  >  NZ_CP009273/3781298‑3781471

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