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

GGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC  >  NZ_CP009273/4558240‑4558409
                                                                                 |                                                                                        
ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGgaaaaaaa                                                                                  >  1:156652/1‑89 (MQ=255)
ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa                                                                                  >  1:224032/1‑89 (MQ=255)
ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa                                                                                  >  2:276533/1‑89 (MQ=255)
 gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGgaaaaaaaa                                                                                 >  1:19327/1‑88 (MQ=255)
 gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaaa                                                                                 >  1:88301/1‑88 (MQ=255)
     aGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGc                                                                             <  1:502790/90‑1 (MQ=255)
      gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGTgggg                                                                            >  1:571675/1‑88 (MQ=255)
      gATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAATGGCa                                                                            >  1:89637/1‑90 (MQ=255)
               aCAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaaca                                                                   <  1:276533/90‑1 (MQ=255)
               aCAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaaca                                                                   <  2:93740/90‑1 (MQ=255)
                cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAacaaacaa                                                                  >  2:200570/1‑90 (MQ=255)
                     aCGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAAATCCt                                                             >  2:191625/1‑90 (MQ=255)
                     aCGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAACAAAAAACTCCt                                                             >  2:362470/1‑90 (MQ=255)
                         agcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAACAAAAAACTCCTTGTc                                                         >  2:35070/1‑90 (MQ=255)
                          gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGGATAAAAAAAAAAACCTTGTcc                                                        >  1:191566/1‑90 (MQ=255)
                                  aTACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCAAAAAAAAAAACCCCTTGTCCGGGGaaac                                                >  1:390426/1‑88 (MQ=255)
                                    aCAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACGTGGCATAACAAAAAACTCCTTGTCCGGAGAAcccc                                              >  2:222009/1‑90 (MQ=255)
                                      aaaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCaa                                            <  2:106726/90‑1 (MQ=255)
                                       aaaTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAGGGCATAAAAAAAAAAACCCTTTCCGGGGAAACCCAAc                                           >  1:93509/1‑90 (MQ=255)
                                                  aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  2:390426/90‑1 (MQ=255)
                                                  aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  2:89637/90‑1 (MQ=255)
                                                  aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                                <  1:306338/90‑1 (MQ=255)
                                                  aTCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAACAAATCCTTGTCCCGAGAACCCCAACGGGGGAAAAGa                                >  2:306338/1‑90 (MQ=255)
                                                                         aTGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGAc         <  2:395592/90‑1 (MQ=255)
                                                                            gCGGGGAAAAAAAAGGGGCATAACAAACAACACCTTGTCCGGGGAACCCCCAAGGGGGAAAAGAAAATCAAAAACATCGGGGATGaagga      >  1:542046/1‑86 (MQ=255)
                                                                                gTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc  <  2:80220/90‑1 (MQ=255)
                                                                                gTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGc  <  2:88301/90‑1 (MQ=255)
                                                                                 |                                                                                        
GGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGAAAGGC  >  NZ_CP009273/4558240‑4558409

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