Question: bedtools coverage and paired-end data
0
gravatar for eric.londin
8 weeks ago by
eric.londin20
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.

Help
Access

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