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

GAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTT  >  NZ_CP009273/3781307‑3781466
                                                                              |                                                                                 
gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTGAGGGGCGGGCTTcccc                                                                        >  1:338261/1‑87 (MQ=255)
gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGGCTTCAcc                                                                        >  1:121141/1‑89 (MQ=255)
gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCAct                                                                        >  2:91144/1‑90 (MQ=255)
gagCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCCGGGGCGGCCTTCAct                                                                        >  1:286536/1‑90 (MQ=255)
   cgcACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTa                                                                     >  1:126364/1‑90 (MQ=255)
    gcACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCAGGGGCGGCCTTCActcttt                                                                    >  1:413438/1‑88 (MQ=255)
     cACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTAtt                                                                   >  2:155656/1‑90 (MQ=255)
          gTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTCGaa                                                              >  2:130303/1‑90 (MQ=255)
              gCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTTGAAATTc                                                          >  2:58041/1‑90 (MQ=255)
                   ggTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGGCAGGGGCGGGCCTCCCTCTATTTTGAAATTcgccgc                                                     >  1:456700/1‑90 (MQ=255)
                                  gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc                                      <  1:306054/90‑1 (MQ=255)
                                  gcgGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTc                                      <  2:338261/90‑1 (MQ=255)
                                     ggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTcg                                                         <  1:131026/68‑1 (MQ=255)
                                     ggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTcg                                                         >  2:131026/1‑68 (MQ=255)
                                     ggCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAAATTGCCGCCCGGGTCAGGCGGTCAAt                                   >  1:43219/1‑90 (MQ=255)
                                        cAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCca                                <  2:283706/90‑1 (MQ=255)
                                                      ggtggGCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTCGCCGCCAGGGGCAGGGGGGCAATCCACAGCCGTGGTTggg                  >  1:116572/1‑90 (MQ=255)
                                                         gggCAGCAGTGGCGGTCAGGGGCGGCCTTCACTCTATTTTGAAATTCGCCGCCAGGGGCAGGCGGTCAATCCACAGCCGTGGGTGGGGAg               >  2:341488/1‑90 (MQ=255)
                                                              gcagTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAg          >  2:126439/1‑90 (MQ=255)
                                                                agTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGTCGATCCCCAGCCGTGGTTGGGGAGATAAGtt        >  2:396660/1‑90 (MQ=255)
                                                                agTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGGCAGGCGGGCAATCCACAGCCGGGGTTGGGGAGATAAGtt        >  2:371527/1‑90 (MQ=255)
                                                                      ggTCAGGGGCGGCCTTCACTCTATTTTGAAATTTGCCGCCAGGGGCAGGGGGGCAATCCACAGCCGGGGGTGGGGAGATAAGGTtttttt  >  1:410480/1‑90 (MQ=255)
                                                                              |                                                                                 
GAGCGCACTGGTGAGCGTTGGTTCGCAGGTTCGCGCGGGCCAGCCGATTGCACTGGTGGGCAGCAGTGGCGGTCAGGGTCGGCCTTCACTCTATTTCGAAATTCGCCGCCAGGGTCAGGCGGTCAATCCACAGCCGTGGTTGGGAAGATAAGTTTTGTTT  >  NZ_CP009273/3781307‑3781466

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: