Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*KT24401,537,4150CG38.1% 20.4 / 11.4 21A388G (GCG→GGG) mdfAMultidrug transporter MdfA
Reads supporting (aligned to +/- strand):  ref base C (12/1);  new base G (0/8);  total (12/9)
Fisher's exact test for biased strand distribution p-value = 4.42e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTTCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT  >  KT2440/1537280‑1537541
                                                                                                                                        |                                                                                                                                
tGTTCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGgcggcg                                                                                                                                 >  1:121861/1‑137 (MQ=255)
tGTTCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGgcggcg                                                                                                                                 >  5:120773/1‑137 (MQ=255)
          cGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCgg                                                                                                                       >  1:236074/1‑137 (MQ=255)
          cGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCgg                                                                                                                       >  3:155933/1‑137 (MQ=255)
               agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑Cg                                                                                                                 >  7:366853/1‑137 (MQ=255)
               agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑Cg                                                                                                                 >  1:386919/1‑137 (MQ=255)
                      ggggCTGGTTTCGGCCATGATGGTGATCAGACCGATGGCCTTGGTGGGTTTGGGCGGCGCGCTGCTGGCGATGATTGGCGCGGTTGCCAGTTTGATGAATTTTGCCT‑TGGCGGGGGGGGTGGCGGGGG‑CGTTGgc‑gc                                                                                                         <  2:379277/137‑1 (MQ=255)
                              tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCg                                                                                                 >  6:278625/1‑137 (MQ=255)
                                     cATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCGCTGTGGt                                                                                          >  8:165596/1‑137 (MQ=255)
                                      aTGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCAATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCGCTGTGGtt                                                                                         >  8:273184/1‑137 (MQ=255)
                                          tGGGGATGATCACGATAGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCGCTGTGGTTtgtg                                                                                     >  1:28891/1‑137 (MQ=255)
                                               atgatCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCGCTGTGGTTTGTggtggg                                                                                >  6:224176/1‑137 (MQ=255)
                                                    cACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCGCTGTGGTTTGTGGTGGGTGggc                                                                           >  4:186386/1‑137 (MQ=255)
                                                              ttGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGa                                                                 >  1:38001/1‑137 (MQ=255)
                                                                    gggTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAggggg                                                           >  1:51093/1‑136 (MQ=255)
                                                                                                                   gctgatttttGCCG‑GGGGGGGGGGGGGGGGGGGGG‑CGTTGGC‑GCGGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  2:3164/131‑1 (MQ=255)
                                                                                                                    tttttttttGCCCGGGGGGGGGGGGGGGGGGGGGG‑CGTGGGC‑GGGGTGGCGGGTGTGGTTTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  5:278625/132‑1 (MQ=255)
                                                                                                                    ggtgcttttGCCCGGGGGGGGGGGGGGGGGGGGGG‑GGTTGGC‑GCGGTGGCCGCGGTGGTTTGGGGGGGGGGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  7:320331/132‑1 (MQ=255)
                                                                                                                    ggtgaTTTTGCCCGGGGGGGGGGGGGGGGGGGGGG‑GGTGGGC‑GCGGGGGCGGGGGTGGTTTGGGGGGGGGGGGGGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  3:269442/133‑1 (MQ=255)
                                                                                                                    ggggaTTTTGCCC‑CGGGGGGGGGGGGGGGGGGGG‑GGGTGGCGGCGGGGGCGGCTGTGGTTTGTGGTGGGGGGCCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  1:367271/133‑1 (MQ=255)
                                                                                                                         ttttGCCCGGGGGGGGGGGGGGGGGGGGGG‑CGGTGGC‑GCGGTGGCCGGTGTGGTTTGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCg       <  4:8757/137‑1 (MQ=255)
                                                                                                                         ttttGCCC‑GGGGGGGGGGGGGGGCGGGGGCCGTTGGC‑GCGGTGGCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCg       <  6:161346/137‑1 (MQ=255)
                                                                                                                              cccgtGGGGGCGGGTGGGGGGGGGG‑CGTTGGC‑GGTGGGGCGGGTGTGGGTGGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGGAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCt  <  7:335247/133‑1 (MQ=25)
                                                                                                                                        |                                                                                                                                
TGTTCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑TGGCGGCGGGGGTGGCGGGGG‑CGTTGGC‑GCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT  >  KT2440/1537280‑1537541

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑