Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,290,9790TG58.3% 21.1 / 27.3 24T639T (ACA→ACCmdtOmultidrug efflux transporter permease subunit MdtO
Reads supporting (aligned to +/- strand):  ref base T (2/8);  new base G (14/0);  total (16/8)
Fisher's exact test for biased strand distribution p-value = 6.12e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGA  >  NZ_CP009273/4290890‑4291063
                                                                                         |                                                                                    
ggggCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGgt                                                                                      >  2:850493/1‑90 (MQ=255)
      gTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTctct                                                                                >  2:168886/1‑90 (MQ=255)
      gTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGctct                                                                                >  2:420361/1‑90 (MQ=255)
       tCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGCTCtt                                                                               >  2:525745/1‑90 (MQ=255)
           gTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGCTCtttcag                                                                           >  2:643277/1‑86 (MQ=255)
             cTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCg                                                                         <  1:90760/90‑1 (MQ=255)
             cTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCg                                                                         <  2:227363/90‑1 (MQ=255)
                gCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGCTCTTTAAGCggtt                                                                      >  2:112871/1‑87 (MQ=255)
                gCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGCTCTTCAAGCGGTa                                                                      >  2:419951/1‑90 (MQ=255)
                      gggCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTgag                                                                <  2:352808/90‑1 (MQ=255)
                       ggCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTgagg                                                               <  2:57019/90‑1 (MQ=255)
                               gCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGc                                                       <  2:853705/90‑1 (MQ=255)
                                                  cTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGAcc                                    <  1:419951/90‑1 (MQ=255)
                                                       ttAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAg                               <  1:524680/90‑1 (MQ=255)
                                                         aataaGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGGTATTTTAGGAGAACGGGTGAGTGCGGGTGCCCGACCGGCAGCg                             >  2:129860/1‑90 (MQ=255)
                                                              ggTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGTTTTTTTGGGGGAACGGCTGAGGGGGGGTGCCAGACCGGCAGCGGAttt                        >  2:524680/1‑90 (MQ=255)
                                                              ggTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGTTTTTTTAGGAGAACGGCTGGGGGGGGGTGCCAGAACGGCAGCGGAttt                        >  2:760003/1‑90 (MQ=255)
                                                              ggTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGTTATTTTAGGGGAACGGGTGGGGGCGGGTGGCCGAACGGGAGCGGAttt                        >  2:446250/1‑90 (MQ=255)
                                                              ggTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGGTATTTTAGGAGAAAGGCTGAGGGGGGGTGCCAGACCGGCAGCGGAttt                        >  2:76719/1‑90 (MQ=255)
                                                              ggTGGGCAGGATGGCCTGAGAGGCGGGGGGCTCTTCAAGCGTTATTTGAGGAGAAAGGCTGAGTGCGGGTGCCAGACCGGCAGCGGAttt                        >  2:625188/1‑90 (MQ=255)
                                                                      ggATGGCCTGAGAGGCGGGGGGCTCTTCAAGCGGTATTTTAGGAGAACGGCTGAGTGCGGGTGCCAGACCGGCAGCGGATTTTTCCAggg                >  2:420557/1‑90 (MQ=255)
                                                                       gATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGc               <  1:521513/90‑1 (MQ=255)
                                                                        aTGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTTTTTTAGGGGAACGGCTGAGGGCGGGTGCCAGACCGGCAGCGGATTTCTCCCGGGCg              >  2:514435/1‑90 (MQ=255)
                                                                        aTGGCCTGAGAGGCGGGGGGCTCTTTAAACGGTATTTTAGGGGAAAGGCTGAGGGCGGGTGCCAGACCGGCAGCGTATTTTTCCAggggg              >  2:569432/1‑88 (MQ=255)
                                                                                    gCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTTTCCAGGGCGTCGGCAAACTGa  >  2:206441/1‑90 (MQ=255)
                                                                                         |                                                                                    
GGGGCTGTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGA  >  NZ_CP009273/4290890‑4291063

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