Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273700,0310TG52.9% 11.1 / 14.6 17V211G (GTT→GGT) nagEPTS N‑acetyl glucosamine transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (5/3);  new base G (9/0);  total (14/3)
Fisher's exact test for biased strand distribution p-value = 8.24e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.71e-02
Rejected as consensus: Frequency below/above cutoff threshold.

GTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCCGATCATGATGTTCGGTC  >  NZ_CP009273/699944‑700120
                                                                                         |                                                                                         
gTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGgtt                                                                                         >  2:210302/1‑90 (MQ=255)
     ctgATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCAc                                                                                    <  1:494563/90‑1 (MQ=255)
       gATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTTCACgg                                                                                  >  2:494563/1‑90 (MQ=255)
                    cTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGtt                                                                                          <  2:864185/69‑1 (MQ=255)
                           aGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGTTTTTCACGGTGACATTAACCGCTTCTATg                                                              >  1:358857/1‑90 (MQ=255)
                           aGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTCCACGGGGACATTTACCGCTTCTATg                                                              >  1:835919/1‑90 (MQ=255)
                           aGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTCCACGGGGACATTAACCGCTTCTATg                                                              >  2:311155/1‑90 (MQ=255)
                                cTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCGCCAACGCGGCGGGGACCGGTTTTCCCGGGGAAATTTACCGCTTCTATTCCggg                                                         >  1:876854/1‑89 (MQ=255)
                                cTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTTCCCCGGGACATTAACCCCTTCTATGCCggg                                                         >  1:666597/1‑89 (MQ=255)
                                cTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTCCACCGGGACATTAACCGCTTCTATGCCggg                                                         >  2:691078/1‑89 (MQ=255)
                                    acacCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTCCACCGGGACATTAACCGCCTCTATGCCGGGGACg                                                     >  1:887726/1‑90 (MQ=255)
                                          tCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTCCACGGTGACATTAACCGCTTCTATGCCGGGGACGGCACcg                                               >  2:854368/1‑90 (MQ=255)
                                            gCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGAGGGTTTTCCCCGGGGACATTAACCGCCTCTTTTTCGGTGGCGGCTCcacg                                             >  2:925546/1‑90 (MQ=21)
                                                ggTTCCAGATTGGTGAATTCACCCACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGa                                         <  1:250569/90‑1 (MQ=255)
                                                 gTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTCCCCGGGGACATTAACCGCTTCTATGCCGGGGACGGCACCCCgggggt                                        >  2:581216/1‑88 (MQ=255)
                                                           ggTGAATTCACCAACGCGGCGGGTACGGGTTTCCACGGGGACATTAACCGCTTCTATGCCGGGGACGGCACCCCGGGGGTGTTCATGTcc                              >  2:792101/1‑90 (MQ=255)
                                                                aTTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATTCCGGAGACGGCACCGCGGGGAAGTTCCTGTTCAGctt                         >  2:818481/1‑90 (MQ=255)
                                                                  tCACCAACGCGGCGGGGACGGGTTTCCACCGGGACATTAAACGCTTCTATGCCGGGGGCGGGACCGCGGGGGTGTTTATGTCCGGcttct                       >  1:411434/1‑90 (MQ=255)
                                                                      cAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTccc                   <  1:581216/90‑1 (MQ=255)
                                                                                         ttttCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCCGATCATGATGTTCGGTc  <  2:193524/90‑1 (MQ=255)
                                                                                         |                                                                                         
GTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCCGATCATGATGTTCGGTC  >  NZ_CP009273/699944‑700120

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: