| Read alignment evidence... | |||||||||||
|---|---|---|---|---|---|---|---|---|---|---|---|
| seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
| * | ACB143 | 1,537,415 | 0 | C | G | 34.8% | 25.6 / 14.4 | 25 | A388G (GCG→GGG) | mdfA | Multidrug transporter MdfA |
| Reads supporting (aligned to +/- strand): ref base C (13/2); new base G (0/8); total (13/12) | |||||||||||
| Fisher's exact test for biased strand distribution p-value = 9.18e-05 | |||||||||||
| Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00 | |||||||||||
| Rejected as consensus: Frequency below/above cutoff threshold. | |||||||||||
TCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTT‑TGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT > ACB143/1537283‑1537541 | tctcCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑Gc‑ggcgg > 6:366106/1‑135 (MQ=255) agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACTATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGc > 3:354315/1‑136 (MQ=255) agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCggggg > 7:375882/1‑135 (MQ=255) aaGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTgg > 3:356301/1‑136 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGcc > 4:160548/1‑136 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGc > 1:225210/1‑135 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGc > 7:117197/1‑135 (MQ=255) tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGc > 6:94189/1‑135 (MQ=255) ttCGGCCTTGCTGGGGATGATAACGATTGCCGTGCTGGGTTTTGTCGGCGCGCTGCTGGCTATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGcc < 2:100031/135‑1 (MQ=255) ttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGcc > 1:148704/1‑135 (MQ=255) tCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCGGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGcc < 8:205379/134‑1 (MQ=255) aTGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGt > 4:84564/1‑136 (MQ=255) tGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTT‑tgt > 3:276043/1‑136 (MQ=255) aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTT‑TGTGGtgggtggg > 5:395363/1‑135 (MQ=255) aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTT‑TGTGGTGGGTGggc > 3:380139/1‑136 (MQ=255) gATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATCGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTT‑TGTGGTGGGTGggcgg > 3:6925/1‑136 (MQ=255) ccccgggTGCCCGGTTGTGGCTTTTGG‑CCG‑G‑G‑GC‑GGGGGGGGGGGGGGGGGCTTTGGCGGTGTGGCGGCTTTGGTT‑TGGGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCgggg < 8:119256/129‑1 (MQ=255) tgatttttG‑CCC‑G‑G‑GG‑GGGGGGGGGGGGGGGGGCGTTGGCGCGGTGGCCGCTGTGGTT‑TGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 1:256620/131‑1 (MQ=255) tgatttttC‑CCC‑T‑GCGC‑GG‑GGGGGGGGCGGGGCCGAGGGCGCGGGGGCCGCTGTGGTT‑GGGGGTGGGGGGGCGGCCTGGGGACGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 8:176466/131‑1 (MQ=255) gtgaTTTTG‑GCC‑TGG‑GG‑GGGGGGGGGGGGGGGGGGGTGGGCGC‑GGGGGCGGTGGGGTTGTGGGGTGGGGGGGCGGCCTGGGGACGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 3:68966/133‑1 (MQ=255) cgcATTTGG‑CCCGG‑G‑GG‑GGGGGGGGGGGGGGGGGCGTGGGCGCGGGGGCCGGTGGGTTT‑TGGGGGGGGGGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 3:235131/135‑1 (MQ=255) tttttttG‑CCG‑G‑G‑GCGGGGGGGGGGGGGGGGGGGGTGGGCGGTGTGGCCGGTGTGTTT‑TGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 8:268409/132‑1 (MQ=255) ttttttgG‑CGC‑T‑G‑GG‑GGTGGGGGGGGGGGGGGGGGTGGGGCTGTGTCCGCGGTGTTT‑TGTGGGTGGTGGGTGGCGGGGGTAAGGGTCTGTAGCTCGATAATCGGTATTGCAGGGGCCGTCCTCTTCGCGTGTaa < 2:365952/131‑1 (MQ=255) tgcttttGCCCC‑G‑G‑GG‑GGGGGGGGGGGGGGGGGCGTTGGGGCTGGGGCGGGTGTGGTT‑TGTGGTGGGGGGGCGGCGGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 6:112917/132‑1 (MQ=255) gctttttG‑CCCGG‑G‑GG‑GGGGGGGGGGGGGGGGGGGTTGGCGGGGGGGCGGGTGGGGTT‑GGGGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa < 2:365957/132‑1 (MQ=255) ccgt‑G‑GG‑GGGGGGGGGGGGGGGGGGGTTGGCGCTGGGGCCGCTGTGTTT‑TGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCt < 4:356768/133‑1 (MQ=37) | TCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CCC‑T‑G‑GC‑GGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTT‑TGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT > ACB143/1537283‑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: ‑ |