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

TCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCAT  >  NZ_CP009273/4292388‑4292556
                                                                                       |                                                                                 
tcaGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGttt                                                                                 <  2:58407/90‑1 (MQ=255)
  aGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGc                                                                               <  2:26454/90‑1 (MQ=255)
           ggATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGGAAATa                                                                      >  2:321862/1‑90 (MQ=255)
                gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTTCCGTAAAAAGCggc                                                                 >  2:310502/1‑88 (MQ=255)
                gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTGCCGGAAAAAGCggc                                                                 >  1:177238/1‑88 (MQ=255)
                                  ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCCTAAAAAGCGACAAAGGGGGCGGCGaaaaa                                               >  2:231593/1‑90 (MQ=255)
                                    acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAAAAACGACAAAGGCGGCGGCGgaaaaaa                                             >  2:379110/1‑88 (MQ=255)
                                    acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAAAAACGACAATGGGGACGGCGAAAAAGa                                             >  2:55666/1‑90 (MQ=255)
                                      aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGTTTGCCCGGAAATAACGACAATGGCGGCGGCGGAAAAAAcc                                           >  1:53643/1‑90 (MQ=255)
                                      aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGGCCCGAAAAAGCGACAAAGGGGACGGGGAAAAAGAcc                                           >  2:244100/1‑90 (MQ=255)
                                      aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc                                           <  1:69856/90‑1 (MQ=255)
                                      aGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc                                           <  2:278550/90‑1 (MQ=255)
                                       gCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAAAAACGACAATGGCGGCGGCGAAAAAAACCa                                          >  2:75516/1‑90 (MQ=255)
                                          ttAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGCAAAAAACGACAAAGGGGACGGCGGAAAAAACCAGcc                                       >  1:282040/1‑90 (MQ=255)
                                                          ttCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGGAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTgcg                       <  1:168017/90‑1 (MQ=255)
                                                                   gTCGAGCATGGCGGGGGAGGGTTGCCCGTAAAAAACGACAATGGCGACGGCGAAAAAAAACAGCCCCCAACGATGGGGGCGCCTCaaaaa              >  2:349569/1‑90 (MQ=255)
                                                                   gTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGGGGCGGCGAAAAAAACCACCCCCAAACGGTGGGGGCGCAGCaaaaa              >  2:192254/1‑90 (MQ=255)
                                                                           tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc      <  1:349569/90‑1 (MQ=255)
                                                                           tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc      <  1:310502/90‑1 (MQ=255)
                                                                           tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc      <  1:244100/90‑1 (MQ=255)
                                                                               ggggAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGCCCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCAt  <  2:205009/90‑1 (MQ=255)
                                                                                       |                                                                                 
TCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCAT  >  NZ_CP009273/4292388‑4292556

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: