Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273700,9960CA63.6% ‑6.4 / 16.6 22P533T (CCG→ACG) nagEPTS N‑acetyl glucosamine transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base C (1/7);  new base A (14/0);  total (15/7)
Fisher's exact test for biased strand distribution p-value = 4.69e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.87e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTG  >  NZ_CP009273/700921‑701073
                                                                           |                                                                             
gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtc                                                                 <  2:1148224/90‑1 (MQ=255)
  tGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAacgt                                                               >  2:351275/1‑90 (MQ=255)
           aCTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCAcc                                                      <  2:1004205/90‑1 (MQ=255)
            cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAACGTCGTATCACCa                                                     >  2:1207263/1‑90 (MQ=255)
                aTCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCg                                                 <  2:981349/90‑1 (MQ=255)
                   aGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAg                                              >  2:448188/1‑90 (MQ=255)
                        ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGGAAACCACAAAAAAAAAACTCGGATCACCACCCGCAGGGGCa                                         >  1:699129/1‑90 (MQ=255)
                        ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAACGTCGTATCACCAGCCGCCGGGACa                                         >  2:1209030/1‑90 (MQ=255)
                        ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACCACACATAAAAACGCCGTATTACCCACCGCCGGGGCa                                         >  2:971547/1‑90 (MQ=255)
                              gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGATAAAAACGTCGGATTACCAACCCCCGGGGCAAACggg                                   >  1:200180/1‑88 (MQ=255)
                              gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAACGTCGTATCACCAGCCCCCGGGGCAATCGTg                                   >  1:165193/1‑90 (MQ=255)
                              gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGAAAGATAAAAAAGTCGCATCCCCAGCCGCAGGGACAAACGTg                                   >  1:628919/1‑90 (MQ=255)
                              gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACCACAGATAAAAAAGTCGTATCACCCGCCGCAGGGGCAAACGTg                                   >  1:1023186/1‑90 (MQ=255)
                               aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAATCGTGa                                  >  2:429647/1‑90 (MQ=255)
                                  cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa                               <  2:367367/90‑1 (MQ=255)
                                  cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa                               <  2:443957/90‑1 (MQ=255)
                                  cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAATCGTGgaaa                               >  2:246008/1‑90 (MQ=255)
                                           gCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACCACAGATAAAAACGTCGTATTACCCGCCGCCGGGGCAAACGTGGAAAACTTTaaaa                      >  1:352383/1‑88 (MQ=255)
                                                        tggcGGTGTGGCGGTGAAACCGACAGATAAAACCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGtt         <  1:407168/87‑1 (MQ=255)
                                                          aCGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCt       <  2:902565/90‑1 (MQ=255)
                                                             gtgtGGCGGTGAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGTGAAAAACTTTAACACCAACCACCCGGTCTGcc    >  1:499475/1‑90 (MQ=255)
                                                               gtgGCGGTGAAAACGACAGAAAAAAATGTCGTATCACCAGCCGCAGGGACAAACGTGAAAAACTTTAACACCAACCACCCGGTCTGCCTg  >  2:1233024/1‑90 (MQ=255)
                                                                           |                                                                             
GATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTG  >  NZ_CP009273/700921‑701073

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