Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273701,0090TA50.0% 12.8 / 12.1 22I537N (ATC→AAC) nagEPTS N‑acetyl glucosamine transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (1/10);  new base A (11/0);  total (12/10)
Fisher's exact test for biased strand distribution p-value = 3.40e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGATCG  >  NZ_CP009273/700925‑701098
                                                                                    |                                                                                         
tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGAAAAAAacgtcg                                                                                      >  2:837728/1‑90 (MQ=255)
tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGAAAAAAacgtcg                                                                                      >  1:824638/1‑90 (MQ=255)
tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAaagtcg                                                                                      >  2:25969/1‑90 (MQ=255)
 cGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGAAAAAAAcgacgt                                                                                     >  1:860376/1‑90 (MQ=255)
        cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCa                                                                              <  1:293872/90‑1 (MQ=255)
        cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCa                                                                              <  1:865594/90‑1 (MQ=255)
          ggATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGATAAAAACGTCGTATTACCCGc                                                                            >  1:563238/1‑90 (MQ=255)
             tCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGAAAAAAACGTCGTATCACCAGCCGc                                                                         >  2:570800/1‑90 (MQ=255)
               aGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAACCGACAGAAAAAAAAGTCGGATCACCcgccgcgg                                                                       >  1:85034/1‑88 (MQ=255)
                  ttCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAACGTCGTATCACCAGCCGCAgggg                                                                    >  2:602967/1‑89 (MQ=255)
                     cTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACGACAGAAAAAAAAGGGGGATCACCACCCCCAGGGGCaa                                                                 >  1:1340209/1‑90 (MQ=255)
                          gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACCACAGAAAAAAAAATCGTATCACCAGCCGCCGGGGCAAACggg                                                            >  1:1289136/1‑88 (MQ=255)
                              cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa                                                        <  1:736115/90‑1 (MQ=255)
                              cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa                                                        <  2:1289136/90‑1 (MQ=255)
                                           aGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaa                                           <  2:824638/90‑1 (MQ=255)
                                           aGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaa                                           <  1:570800/90‑1 (MQ=255)
                                               gtgggtgACGGTGTGGCGGTGAAAACGACAGAAAAAAAAGTCGTATCACCACCCCCAGGGGCAAACGGGGAAAAATTCAACacccaccaa                                       >  1:385050/1‑89 (MQ=255)
                                                     gACGGTGTGGCGGGGAAACCCACAGATAAAAACGTCGCATCACCCACCCCAGGGACAAACGGGGAAAACTTTAACACCAACCACCCGTTc                                 >  2:594519/1‑90 (MQ=255)
                                                      aCGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCt                                <  1:1005043/90‑1 (MQ=255)
                                                         gtgtGGCGGGGAAAACGACAGATAAAAAAGGCGTATCACCAGCCGCAGGGGCAAACGGGGAAAAATTTAAAACCAACCCCGCGCTCTGcc                             >  1:1074799/1‑90 (MQ=255)
                                                           gtgGCGGTGGAACCGACAGATAAAAAAGTCGTATCACCCGCCGCAGGGGGAAACGGGAAAAAATTCAAAACCAACCACGCCGTCTGcccg                           >  1:905662/1‑88 (MQ=255)
                                                                     aaCCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGaaa                 >  2:526877/1‑90 (MQ=255)
                                                                      aCCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGaaaa                <  1:1109925/90‑1 (MQ=255)
                                                                                aaaaTCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGgag      <  2:860376/90‑1 (MQ=255)
                                                                                   aTCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtc   <  2:139807/90‑1 (MQ=255)
                                                                                    tcgtcgTATCACCATCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtcg  <  1:837728/90‑1 (MQ=255)
                                                                                    tcgtcgTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtcg  <  1:594519/90‑1 (MQ=255)
                                                                                    |                                                                                         
TCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGATCG  >  NZ_CP009273/700925‑701098

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