Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,964,4660TG53.8% 11.0 / 16.2 13V325G (GTA→GGA) wecBUDP‑N‑acetylglucosamine 2‑epimerase (non‑hydrolyzing)
Reads supporting (aligned to +/- strand):  ref base T (2/4);  new base G (7/0);  total (9/4)
Fisher's exact test for biased strand distribution p-value = 2.10e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACC  >  NZ_CP009273/3964378‑3964555
                                                                                        |                                                                                         
gaagCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTa                                                                                          <  2:299684/90‑1 (MQ=255)
                        gtgCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGGGGCTGGGAGGGACGGGGAAGCAGCGGAATGTc                                                                  >  2:119511/1‑90 (MQ=255)
                                       gATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGGGGGTCTGGGAGGCACGGATAAGCAGCGAATTTTCGAGGGAGAGAcgcgg                                                   >  1:210585/1‑89 (MQ=255)
                                             aCTGAGCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGGCTGGGGGGGACGGAAAAGCAGGGGAATGTGGGGGGAGGGAAGCGGCttttt                                             >  1:519649/1‑89 (MQ=255)
                                                gAGCGTCCGGAAGCGGTGACGGCGGGGTGGGGGTGTCTGGGGGGCCAGGGTTACCAGGGAAATGTGGAGGGAGGGGGGCGGCTTTTtaaa                                          >  2:511142/1‑90 (MQ=255)
                                                  gCGTCCGGAAGCGGTGACGGCGGGGACGGGGCGTCTGGGAGGGACGGGTAAAAAGCGAAATTTCGAGGGAGGGAAGCGTCTTTTtaaaaa                                        >  2:114463/1‑88 (MQ=255)
                                                      ccGGAAGCGGTGACGGCGGGGACGGGGCGGCTGGGAGGGACGGGTAAGCAGCGGATTGTCGGGGGAGGGGCGCGGGTTTTTAAAAAAGaa                                    >  2:9118/1‑90 (MQ=255)
                                                            gCGGTGACGGCGGGGACGGGGCGGCTGGGAGGGACGGGTAAACAGCGAAATTTCGGGGGAGGGGAGCGGCTTTTTAAAaaagaaaaagaa                              >  1:512068/1‑90 (MQ=255)
                                                                              gTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCcgc            <  1:279302/90‑1 (MQ=255)
                                                                              gTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCcgc            <  1:391815/90‑1 (MQ=255)
                                                                                 cGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGcc         <  1:98139/90‑1 (MQ=255)
                                                                                  gTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGccc        >  1:372491/1‑90 (MQ=255)
                                                                                  gTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGccc        >  2:447738/1‑90 (MQ=255)
                                                                                        tAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAAcc  >  2:178602/1‑90 (MQ=255)
                                                                                        |                                                                                         
GAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGTCTGGTAGGCACGGATAAGCAGCGAATTGTCGAGGAAGTGACGCGTCTTTTAAAAGACGAAAACGAATATCAAGCTATGAGCCGCGCCCATAACC  >  NZ_CP009273/3964378‑3964555

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