Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,781,3850TG54.5% 16.4 / 27.3 22G396G (GGT→GGGenvCmurein hydrolase activator EnvC
Reads supporting (aligned to +/- strand):  ref base T (3/7);  new base G (12/0);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 7.04e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-04
Rejected as consensus: Frequency below/above cutoff threshold.

CTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTT  >  NZ_CP009273/3781298‑3781472
                                                                                       |                                                                                       
cTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCg                                                                                       >  1:141647/1‑90 (MQ=255)
        agagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCAc                                                                               <  1:276818/90‑1 (MQ=255)
         gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCAct                                                                              >  1:182357/1‑90 (MQ=255)
         gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCAct                                                                              >  1:392537/1‑90 (MQ=255)
         gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGTAGGGGCGGCCTTCCct                                                                              >  1:605088/1‑90 (MQ=255)
                cTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTc                                                                       <  1:630103/90‑1 (MQ=255)
                            ggTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTcgccgc                                                           >  1:79918/1‑90 (MQ=255)
                                           gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc                                            <  1:103869/90‑1 (MQ=255)
                                           gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc                                            <  2:909038/90‑1 (MQ=255)
                                           gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc                                            <  2:354483/90‑1 (MQ=255)
                                                       aTTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCg                                <  1:874988/90‑1 (MQ=255)
                                                        ttGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCCCTCTTTTTTGAAAATTGCCCCCCGGGGCAGGCGGGCAATCCCCAGCCgg                               >  2:953650/1‑89 (MQ=255)
                                                            aCTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCGTTACTCTATTTTGGAAATGGCCGGCAAGGGCAGGGGGGCAATACAGAGCCgggggt                           >  1:944273/1‑88 (MQ=255)
                                                                gtggGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGTCAATCCACAGCCGGGGGTgggg                       >  1:824456/1‑89 (MQ=255)
                                                                    gcagcaGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGGCAATCCACAGCCGTGGGTGGGGAGAt                   >  2:350658/1‑90 (MQ=255)
                                                                    gcagcaGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCCCCAGGGGCAGGGGGGCAATCCACAGCCGTGGGTGGGGAGAt                   >  1:664686/1‑90 (MQ=255)
                                                                       gcagTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAg                >  1:158875/1‑90 (MQ=255)
                                                                         agTGGCGGTCAGGGTCGGCCTTCACTCTATTTTGAAATTCGCCGCCAGGGGCAGGCGGGCAATCCCCAGCCGTGGGTGGGGAGATAAggt              >  2:861879/1‑88 (MQ=255)
                                                                         agTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGAAATTCGCCGCCCGGGGCAGGGGGGCAAACCCCAGCCGTGGTTGGGGAGATAAGtt              >  2:633426/1‑90 (MQ=255)
                                                                               ggTCAGGGTCGGCCTTCACTCTATTTCGAAAATCGCCGCCAGGGGCAGGCGGGCAATCCACAGCCGTGGGTGGGAAGATAAGTTTTGttt        >  1:615537/1‑90 (MQ=255)
                                                                               ggTCAGGGGCGGCCCTCACTCTTTTTTGAAATTCGCCGCCAGGGGCAGGGGGGCAATCCCCAGCCGGGGGTGGGGAGATAAGATTTGttt        >  2:592020/1‑90 (MQ=255)
                                                                                     ggTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCAttt  <  2:182357/90‑1 (MQ=255)
                                                                                       |                                                                                       
CTATAATCAGAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTTCCATTT  >  NZ_CP009273/3781298‑3781472

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