Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,558,3210TG56.4% 35.7 / 46.5 39T6P (ACC→CCC) mdtMmultidrug efflux MFS transporter MdtM
Reads supporting (aligned to +/- strand):  ref base T (2/15);  new base G (22/0);  total (24/15)
Fisher's exact test for biased strand distribution p-value = 5.41e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-08
Rejected as consensus: Frequency below/above cutoff threshold.

GGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC  >  NZ_CP009273/4558236‑4558409
                                                                                     |                                                                                        
ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa                                                                                      >  2:121686/1‑90 (MQ=255)
    ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGgaaaaaaa                                                                                  >  2:401911/1‑89 (MQ=255)
    ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGgaaaaaaa                                                                                  >  1:621581/1‑89 (MQ=255)
    ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa                                                                                  >  2:279793/1‑89 (MQ=255)
     gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaag                                                                                 >  1:201348/1‑88 (MQ=255)
          gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaa                                                                                       >  2:629820/1‑79 (MQ=255)
          gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTaaa                                                                                       <  1:629820/79‑1 (MQ=255)
          gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTggga                                                                            >  1:33427/1‑88 (MQ=255)
          gATCCGTCCACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTgg                                                                              <  2:202596/88‑1 (MQ=255)
          gATCCGTCCACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTgg                                                                              >  1:202596/1‑88 (MQ=255)
              cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATaac                                                                        >  1:162030/1‑90 (MQ=255)
              cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGGATaaa                                                                        >  1:343603/1‑89 (MQ=255)
              cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATaaa                                                                        >  2:501735/1‑89 (MQ=255)
               gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATaaaa                                                                       >  1:56462/1‑88 (MQ=255)
               gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATaaaa                                                                       >  2:340593/1‑88 (MQ=255)
                    cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAacaaacaa                                                                  <  1:340593/90‑1 (MQ=255)
                     aGATACGCAGCAAAGTCATACAAGATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAaacaac                                                                 <  2:276730/90‑1 (MQ=255)
                      gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCGTAaaaaaaaaaa                                                                >  1:493799/1‑88 (MQ=255)
                      gATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCAaaaaaaaaaaaa                                                                >  2:371844/1‑88 (MQ=255)
                              gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAACACCTTTTcc                                                        >  2:270565/1‑90 (MQ=255)
                              gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGGATAAAAAAAAACCCCCTTTcc                                                        >  2:83735/1‑90 (MQ=255)
                                aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAAAAAAAAAAACCCTTTCCgg                                                      >  2:286082/1‑90 (MQ=255)
                                 aaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGGATAAAAAAAAAAACCCTTTCCGga                                                     >  2:7569/1‑90 (MQ=255)
                                     cATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAAc                                                 <  1:121686/90‑1 (MQ=255)
                                           aaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAAc                                           <  2:615985/90‑1 (MQ=255)
                                             aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTgg                                                                              <  2:61055/53‑1 (MQ=255)
                                             aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTgg                                                                              >  1:61055/1‑53 (MQ=255)
                                             aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGGATAAAAAAAAACTCCTTTTCCGGAGAACCCCAACgg                                         >  1:132377/1‑90 (MQ=255)
                                             aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAAAAAAAAAAACCCTTTCCCGGGGACCCCCACgg                                         >  2:436670/1‑90 (MQ=255)
                                                 aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAACTCCTTGTCCGGAGAACCCCAACggggga                                     >  1:172424/1‑88 (MQ=255)
                                                 aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCATAAAAAAAAACTCCTTTTCCGGGGGACCCCCACggggga                                     >  1:506062/1‑88 (MQ=255)
                                                      aTCGGGAAAAACAGCGTGGCTGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  1:691453/90‑1 (MQ=255)
                                                      aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  2:518179/90‑1 (MQ=255)
                                                      aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  2:544421/90‑1 (MQ=255)
                                                      aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  1:286082/90‑1 (MQ=255)
                                                      aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  2:33427/90‑1 (MQ=255)
                                                               aaCAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATa                       <  2:599788/90‑1 (MQ=255)
                                                                                  gggTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAg    <  1:114454/90‑1 (MQ=255)
                                                                                    gTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc  <  1:220411/90‑1 (MQ=255)
                                                                                     |                                                                                        
GGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC  >  NZ_CP009273/4558236‑4558409

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