Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,558,3210TG52.6% 21.8 / 24.2 19T6P (ACC→CCC) mdtMmultidrug efflux MFS transporter MdtM
Reads supporting (aligned to +/- strand):  ref base T (1/8);  new base G (10/0);  total (11/8)
Fisher's exact test for biased strand distribution p-value = 1.19e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.68e-05
Rejected as consensus: Frequency below/above cutoff threshold.

GGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC  >  NZ_CP009273/4558240‑4558409
                                                                                 |                                                                                         
ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGgaaaaaaa                                                                                   >  1:326213/1‑89 (MQ=255)
 gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaaa                                                                                  >  1:397652/1‑88 (MQ=255)
      gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCa                                                                             >  1:174907/1‑90 (MQ=255)
         ccGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaa                                                                          <  1:274896/90‑1 (MQ=255)
          cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATaaa                                                                         >  1:5377/1‑89 (MQ=255)
          cGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATaaa                                                                         >  1:208125/1‑89 (MQ=255)
           gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATaaca                                                                        >  1:157538/1‑90 (MQ=255)
               aCAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaaca                                                                    <  2:116817/90‑1 (MQ=255)
                                  aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATAAAAAAAAACTCCCTTTCCGGAGAAcc                                                >  1:339466/1‑90 (MQ=255)
                                    aCAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAAAAAAAAACTCCTTGTCCGGAGAAaccc                                              >  1:360856/1‑90 (MQ=255)
                                       aaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGgaga                                                   >  1:127632/1‑82 (MQ=255)
                                       aaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGgaga                                                   <  2:127632/82‑1 (MQ=255)
                                             aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAAAAAAAAAATCCTTGTCCGGGGAAACCCCACggggga                                     >  2:44052/1‑88 (MQ=255)
                                             aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCAAAAAAAAAAACACCCTTTCCGGGGGAACCCCAAAGGGGa                                      >  2:422229/1‑90 (MQ=255)
                                                  aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  2:237753/90‑1 (MQ=255)
                                                  aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  2:313663/90‑1 (MQ=255)
                                                  aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  1:149515/90‑1 (MQ=255)
                                                                                gTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc  <  1:293377/90‑1 (MQ=255)
                                                                                gTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc  <  2:157538/90‑1 (MQ=255)
                                                                                 |                                                                                         
GGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC  >  NZ_CP009273/4558240‑4558409

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