Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,101,6750GT46.7% 10.5 / 10.3 15N156K (AAC→AAAyiiQYiiQ family protein
Reads supporting (aligned to +/- strand):  ref base G (8/0);  new base T (0/7);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 1.55e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCC  >  NZ_CP009273/4101600‑4101762
                                                                           |                                                                                       
gAAGGTCGCCCCTTTTTCGCCTTTGTCCGCGAAAACATACCGCAGTGCACCTTTTGTCTGGGTGTAATAACGTTTTTTTTTCCCCGCCGt                                                                           <  1:55663/90‑1 (MQ=255)
gAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGt                                                                           >  2:346879/1‑90 (MQ=255)
        gCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCg                                                                   >  2:55663/1‑90 (MQ=255)
                  gCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGttt                                                         >  1:198271/1‑90 (MQ=255)
                                    aTACCGCAGTCCCCCTTCTGTCTGGGTTTAATACCTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGtttt                                       <  1:258596/90‑1 (MQ=255)
                                      aaCGCAGTCCCCTTTTTGTCTGGGTTTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTgg                                     <  2:49392/90‑1 (MQ=255)
                                      aaCGCAGTCCCCCTTCTGTCTCGGTGTAATAACTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTgg                                     <  1:175360/90‑1 (MQ=255)
                                            gTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTa                               >  1:127419/1‑90 (MQ=255)
                                            gTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTa                               >  1:571939/1‑90 (MQ=255)
                                              cccccTTCTGTCTCGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAAt                             <  2:488285/87‑1 (MQ=255)
                                                    tctgtctgGGTGTAATATTTTTTTTTTTCCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTg                       <  1:456987/90‑1 (MQ=255)
                                                    tctgtctgGGTGTAATAACTTTTTTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTg                       <  1:181806/90‑1 (MQ=255)
                                                      tGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAg                     >  2:521962/1‑90 (MQ=255)
                                                           cGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGc                >  1:56960/1‑90 (MQ=255)
                                                             gtgtAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCgg              >  2:577861/1‑90 (MQ=255)
                                                                       gtttttttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                                        <  2:329659/54‑1 (MQ=255)
                                                                       gtttttttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTc                                   <  1:307578/59‑1 (MQ=255)
                                                                       gtttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGttt                                        >  1:329659/1‑54 (MQ=255)
                                                                       gtttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTc                                   >  2:307578/1‑59 (MQ=255)
                                                                       gtttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg    >  2:220515/1‑90 (MQ=255)
                                                                       gtttgtttTTCCCCACCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTg    >  2:522322/1‑90 (MQ=255)
                                                                        ttttttttTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGc   <  1:521962/90‑1 (MQ=255)
                                                                         ttgtttTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGcc  >  1:567247/1‑90 (MQ=255)
                                                                           |                                                                                       
GAAGGTCAGCCCCTTTTCGCCGTTGTCCGCGACAACATAACGCAGTGCACCTTCTGTCTCGGTGTAATAACGTTTGTTTTTCCCCGCCGTTAGTAGCGACTGCAGTTTTTTCTGGCTTTGTGTTTTGGTCATTAATGGGGTGAGTGTGCGGATCACCGCTGCC  >  NZ_CP009273/4101600‑4101762

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