Question: how can I get the coverage over introns only in a RNA-seq experiment?
0
gravatar for Lila M
4 days ago by
Lila M 390
UK
Lila M 390 wrote:

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!

rna-seq coverage introns • 75 views
ADD COMMENTlink written 4 days ago by Lila M 390

There are plenty of posts out there about getting exon coverage. You can do the same just with introns as reference. I am not sure if UCSC has intron lists ready, but you can follow my tutorial to get them from any gtf/gff file.

ADD REPLYlink written 4 days ago by ATpoint9.1k

Do you want to create coverage plots or count the number of reads mapping to exons?

ADD REPLYlink written 3 days ago by kristoffer.vittingseerup700

both , coverage plot and number of reads mapping to INTRONS will be nice :)

ADD REPLYlink written 3 days ago by Lila M 390
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1528 users visited in the last hour