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

CAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCG  >  NZ_CP009273/4558234‑4558394
                                                                                       |                                                                          
cAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaa                                                                          >  1:126409/1‑90 (MQ=255)
  ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa                                                                        >  2:402725/1‑90 (MQ=255)
  ggCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaa                                                                        >  2:380553/1‑90 (MQ=255)
      ggATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaa                                                                    >  2:326066/1‑89 (MQ=255)
       gATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGaaaaaaaaa                                                                   >  2:51650/1‑88 (MQ=255)
             aTCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCAt                                                             <  1:51650/90‑1 (MQ=255)
               ccGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaa                                                           <  1:318498/90‑1 (MQ=255)
                 gTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATaaaa                                                         >  1:80473/1‑88 (MQ=255)
                           aCGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCt                                               <  1:91266/90‑1 (MQ=255)
                           aCGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCt                                               <  1:128662/90‑1 (MQ=255)
                                gcaAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGAGGCCAAAAAAAAAACACCtttctc                                          >  1:195480/1‑86 (MQ=255)
                                               aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTGAAAAAAATTGGCCTAAAAAACAAACCCTTGTCCGGGGAACCCCCACgg                          >  2:217027/1‑90 (MQ=255)
                                               aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAATTTGCCCAACAAAAAAACCCTTTTCCCGCGGACCCCaaacg                          >  2:264593/1‑87 (MQ=255)
                                               aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGGATAAAAAAAAACTCCTTGTCCGGAGAACCCCCACgg                          >  2:317898/1‑90 (MQ=255)
                                                   aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAAAAAAAAACACCTTGTTCGGAGAAGCCCCACCGGGGa                       >  1:345874/1‑90 (MQ=255)
                                                        aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                 <  1:380553/90‑1 (MQ=255)
                                                        aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                 <  1:402725/90‑1 (MQ=255)
                                                        aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                 <  1:125900/90‑1 (MQ=255)
                                                          cGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGaaa               <  1:264729/90‑1 (MQ=255)
                                                                      cGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTc   <  2:345874/90‑1 (MQ=255)
                                                                       gTGGCGGCATGGCGGGGAAAAAAAAGGGGCATAAAAAAAAACCCCTTGTCCGGGGAAACCCCACGGGGGAAAAGAAAATAAAAAACGTCg  >  1:10802/1‑90 (MQ=255)
                                                                                       |                                                                          
CAGGCTGGATCAGATCCGTCGACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCG  >  NZ_CP009273/4558234‑4558394

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: