View Single Post
Old 01-22-2015, 11:38 PM   #2
sarvidsson
Senior Member
 
Location: Berlin, Germany

Join Date: Jan 2015
Posts: 137
Default

You can extract the overlapping coverage data with bedtools intersect, e.g.

Code:
bedtools intersect -wb -a coverage.bdg -b transcripts.bed
(-wb is helpful to see which transcript in the "b" file overlapped the features in the "a" file - you could add transcript identifiers to the "b" for easier parsing)

But if you are actually interested in the actual number of reads overlapping with transcript features, you're probably better of using the original BAM with a GFF3/GTF file and HTSeq-count (http://www-huber.embl.de/users/ander...doc/count.html) - it has some further interesting features for RNA-Seq (if that's what you're into).
sarvidsson is offline   Reply With Quote