Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,623,1510AG67.6% 35.6 / 33.9 37I271T (ATC→ACC) yhiIHlyD family efflux transporter periplasmic adaptor subunit
Reads supporting (aligned to +/- strand):  ref base A (4/8);  new base G (14/11);  total (18/19)
Fisher's exact test for biased strand distribution p-value = 2.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.72e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TTGACGCGGAACATCAGTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCGGCAA  >  NZ_CP009273/3623028‑3623286
                                                                                                                           |                                                                                                                                       
ttGACGCGGAACATCAGTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCAc                                                                                                                          >  2:339691/1‑139 (MQ=255)
    cgcgGAACATCAGTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCGTCGAGGATCAGGCGGGCTTCACCGcc                                                                                                                      <  1:152869/139‑1 (MQ=255)
    cgcgGAACATCAGTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGcc                                                                                                                      <  2:540036/139‑1 (MQ=255)
           cATCAGTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGttt                                                                                                                >  1:586271/1‑138 (MQ=255)
                 tttCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTg                                                                                                         >  2:527774/1‑139 (MQ=255)
                   tCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGcc                                                                                                       >  1:434764/1‑139 (MQ=255)
                     aGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGccc                                                                                                      <  1:416101/138‑1 (MQ=255)
                          ttCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGt                                                                                                >  2:321653/1‑139 (MQ=255)
                          ttCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGt                                                                                                >  1:690016/1‑139 (MQ=255)
                           tCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGtt                                                                                               >  2:201109/1‑139 (MQ=255)
                            cATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTc                                                                                              >  1:593547/1‑139 (MQ=255)
                                    gTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTgg                                                                                        <  2:136431/137‑1 (MQ=255)
                                         gACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGgaa                                                                                  <  2:451503/138‑1 (MQ=255)
                                          aCGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGgaaga                                                                                <  2:85644/139‑1 (MQ=255)
                                          aCGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGgaaga                                                                                <  1:201109/139‑1 (MQ=255)
                                          aCGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGgaaga                                                                                <  1:472466/139‑1 (MQ=255)
                                             gTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAg                                                                             <  1:246088/139‑1 (MQ=255)
                                                        gAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACg                                                                  >  1:674336/1‑139 (MQ=255)
                                                        gAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACg                                                                  >  2:567958/1‑139 (MQ=255)
                                                              ggCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCt                                                             >  1:91977/1‑138 (MQ=255)
                                                              ggCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCt                                                             >  1:567164/1‑138 (MQ=255)
                                                                  acacTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGAt                                                        <  1:567958/139‑1 (MQ=255)
                                                                   cacTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCt                                                             >  2:451318/1‑133 (MQ=255)
                                                                       ggCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGAcc                                                   <  1:74789/139‑1 (MQ=255)
                                                                        gCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCa                                                  >  2:663958/1‑139 (MQ=255)
                                                                        gCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCa                                                  >  2:689625/1‑139 (MQ=255)
                                                                        gCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCa                                                  >  2:359162/1‑139 (MQ=255)
                                                                                aCTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGc                                          <  2:586271/139‑1 (MQ=255)
                                                                                           cAGGAATACGCAGATCTGGCGCGGCATCGAGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGc                                <  1:339691/138‑1 (MQ=255)
                                                                                            aGGAATACCCAGATCTGGCGCGGCACCGGGGGTCGGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCt                              <  2:567164/139‑1 (MQ=255)
                                                                                                       gATCTGGCGCGGCATCGAGGGTCGGCCGGGCTTCACCCCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCAc                   <  2:399594/139‑1 (MQ=255)
                                                                                                           gggCGCGGCATCGGGGATCAGCCGGGTTTCACCGCCCAGTTTCAGCGTGCCCGCCTTTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCg               <  2:546185/138‑1 (MQ=255)
                                                                                                                     tCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCgg     <  1:656280/139‑1 (MQ=255)
                                                                                                                     tCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCgg     <  1:359162/139‑1 (MQ=255)
                                                                                                                        aGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGAt                                                        <  2:130853/85‑1 (MQ=255)
                                                                                                                        aGGGTCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGAt                                                        >  1:130853/1‑85 (MQ=255)
                                                                                                                        aGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCGGCaa  >  1:617475/1‑139 (MQ=255)
                                                                                                                           |                                                                                                                                       
TTGACGCGGAACATCAGTTTCAGCCGTTCATCGCTGGTTTCGACGGTTTTTGGCGTGAACTGGGCGACACTGGCGACAAAACTGATGGTTGCAGGAATACGCAGATCTGGCGCGGCATCGAGGATCAGCCGGGCTTCACCGCCCAGTTTCAGCGTGCCCGCCTGTTCGGTTGGCAGGAAGAAAGTCATATAGACGTCGCTGAGATCGACCATATTCAGCACCCGACCGCCTGCCGCCAGCACTTCGCCTGGCTCGGCAA  >  NZ_CP009273/3623028‑3623286

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑