Splitting paired-end RNA-seq data in bedgraph format according to strand
0
0
Entering edit mode
6.6 years ago

I currently have several paired-end RNA-seq datasets scaled to rpm in bedgraph format (thanks to a former biostar question post reply that was very helpful! :)).

Is it possible to work from these files as they are, and split them into two tracks, one with only positive strand reads and one with only negative strand reads? Or is it necessary to return to the original bam files to do this? Either way, how might the command lines look, and how will the fact that these are paired-end datasets influence the commands I have to use?

RNA-Seq • 1.3k views
ADD COMMENT
0
Entering edit mode

Could you show the first 5 rows of your bedgraph-formatted files?

ADD REPLY

Login before adding your answer.

Traffic: 1383 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6