Hi everybody, I've just wondering the best way to get the coverage over introns (only) in RNA-seq experiments (human genome) My first approach is to download the complete intron list from UCSC. Then I could get the overlaps (using bedtools intersect or samithing similar) with each bed file. Once I get that, I may create a coverage profile using computeMatrix (deepTools). But I'm not 100% if this is totally right. I appreciate if someone has done this before and would like to share! Thanks!