Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273700,0260TG54.5% 16.1 / 24.3 22G209G (GGT→GGGnagEPTS N‑acetyl glucosamine transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (5/5);  new base G (12/0);  total (17/5)
Fisher's exact test for biased strand distribution p-value = 9.57e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.77e-03
Rejected as consensus: Frequency below/above cutoff threshold.

ATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCC  >  NZ_CP009273/699937‑700103
                                                                                         |                                                                               
atcaACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGt                                                                                 <  2:536863/90‑1 (MQ=255)
       gTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGgtt                                                                        >  2:210302/1‑90 (MQ=255)
            ctgATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCAc                                                                   <  1:494563/90‑1 (MQ=255)
              gATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGGACGGGTTTTCACgg                                                                 >  2:494563/1‑90 (MQ=255)
                           cTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACgggt                                                                         >  1:864185/1‑67 (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)
                                   ggTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTAc                                                                              >  1:549618/1‑57 (MQ=255)
                                   ggTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTAc                                                                              <  2:549618/57‑1 (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)
                                                                                         |                                                                               
ATCAACCGTCTGCTGATCCCAACCGGTCTGCATCAGGTACTGAACACCATCGCCTGGTTCCAGATTGGTGAATTCACCAACGCGGCGGGTACGGTTTTCCACGGTGACATTAACCGCTTCTATGCCGGTGACGGCACCGCGGGGATGTTCATGTCCGGCTTCTTCCC  >  NZ_CP009273/699937‑700103

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: