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

TAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATC  >  NZ_CP009273/4292402‑4292563
                                                                         |                                                                                        
tAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCg                                                                          <  1:462780/90‑1 (MQ=255)
  gagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGAc                                                                        >  2:14338/1‑90 (MQ=255)
    ggCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGGCCGTAAATAACGACaa                                                                      >  2:462780/1‑90 (MQ=255)
       caacaaCGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCGAAAAAACGACAAAgg                                                                   >  1:158286/1‑90 (MQ=255)
                    ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCGCAAAAAACGGCAAAGGCGGCGGCGgaaaa                                                      >  1:342960/1‑90 (MQ=255)
                      acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGGAAAAAAGGACAATGGCGACGGCGaaaaaaa                                                    >  2:387687/1‑88 (MQ=255)
                         gCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCCTAAAAAACGACAAAGGCGGCGGCGGAAAAGACCa                                                 >  1:282895/1‑90 (MQ=255)
                            ttAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCTAAAAAACGACAAAGGCGGCCGCGGAAAAAAACAGcc                                              >  1:268423/1‑90 (MQ=255)
                              aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGcccc                                            <  1:209905/90‑1 (MQ=255)
                              aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGTAAAAAACGACAAAGGCGGCCGCGGAAAAAACCAGcccc                                            >  2:44865/1‑90 (MQ=255)
                                   cACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAATAACGGCAAAGGCGGCGGGGGAAAAAAACCGCCCCAAGCg                                       >  1:429515/1‑90 (MQ=255)
                                     cGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGAt                                     <  1:288585/90‑1 (MQ=255)
                                                ggATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCAAGCGATGGGGGCGCATc                          >  2:245360/1‑90 (MQ=255)
                                                   tAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCaaa                       <  2:420106/90‑1 (MQ=255)
                                                     gTCGAGCATGGCGGGGGAGGGTTTCCCCTAAAAAAGGACAAAGGCGACGGGGAAAAAAACCAGCCCCCAACGGAGGGGGGGCGTCaaaaa                     >  1:250223/1‑90 (MQ=255)
                                                     gTCGAGCATGGCGGGGAAGGTTTGCCCCGAAAAAGCGGCAATTGCGACGGCGAAAAAAACCCGCCCCAAGGGATGGGGGCGCCTCaaaaa                     >  2:291419/1‑90 (MQ=255)
                                                     gTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCAAGCGATGGGGGCGCATCaaaaa                     >  2:226885/1‑90 (MQ=255)
                                                             tGGCGGGGGAGGGTTGCCCGGAAAAAGCGACAAAGGCGACGGCGAAAAAAAACAGCCCCCAGCGATGGGGGCGCATCAAAAAAAagcagc             >  1:25069/1‑90 (MQ=255)
                                                             tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc             <  1:14338/90‑1 (MQ=255)
                                                             tGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGc             <  2:214262/90‑1 (MQ=255)
                                                                   ggAAGGGTTGCCCGTAAAAAGCGACAAAGGCGACGGCGAAAAAGACCAGCCCCCAGCGATGGGGGGGCATCAAAAAAAAGCAGCCCATCa       >  1:48947/1‑90 (MQ=255)
                                                                        gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc  <  1:492843/90‑1 (MQ=255)
                                                                        gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATc  <  1:291419/90‑1 (MQ=255)
                                                                         |                                                                                        
TAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATC  >  NZ_CP009273/4292402‑4292563

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: