Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273699,6270TG56.2% 12.8 / 20.5 16G76G (GGT→GGG) nagEPTS N‑acetyl glucosamine transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (0/7);  new base G (9/0);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 8.74e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.80e-04
Rejected as consensus: Frequency below/above cutoff threshold.

TTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAACCCAGAAATTAACATGGGTGTACTGGCGGG  >  NZ_CP009273/699565‑699701
                                                              |                                                                          
ttCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCgg                                                                        <  1:542500/67‑1 (MQ=255)
ttCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGGGCgg                                                                        >  2:542500/1‑67 (MQ=255)
    ccATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGa                                             <  2:581830/90‑1 (MQ=255)
      aTCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGGGCGGGAGGGTACTTTTTGTTTAACAAAACGATg                                           >  2:511279/1‑90 (MQ=255)
              ggCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGGGGGGGAGGGTACTTTTTGTTTAACAAAACGGTGGGGGCCAt                                   >  2:14831/1‑90 (MQ=255)
                    cAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGc                                                                          <  2:128318/45‑1 (MQ=255)
                    cAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGGGc                                                                          >  1:128318/1‑45 (MQ=255)
                         ggTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGGGCGGGAGGTTACTTTTTGTTTAACAAAACGATGGGGACCATCAACCCCGaaa                        >  2:217112/1‑90 (MQ=255)
                         ggTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGGGCGGGAGGGTACTTTTTGTTAACCAAAACGATGGGGACCCTCAACCCCGaaa                        >  2:90468/1‑90 (MQ=255)
                              aaaGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAACCCAGAAATTAAc                   <  1:90468/90‑1 (MQ=255)
                                 gACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAACCCAGAAATTAACATg                <  2:157728/90‑1 (MQ=255)
                                      cgcTGGTGCGGCGGCGCTGGCGGGGGCGGGAGGGTTCTTTTTGTTTACCAAAACGATGGGGGCCCTCAACCCACAAAATAAAAAGGgggt           >  2:281411/1‑90 (MQ=255)
                                        cTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAACCCAGAAATTAACATGGGTGTAc         <  2:590925/90‑1 (MQ=255)
                                          ggTGCGGCGGCGCTGGCGGGGGCGGGGGGGTTCTTTTTGTTTAACAAAACGGTGGGGGCCCTCAACCCCAAAAATAAAAAGGGGGTACTg       >  1:226273/1‑90 (MQ=255)
                                              cggcggcgCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAACCCAGAAATTAACATGGGTGTACTGGCgg   <  1:328004/90‑1 (MQ=255)
                                               ggcggcgCTGGCGGGGGCGGGAGGGTTCTTTTTGTTTACCCAAACGGGGGGGGCCCTCAACCCCGAAAATAAAAAGGGGGTGATGGCggg  >  1:139429/1‑90 (MQ=255)
                                                              |                                                                          
TTCGCCATCGGTGTGGCATCCAGCTGGTCGAAAGACAGCGCTGGTGCGGCGGCGCTGGCGGGTGCGGTAGGTTACTTTGTGTTAACCAAAGCGATGGTGACCATCAACCCAGAAATTAACATGGGTGTACTGGCGGG  >  NZ_CP009273/699565‑699701

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