Question: bedtools coverage and paired-end data
gravatar for eric.londin
8 weeks ago by
eric.londin20 wrote:

I have a question using bedtools with paired-end stranded RNA-seq data. Since one read is on the forward direction, and the second read is in the reverse direction, how does bedtools coverage handle these paired-reads when running it with the "-s" flag. For example, I am interested in getting coverage of a gene in a strand specific manner. Does bedtools using the paired read information. from what I can tell, it is handling each read separately, and only considering the read that is on the same strand as the gene of interest. Am I correct in assuming this? If so, then it would only count half the reads

for clarity, here is my command line that I have run:

./bedtools coverage -a MyGenes.bed -b MyMappings.bam -s -split > outfile.txt

Any thoughts would be appreciated.

rna-seq paired-end bedtools • 94 views
ADD COMMENTlink written 8 weeks ago by eric.londin20
Please log in to add an answer.


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