Predicted mutation | |||||||
---|---|---|---|---|---|---|---|
evidence | seq id | position | mutation | freq | annotation | gene | description |
RA | minE | 1,729,738 | C→T | 17.5% | L372F (CTT→TTT) | yfhD → | predicted transglycosylase |
Read alignment evidence... | |||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|
seq id | position | ref | new | freq | score (cons/poly) | reads | annotation | genes | product | ||
* | minE | 1,729,738 | 0 | C | T | 17.5% | 129.3 / 19.4 | 58 | L372F (CTT→TTT) | yfhD | predicted transglycosylase |
Reads supporting (aligned to +/- strand): ref base C (20/27); new base T (6/4); total (27/31) | |||||||||||
Fisher's exact test for biased strand distribution p-value = 4.86e-01 | |||||||||||
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01 |
GCGCCCTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCTTA > minE/1729670‑1729805 | gcgcCCTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACtt < 1:3363953/71‑1 (MQ=255) gcgcCCTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACt > 1:971579/1‑70 (MQ=255) gcgcCCTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACt > 1:1810363/1‑70 (MQ=255) cTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGc > 1:559325/1‑69 (MQ=255) tGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAg > 1:3396385/1‑67 (MQ=255) tGACGGCAAAAACCAAAGGAAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGc > 1:576264/1‑68 (MQ=255) cGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGaaa > 1:2657924/1‑71 (MQ=255) ggCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAAc > 1:3133349/1‑71 (MQ=255) gCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAAcc < 1:2590229/71‑1 (MQ=255) cAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAAccc > 1:586978/1‑71 (MQ=255) cAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAAc > 1:3370808/1‑69 (MQ=255) aaaaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCt < 1:3387918/71‑1 (MQ=255) aaaaCCTAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCtttattta > 1:4322/1‑58 (MQ=255) aaaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTa > 1:3202518/1‑71 (MQ=255) aaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTa > 1:1921448/1‑70 (MQ=255) aaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTAt > 1:2280563/1‑71 (MQ=255) aaaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTAt > 1:2346819/1‑71 (MQ=255) aaCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTAATTAGCCAGaa > 1:520654/1‑63 (MQ=255) aCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTa < 1:1475850/71‑1 (MQ=255) aCCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTa > 1:3187880/1‑71 (MQ=255) aaaGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAg < 1:1443549/71‑1 (MQ=255) aaaGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAGCCCTATTACAg > 1:3593630/1‑71 (MQ=255) aaGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAAcc > 1:320257/1‑61 (MQ=255) aGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTATTTAGCCAGAAACCCTATTACAGCa < 1:1814251/71‑1 (MQ=255) ggAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCtttattta < 1:2823997/48‑1 (MQ=255) ggAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTa < 1:2013396/64‑1 (MQ=255) aaTCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCt < 1:3540318/71‑1 (MQ=255) aaTCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCt < 1:2666413/71‑1 (MQ=255) aTCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCt > 1:2927594/1‑70 (MQ=255) aTCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTg < 1:1045750/71‑1 (MQ=255) tGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGaa < 1:1386099/48‑1 (MQ=255) aCAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACtt < 1:2194069/69‑1 (MQ=255) aCAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACtt < 1:1357992/69‑1 (MQ=255) ttGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGAc > 1:1353495/1‑63 (MQ=255) ttGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACg > 1:717441/1‑68 (MQ=255) ttGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGc > 1:3751542/1‑70 (MQ=255) gggCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTa < 1:1263801/70‑1 (MQ=255) tGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTc > 1:1722216/1‑71 (MQ=255) gACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCg > 1:2137520/1‑71 (MQ=255) aCGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCGAGCTGACTTACGGCTAc > 1:1782010/1‑65 (MQ=255) aCGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCt < 1:2517391/52‑1 (MQ=255) aCGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCt < 1:2207265/52‑1 (MQ=255) aCGTAAAACAGCATCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTc < 1:2554697/69‑1 (MQ=255) gTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTAt < 1:360380/39‑1 (MQ=255) tAAAACAGCGTCTGCCTTTACTTAGCCAGAAAccc < 1:1002962/35‑1 (MQ=255) tAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGa < 1:91322/71‑1 (MQ=255) aaaaCAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGa > 1:2802000/1‑70 (MQ=255) aaaCAGCGTCTGCCTTTACTTAGCCAGAAACCCTAtt < 1:2861832/37‑1 (MQ=255) aaaCAGCGTCTGCCTTTACTTAGCCAGAAACCCTAtt < 1:2691227/37‑1 (MQ=255) aaCAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCg < 1:2172821/64‑1 (MQ=255) aaCAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACAt < 1:614596/71‑1 (MQ=255) aCAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCg < 1:1112959/63‑1 (MQ=255) aCAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATg < 1:34255/71‑1 (MQ=255) aCAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATg < 1:3650280/71‑1 (MQ=255) gTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCaa < 1:2623503/37‑1 (MQ=255) gTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGcc < 1:3264689/71‑1 (MQ=255) tGCCTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGaa > 1:1125561/1‑65 (MQ=255) cccttACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTAc < 1:1254576/66‑1 (MQ=255) ccTTTATTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGt < 1:1413838/54‑1 (MQ=255) ttACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCTTa > 1:3029219/1‑71 (MQ=255) | GCGCCCTGACGGCAAAAACCAAAGGGAATCCTGACAGTTGGGCTGACGTAAAACAGCGTCTGCCTTTACTTAGCCAGAAACCCTATTACAGCAAGCTGACTTACGGCTACGCTCGTGGACATGAAGCCTACGCTTA > minE/1729670‑1729805 |
Alignment Legend |
---|
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG |
Unaligned base: atcg Masked matching base: atcg Alignment gap: ‑ Deleted base: ‑ |