Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092734,558,3390CA57.1% 8.0 / 17.1 14intergenic (‑3/‑476)mdtM/BW25113_RS25845multidrug efflux MFS transporter MdtM/Rpn family recombination‑promoting nuclease/putative transposase
Reads supporting (aligned to +/- strand):  ref base C (1/5);  new base A (8/0);  total (9/5)
Fisher's exact test for biased strand distribution p-value = 3.00e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-03
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

ACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGA  >  NZ_CP009273/4558255‑4558404
                                                                                    |                                                                 
aCAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATaacaaaca                                                              <  1:409481/90‑1 (MQ=255)
 cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAacaaacaa                                                             <  1:333508/90‑1 (MQ=255)
 cAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCATAaaaaaaaa                                                             >  1:704758/1‑90 (MQ=255)
      aCGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGGGGCCAAAAAAAAAACTCCt                                                        >  2:604499/1‑90 (MQ=255)
             aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAACCGGGCATAAAAAACAACACCCTTTcccg                                                 >  2:521722/1‑88 (MQ=255)
             aaaGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAAAGGCAAAAAAAAAAAACCCTTTTcccg                                                 >  1:687496/1‑88 (MQ=255)
                          aTCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAAAGGGGCATAAAAAAAAACTCCTTTTCCCGGGAACCCCAACgg                                    >  2:334998/1‑90 (MQ=255)
                              aCGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGCATAAAAAAAAAATCCCTGTCCGGAGAACCCCCACggggga                                >  2:367192/1‑88 (MQ=255)
                               cGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAAAAAAAAACCCCTTTTCCGGGGAACCCCCACGGGGgaa                               >  2:481353/1‑90 (MQ=255)
                                   caggggAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                           <  2:485255/87‑1 (MQ=255)
                                   aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                           <  2:105725/90‑1 (MQ=255)
                                   aTCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGa                           <  1:363031/90‑1 (MQ=255)
                                     cGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGaaa                         <  1:555108/90‑1 (MQ=255)
                                              cAGCGTGGCGGCATGGCGGGGAAAAAAAAGTGGGATAAAAAAAAACTCCTTTTCCGGGGAAACCCCACGGGGGAAAAGAAAATCaaaaac                >  2:194270/1‑86 (MQ=255)
                                                 cGTGGCGGCATGGCGGGTAAAAAAAAGGGGCATAAAAAAAAACTCCTTGTCCGGGGAACCCCCAACGGGAGAACGAAAGTCAATAACgac             >  1:727140/1‑88 (MQ=255)
                                                 cGTGGCGGCATGGCGGGGAAAAAAAAGGGGAGGAACAAACAAAATCTTGTACCGGGAAACCCCACGGGGGAAAAGAAAATTAATACCgcg             >  2:582204/1‑88 (MQ=255)
                                                            ggCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGa  <  1:367115/90‑1 (MQ=255)
                                                                                    |                                                                 
ACAGATACGCAGCAAAGTCATACAAAATCAACGCCATCGGGAAAAACAGCGTGGCGGCATGGCGGGTAAAAAAACGTGGCATAACAAACAACTCCTTGTCCGGAGAACCCCAACGGGGAAAACGAAATTCAATACCGTCGGTAATGACGA  >  NZ_CP009273/4558255‑4558404

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