| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | KT2440 | 1,537,415 | 0 | C | G | 34.3% | 49.8 / 18.4 | 36 | A388G (GCG→GGG) | mdfA | Multidrug transporter MdfA |
| Reads supporting (aligned to +/- strand): ref base C (22/1); new base G (0/12); total (23/13) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 1.56e-08 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTC > KT2440/1537285‑1537539 | tcCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑Gt > 1:266504/1‑137 (MQ=255) cGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑gg > 1:118155/1‑137 (MQ=255) cGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑gg > 5:207717/1‑137 (MQ=255) cGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑gg > 6:443476/1‑137 (MQ=255) agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCg > 2:320921/1‑137 (MQ=255) agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCg > 7:213118/1‑137 (MQ=255) agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCg > 7:276473/1‑137 (MQ=255) ggggCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGgcgc > 7:351181/1‑137 (MQ=255) gTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGcc > 3:413996/1‑137 (MQ=255) tttCGGCCATGCTTGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCg > 3:230557/1‑137 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCg > 5:388893/1‑137 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCg > 5:289695/1‑137 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCg > 4:430440/1‑137 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGACGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCg > 5:351073/1‑137 (MQ=255) ggggtttgatCACGTGTGCGTGGCTGGGGTTGGGCGGGCGGGTGGTGGC‑TTGGTTGGGCGCGGTTCCACCGTTGC‑GCATTTTGGCCG‑GGGGGGGGGGG‑GGGGG‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTtgtg < 8:188107/131‑1 (MQ=255) ttctcccGATGGCCTTGGTGGGGTTTGGCGGCGCGCGGCGGGCTATGATTGGCGCGGTGGCCTGC‑TTGTTGCTTTTT‑GCCG‑TGGGGGGGGGG‑GTGGC‑GGGGGCGGTGGCGCTGTGGCCGCTGTGGTTTGTtgtgggt < 8:337885/131‑1 (MQ=255) cACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggc > 2:369866/1‑137 (MQ=255) cACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggc > 8:240800/1‑137 (MQ=255) aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggcg > 7:298076/1‑137 (MQ=255) aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGAGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGGAGTGGGTTGTTGTTGGTGggcc > 1:224678/1‑136 (MQ=255) ttGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCt > 7:56171/1‑137 (MQ=255) tGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAggg > 3:10766/1‑137 (MQ=255) tGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAg > 6:307954/1‑137 (MQ=255) gctgctGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTATGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAAt > 4:99021/1‑137 (MQ=255) gcttattttt‑GCCT‑GGGCGGGGGGG‑GGGGG‑GGGGGGTTTGGCGGGGTGGCGGGTGTGGTTTGTGGTGGGTGGGGGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑Taa < 6:351073/131‑1 (MQ=255) gctgcATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GAGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑Taa < 3:363987/137‑1 (MQ=255) gctccATTTT‑GCCG‑GGGGGGGGGGG‑GGGGG‑GGGGCCGTTGGCGCGGTGGCGGCTGTGGTTTGGGGTGGGGGGGCGGCCTGGGGACGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑Taa < 1:91280/137‑1 (MQ=255) gctccATTTG‑GCCC‑TGGGGGGGGGG‑GGGGC‑GGGGGGGTTGGCGCTGTGGCGGGTGGTGTTTGTGGTGGGGGGCCGGCCTGGGGAAGGGGCTGTAGCTCTATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑Taa < 7:370699/137‑1 (MQ=255) gcggcttttt‑CCCC‑GGGGGGGGGGG‑GGGGG‑GGGGGCGTGGGCGCGGGGGCCGCTGTGGTTTGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑Taa < 3:252564/131‑1 (MQ=255) gtttaTTTT‑GCCCGGGGGGGGGGGG‑GGGGG‑GGGGGGGTGGGCGCGGGGGCCGCTGTGGTTTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑Taa < 7:240800/133‑1 (MQ=255) ggttctttt‑GCCCGGGGGGGGGGGG‑GGGGG‑GGGGGCGTGGGCGCTGTGGCGGGTGTGGTTTGGGGGGGGTGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑Taa < 4:413996/132‑1 (MQ=255) ggtcttttt‑GCCCGTGGGGGGGGGG‑GGGGG‑GGGGGGGTTGGGGCTGTGGCGGGTGTGGTTTGTGGGGGGTGGGGGGCCTGGGTAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑Taa < 8:56171/132‑1 (MQ=255) gggcATTTT‑GCCCGGGGGGGGGGGG‑GGGGG‑GGGGGGGGTGGCGCGGGGGCGGGTGGGGTTGGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑Taa < 1:369866/135‑1 (MQ=255) cctgtggCCC‑TGGGGGGGGGG‑GGGGCGGGGGGGGTGGGCGTGGTGGCGGCGGGGGTGTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑TAAAcc < 8:159358/131‑1 (MQ=255) tttt‑GCCG‑GGGGGGGGGGGTGGGGG‑GGGGGCGTGGGCGGTGGGGCCGGGGGGGTGGGGGGGGGGGGGGCGGCGGGGGGACGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑TAAACCCg < 5:73867/137‑1 (MQ=255) t‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTTAAAccccaac > 2:347148/1‑133 (MQ=255) | TCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGC‑TTGCTGCATTTT‑GCCC‑TGGCGGCGGGG‑GTGGC‑GGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGG‑TAAACCCGCTC > KT2440/1537285‑1537539 |
| 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: ‑ |