Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273700,9960CA50.0% 9.5 / 10.8 18P533T (CCG→ACG) nagEPTS N‑acetyl glucosamine transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base C (2/7);  new base A (9/0);  total (11/7)
Fisher's exact test for biased strand distribution p-value = 2.26e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAA  >  NZ_CP009273/700910‑701085
                                                                                      |                                                                                         
cGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGa                                                                                        <  2:170895/90‑1 (MQ=255)
     aCCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGAt                                                                                   <  1:864553/90‑1 (MQ=255)
           gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGGAAACGACAGATaaaaaa                                                                             >  2:657098/1‑88 (MQ=255)
           gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGATaaaaac                                                                             >  1:713577/1‑88 (MQ=255)
               tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtcgtcg                                                                         <  1:135556/90‑1 (MQ=255)
                       cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAACGTCGGATCACCa                                                                 >  2:144768/1‑90 (MQ=255)
                        tGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGATAAAAAAGTCGGATCACCAg                                                                >  2:379037/1‑90 (MQ=255)
                        tGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACCACAAAAAAAAAAGTCGTATTACCAg                                                                >  2:804267/1‑90 (MQ=255)
                               ggTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAAAGTCGTATCACCAGCCGCAgg                                                         >  2:289601/1‑90 (MQ=255)
                                   ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAAAGTCGGATCACCAGCCGCCggggga                                                     >  2:370083/1‑87 (MQ=255)
                                   ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAAATAAAAAAGGCGGATCACCAGCCGCAGGGGCa                                                     >  2:114030/1‑90 (MQ=255)
                                      gACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATc                                                  <  1:604775/90‑1 (MQ=255)
                                                       cAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAacac                                 <  1:800092/90‑1 (MQ=255)
                                                       cAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAacac                                 <  1:327623/90‑1 (MQ=255)
                                                           gCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaac                             >  1:123960/1‑90 (MQ=255)
                                                                          gtgGCGGTGAAAACGACAGATAAAAACGTCGTATCACCCGCCGCAGGGGCAAACGGGGAAAACTTTAAAACCAAACACCCGGTCTGCCTg              >  2:310571/1‑90 (MQ=255)
                                                                            ggCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGa            <  1:370083/90‑1 (MQ=255)
                                                                                  gAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGa      >  2:130085/1‑90 (MQ=255)
                                                                                      ccGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGaaaaa  <  2:123960/90‑1 (MQ=255)
                                                                                      ccGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGaaaaa  <  1:379037/90‑1 (MQ=255)
                                                                                      |                                                                                         
CGATTACCGGTGATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAA  >  NZ_CP009273/700910‑701085

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