Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,292,4710AG61.9% 6.3 / 18.2 21F142S (TTC→TCC) mdtOmultidrug efflux transporter permease subunit MdtO
Reads supporting (aligned to +/- strand):  ref base A (2/6);  new base G (13/0);  total (15/6)
Fisher's exact test for biased strand distribution p-value = 5.16e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.81e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC  >  NZ_CP009273/4292387‑4292552
                                                                                     |                                                                                 
atcaGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAgggt                                                                              >  2:229809/1‑88 (MQ=255)
    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)
                                                                                     |                                                                                 
ATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGC  >  NZ_CP009273/4292387‑4292552

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