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.

TGTTCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT  >  KT2440/1537280‑1537541
                                                                                                                                        |                                                                                                                                
tGTTCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGgcggcg                                                                                                                                 >  1:121861/1‑137 (MQ=255)
tGTTCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGgcggcg                                                                                                                                 >  5:120773/1‑137 (MQ=255)
          cGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCgg                                                                                                                       >  1:236074/1‑137 (MQ=255)
          cGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCgg                                                                                                                       >  3:155933/1‑137 (MQ=255)
               agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCg                                                                                                                 >  7:366853/1‑137 (MQ=255)
               agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCg                                                                                                                 >  1:386919/1‑137 (MQ=255)
                      ggggCTGGTTTCGGCCATGATGGTGATCAGACCGATGGCCTTGGTGGGTTTGGGCGGCGCGCTGCTGGCGATGATTGGCGCGGTTGCCAGTTTGATGAATTTTGCCTTGGCGGGGGGGGTGGCGGGGGCGTTGgcgc                                                                                                         <  2:379277/137‑1 (MQ=255)
                              tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCg                                                                                                 >  6:278625/1‑137 (MQ=255)
                                     cATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGt                                                                                          >  8:165596/1‑137 (MQ=255)
                                      aTGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCAATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGtt                                                                                         >  8:273184/1‑137 (MQ=255)
                                          tGGGGATGATCACGATAGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTtgtg                                                                                     >  1:28891/1‑137 (MQ=255)
                                               atgatCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTggtggg                                                                                >  6:224176/1‑137 (MQ=255)
                                                    cACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggc                                                                           >  4:186386/1‑137 (MQ=255)
                                                              ttGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGa                                                                 >  1:38001/1‑137 (MQ=255)
                                                                    gggTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAggggg                                                           >  1:51093/1‑136 (MQ=255)
                                                                                                                   gctgatttttGCCGGGGGGGGGGGGGGGGGGGGGGCGTTGGCGCGGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  2:3164/131‑1 (MQ=255)
                                                                                                                    tttttttttGCCCGGGGGGGGGGGGGGGGGGGGGGCGTGGGCGGGGTGGCGGGTGTGGTTTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  5:278625/132‑1 (MQ=255)
                                                                                                                    ggtgcttttGCCCGGGGGGGGGGGGGGGGGGGGGGGGTTGGCGCGGTGGCCGCGGTGGTTTGGGGGGGGGGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  7:320331/132‑1 (MQ=255)
                                                                                                                    ggtgaTTTTGCCCGGGGGGGGGGGGGGGGGGGGGGGGTGGGCGCGGGGGCGGGGGTGGTTTGGGGGGGGGGGGGGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  3:269442/133‑1 (MQ=255)
                                                                                                                    ggggaTTTTGCCCCGGGGGGGGGGGGGGGGGGGGGGGTGGCGGCGGGGGCGGCTGTGGTTTGTGGTGGGGGGCCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa            <  1:367271/133‑1 (MQ=255)
                                                                                                                         ttttGCCCGGGGGGGGGGGGGGGGGGGGGGCGGTGGCGCGGTGGCCGGTGTGGTTTGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCg       <  4:8757/137‑1 (MQ=255)
                                                                                                                         ttttGCCCGGGGGGGGGGGGGGGCGGGGGCCGTTGGCGCGGTGGCCGCTGTGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCg       <  6:161346/137‑1 (MQ=255)
                                                                                                                              cccgtGGGGGCGGGTGGGGGGGGGGCGTTGGCGGTGGGGCGGGTGTGGGTGGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGGAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCt  <  7:335247/133‑1 (MQ=25)
                                                                                                                                        |                                                                                                                                
TGTTCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT  >  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: