Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273701,0090TA45.8% 32.1 / 34.7 48I537N (ATC→AAC) nagEPTS N‑acetyl glucosamine transporter subunit IIABC
Reads supporting (aligned to +/- strand):  ref base T (8/18);  new base A (22/0);  total (30/18)
Fisher's exact test for biased strand distribution p-value = 3.08e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGATCG  >  NZ_CP009273/700921‑701098
                                                                                         |                                                                                         
gATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtc                                                                                          <  1:330080/90‑1 (MQ=255)
 aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCTACAGATAAAAtcg                                                                                         <  2:1893414/90‑1 (MQ=255)
 aTGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtcg                                                                                         <  1:2592921/90‑1 (MQ=255)
  tGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAtcgt                                                                                        <  2:2169937/90‑1 (MQ=255)
    tCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAacgtcg                                                                                      >  2:745779/1‑90 (MQ=255)
      gTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTa                                                                                    <  1:1207096/90‑1 (MQ=255)
       tGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTAt                                                                                   <  2:1867270/90‑1 (MQ=255)
            cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAACGTCGTATCACCa                                                                              >  2:2803741/1‑90 (MQ=255)
            cTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAACGTCGTATCACCa                                                                              >  2:1538755/1‑90 (MQ=255)
             tGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTTGAAACCGCATACAAAAATGTCGGATCACCac                                                                             >  1:2249144/1‑89 (MQ=255)
                  cAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCa                                                                        <  1:401760/90‑1 (MQ=255)
                    ggTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAACGTCGGATCACCCGCCGCCgg                                                                      >  2:474680/1‑90 (MQ=255)
                      ttCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGGAAACGACAGATAAAATTGTCGTATTACCAGCCGCAgggg                                                                    >  2:2801293/1‑89 (MQ=255)
                       tCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGAc                                                                   <  1:622075/90‑1 (MQ=255)
                        ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGGAAACGACAGATAAAAACGTCGTATTACCAACCGCCGGGGCa                                                                  >  1:1004861/1‑90 (MQ=255)
                        ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAg                                                                             <  1:2361799/79‑1 (MQ=255)
                        ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAAACGTCGTATCACCAg                                                                             >  2:2361799/1‑79 (MQ=255)
                        ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAATCGTCGTATCACCAACCCCAGGGGCa                                                                  >  2:1790609/1‑90 (MQ=255)
                        ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGAAAAAAAAGTCGTATCACCAGCCCCCggggga                                                                  >  1:1227900/1‑87 (MQ=255)
                        ccTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACAACAGGTAAAAAAGGCGGATTACCCGCCGCCggggga                                                                  >  1:2127695/1‑87 (MQ=255)
                         cTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGGAAACCACAAAAAAAAAAATCGGATCACCCGCCGCAGGGACaa                                                                 >  1:1955342/1‑90 (MQ=255)
                          tGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAAt                                                                >  2:1008273/1‑90 (MQ=255)
                              gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAATCGTCGGATCACCAGCCGCAGGGGCAATCGTg                                                            >  2:2614456/1‑90 (MQ=255)
                              gAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAACGTCGGATCACCAGCCGCCGGGGCAAACggg                                                            >  2:890162/1‑88 (MQ=255)
                               aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGGGa                                                           >  2:726194/1‑90 (MQ=255)
                               aaGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACCACAGATAAAAACGTGGTTTCACCAACCCCAGGGGGAATCgaga                                                           >  1:1221656/1‑90 (MQ=255)
                                  cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa                                                        <  2:2564332/90‑1 (MQ=255)
                                  cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa                                                        <  1:79516/90‑1 (MQ=255)
                                  cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa                                                        <  2:819936/90‑1 (MQ=255)
                                  cATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa                                                        <  1:2019625/90‑1 (MQ=255)
                                  cATTCGCCAGCAAAGCGGTCGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAAAGTCGTATCACCAACCGCAGGGGGAAACGTGgaaa                                                        >  1:590388/1‑90 (MQ=255)
                                  acagcGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGaaaa                                                        <  2:2320275/86‑1 (MQ=255)
                                    ttCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAAAGTCGTATCACCAGCCGCAGGGGCAAACGTGaaaaac                                                      >  1:1016937/1‑88 (MQ=255)
                                     tCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAAACCACAGAAAAAAAAGTCGTATCACCAGCCGCAGGGGCAAAAGTGAAAAACt                                                     >  1:462663/1‑90 (MQ=255)
                                      cGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGGGAAAACGACAGATAAAAACGTCGGATCACCCGCCGCAGGGGGAAACGTGGAAAACtt                                                    >  1:525822/1‑90 (MQ=255)
                                       gCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTc                                                   <  1:726194/90‑1 (MQ=255)
                                        ccAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCa                                                  <  2:2512950/90‑1 (MQ=255)
                                        ccAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCa                                                  <  1:2562412/90‑1 (MQ=255)
                                               aGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaa                                           <  2:2237764/90‑1 (MQ=255)
                                               aGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACaccaa                                           <  2:2606353/90‑1 (MQ=255)
                                                gCGGTGGGTGACGGTGTGGCGGGGAAAACGGCAAAAAAAAAAGTCGTATTACCAGCCGCCGGGGCAAACGTGGAAAAATTCAAAacccgg                                          >  1:2320275/1‑87 (MQ=255)
                                                   gtgggtgACGGTGTGGCGGTGAAACCGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGGGAAAAACTTCAACACCAACCAc                                       >  2:2612965/1‑90 (MQ=255)
                                                        tgACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGtt                                  <  2:65601/90‑1 (MQ=255)
                                                          aCGGTGTGGCGGTGAAAACGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAAACTTCAACACCAACCACGCGTTCt                                >  2:481698/1‑90 (MQ=255)
                                                             gtgtGGCGGTGAAACCGACAGATAAAAACGTCGTATCACCAGCCGCAGGGGCAATCGTGAAAATCTTCAACACCAACCACGCGGTCTGcc                             >  2:2312593/1‑90 (MQ=255)
                                                               gtgGCGGTGAAAACCACAAATAAAAAAGTCGTATCACCAGCCGCAGGGGGAAACGTGAAAAAATTTAAAAACCACCACCCGCTCTGCCTg                           >  1:1569865/1‑90 (MQ=255)
                                                                 ggCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGa                         <  2:1291368/90‑1 (MQ=255)
                                                                 ggCGGTGAAAACGACAGATAAAAAAGTCGTATCACCAGCCCCAGGGGCAAACGGGAAAAAATTCAAAACCAACCCCCCGGTCTGCCTggg                         >  1:2216166/1‑89 (MQ=255)
                                                                        gAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGGCAAACGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGa                   >  2:1891294/1‑90 (MQ=255)
                                                                        gAAACCGACAGATAAAATCGGCGTGTCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGa                   >  2:644936/1‑90 (MQ=255)
                                                                        gAAAACGACAGAAAAAAACGTCGTATCACCAGCCGCAGGGGCAAACGTGAAAAAATTTAACACCCACCACCCGTTCTGCCTGGGAACCGa                   >  1:1688698/1‑90 (MQ=255)
                                                                             cGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAg              >  1:95754/1‑90 (MQ=255)
                                                                                        aTCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtc   <  1:926190/90‑1 (MQ=255)
                                                                                         tcgtcgTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtcg  <  1:2624220/90‑1 (MQ=255)
                                                                                         tcgtcgTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGAtcg  <  1:1166662/90‑1 (MQ=255)
                                                                                         |                                                                                         
GATGTCGTGGCACTGGATCAGGTTCCTGACGAAGCATTCGCCAGCAAAGCGGTGGGTGACGGTGTGGCGGTGAAACCGACAGATAAAATCGTCGTATCACCAGCCGCAGGGACAATCGTGAAAATCTTCAACACCAACCACGCGTTCTGCCTGGAAACCGAAAAAGGCGCGGAGATCG  >  NZ_CP009273/700921‑701098

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