Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,292,4750TG55.0% 15.1 / 17.7 20T141P (ACC→CCC) mdtOmultidrug efflux transporter permease subunit MdtO
Reads supporting (aligned to +/- strand):  ref base T (3/6);  new base G (11/0);  total (14/6)
Fisher's exact test for biased strand distribution p-value = 2.17e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.32e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC  >  NZ_CP009273/4292391‑4292552
                                                                                     |                                                                             
gCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTGc                                                                           >  2:309005/1‑90 (MQ=255)
 cAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGccc                                                                         <  2:53468/90‑1 (MQ=255)
  aaGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCg                                                                        >  2:110616/1‑90 (MQ=255)
             gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGGCCGGAAAAAACGAc                                                             >  2:440528/1‑90 (MQ=255)
             gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGTAAAAAGCGAc                                                             >  1:186241/1‑90 (MQ=255)
               ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCTAAATAGCGACaa                                                           >  1:350551/1‑90 (MQ=255)
                   aacaacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGGTGGCCGTAAAAAACGACAAAGGc                                                       >  2:445372/1‑90 (MQ=255)
                                   aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc                                       <  2:21022/90‑1 (MQ=255)
                                         aaGCGCACGACCGCTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGcccc                                 <  2:314020/90‑1 (MQ=255)
                                           gcgcACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCaa                               <  1:44496/90‑1 (MQ=255)
                                              cACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCTAAATAACGACAAAGGCGGCGGCGAAAAAAAACAGCCCCCAGCg                            >  1:418142/1‑90 (MQ=255)
                                                cGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGTAAAAAAGACAATGGGGACGGCGAAAAAAAACCGCCCCCAggggt                          >  1:335053/1‑86 (MQ=255)
                                                     aCTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCGTAAAAAACGACAATGGCGGCGGCGAAAAAAAACAGCCCCAAGCGATggggg                     >  1:428412/1‑88 (MQ=255)
                                                                gTCGAGCATGGCGGGGGGGGGTTTCCCGGAAAAAACGACAAAGGGGGCCGCGGAAAAAACCAACCCCCAACGAAGGGGGCGGCTCaaaaa          >  2:413503/1‑90 (MQ=255)
                                                                gTCGAGCATGGCGGGGGAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGGAAAAAAACAGCCCCCAGCGATGGGGGGGCATCaaaaa          >  2:242458/1‑90 (MQ=255)
                                                                gTCGAGCATGGCGGGGGAGGGTTTCCCGGAAAAAACGACAAAGGGGGCGGCGGAAAAAAACAGCCCCCAACGATGGGGGGGCCTCaaaaa          >  2:398410/1‑90 (MQ=255)
                                                                gTCGAGCATGGCGGGGGAGGGTTGCCCGTAAAAAGCGACAAAGGCGGCGGCGAAAAAAACCAGCCCCCAACGATGGGGGGGCCTCaaaaa          >  1:255727/1‑90 (MQ=255)
                                                                  cGAGCATGGCGGGGGAGGTGTTCCCCTAAAAAGCGACAAAGGCGGTGGCGGAAAAGAACCGCCCCCAGCGGTGGGGGCGCCTGaaaaaaa        >  1:448786/1‑88 (MQ=255)
                                                                     gCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGc     <  2:74649/90‑1 (MQ=255)
                                                                        tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc  <  1:419151/90‑1 (MQ=255)
                                                                                     |                                                                             
GCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC  >  NZ_CP009273/4292391‑4292552

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: