Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,290,9790TG52.9% 13.7 / 20.1 17T639T (ACA→ACCmdtOmultidrug efflux transporter permease subunit MdtO
Reads supporting (aligned to +/- strand):  ref base T (1/7);  new base G (9/0);  total (10/7)
Fisher's exact test for biased strand distribution p-value = 4.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.53e-03
Rejected as consensus: Frequency below/above cutoff threshold.

GTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGA  >  NZ_CP009273/4290896‑4291063
                                                                                   |                                                                                    
gTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTctct                                                                                >  1:822976/1‑90 (MQ=255)
gTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGctct                                                                                >  2:326732/1‑90 (MQ=255)
  ccAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTc                                                                              <  2:397507/90‑1 (MQ=255)
      tCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGGCTCTTCAAGc                                                                          >  2:308268/1‑90 (MQ=255)
          gCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTTAAGCGttt                                                                      >  1:583359/1‑89 (MQ=255)
            aaGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTAtt                                                                    <  2:822976/90‑1 (MQ=255)
                gggCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCCAGCGTTTTTTtag                                                                >  2:211248/1‑90 (MQ=255)
                                                 ttAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAg                               <  1:919249/90‑1 (MQ=255)
                                                      aaGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTAt                          <  2:768960/90‑1 (MQ=255)
                                                      aaGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTAt                          <  2:527481/90‑1 (MQ=255)
                                                       aGGTGGGCAGGATGGCCTGAGAGGCGGGGGTCTCTTCAAGCGTTATTTTAGGAGAACGGCTGAGGGGGGTTGCCAGACCGGCAGCGTAtt                         >  2:106786/1‑90 (MQ=255)
                                                                  aTGGCCTGAGAGGCGGGGGTCTCTTCAAGGGTTATTTGTGGGGAACGGGTGGGGGGGGGTGGCCGACCCGGAGCGTGTTTTTTCAGGGCg              >  1:325333/1‑90 (MQ=255)
                                                                   tGGCCTGAGAGGCGGGGGTCTCTTCAAGCGTTTTTTTAGGAGAACGGGTGAGTGCGGGTGCCAGACCGGCAGCGGATTTTTCCAgggggg             >  2:2080/1‑87 (MQ=255)
                                                                         gagaGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCaa       <  1:658444/90‑1 (MQ=255)
                                                                           gagGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAAc     <  1:308268/90‑1 (MQ=255)
                                                                              gCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGGATTTCTCCAGGGCGTCGGCAAACTGa  >  2:58320/1‑90 (MQ=255)
                                                                              gCGGGGGTCTCTTCAAGCGGTTTTTGAGGGGAACGGCTGAGGGCGGGTGGCCGACCGGCCGGGGATTTTTCCCGGGCGTCGGCAAAATgc  >  1:560906/1‑89 (MQ=255)
                                                                                   |                                                                                    
GTCCAGTCTGGCAAGCGGGCGAAAAGCTGGCAGACGTGTTGCTCCTGTTTTAATAAGGTGGGCAGGATGGCCTGAGAGGCGGGTGTCTCTTCAAGCGTTATTTGAGGAGAACGGCTGAGTGCGGTTGCCAGACCGGCAGCGTATTTCTCCAGGGCGTCGGCAAACTGA  >  NZ_CP009273/4290896‑4291063

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