Closed hargup closed 9 years ago
Won't work if transcript on minus strand. Won't work if splice sequence is lowercase or mixedcase, since you only compare against uppercase GT and AG. In markNonCanonicalSpliceSites
you are marking the entire intron, instead of just the splice site. See - https://github.com/yeban/afra/commit/47f27ae43998f3a5e1270c1b6fbe1a854cc59f4d.
This builds on the work of @yeban.