Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
50 bp paired end reads vs. 100 bp single end reads efoss Bioinformatics 12 01-15-2014 08:05 PM
Mapping paired-end stranded RNA-seq data hubin.keio Bioinformatics 3 03-28-2013 02:10 PM
Can Cuffdiff treat paired-end and single-end reads at the same time? zun RNA Sequencing 3 06-12-2012 05:37 PM
paired-end reads mapped to genome.. gene with only one direction of paired-end reads? danwiththeplan Bioinformatics 2 09-22-2011 02:06 AM
Method to get BigWig of Paired-End Data gibsongenetics Bioinformatics 0 09-19-2011 09:25 AM

Thread Tools
Old 03-24-2014, 01:17 PM   #1
Senior Member
Location: USA

Join Date: Sep 2012
Posts: 130
Default Stranded BedGraph/BigWig for paired-end reads

For stranded data, I use "bedtools genomecov" with -strand to generate coverage plots. For single reads, this works fine, but for paired-end reads, it doesn't actually take mate pairs into account. Thus, half the data (reverse reads) ends up on the "wrong" strand. Is there a way to generate coverage plots for paired-end stranded data?

Is the only solution to split the BAMs, process them separately, and then merge? That seems like a lot of work. Also, very heavy on the I/O.
id0 is offline   Reply With Quote
Old 08-23-2014, 06:06 PM   #2
Location: Canada

Join Date: Dec 2013
Posts: 62

I am about to start RNA-seq and I'm wondering the same thing. When it comes time to generate coverage files from the aligned sam/bam files, will bedtools genomecov (using the -split option) correctly process paired-end sequencing? If not is there another method to generate these coverage files?
biocomputer is offline   Reply With Quote
Old 08-26-2014, 12:27 AM   #3

Join Date: Jun 2010
Posts: 81

you can use igvtools count --strands first for creating a wig file that contains columns for both strands. subsequently use cut/awk and wigToBigWig.
volks is offline   Reply With Quote
Old 06-22-2017, 02:13 PM   #4
Junior Member
Location: CA

Join Date: Jul 2009
Posts: 8

since this is one of the top google hits for this problem.
let me plug

deeptools (the package)
bamCoverage (the command)

it seems to do it without the gd flags and without all the awk

Last edited by yoyoq; 06-22-2017 at 02:18 PM.
yoyoq is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 01:51 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2019, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO