Question: Intersectbed: Return Reads In Fraction In Input Files
0
gravatar for anneleen.beckers
6.5 years ago by
anneleen.beckers70 wrote:

I have a question with respect to intersectBED and multiple input files:

Is it possible to return reads which are present in, say 8/10 input files, without fractioning the reads in smaller intervals?

Thank you

bedtools overlap intersect • 1.2k views
ADD COMMENTlink written 6.5 years ago by anneleen.beckers70
0
gravatar for Alex Reynolds
6.5 years ago by
Alex Reynolds27k
Seattle, WA USA
Alex Reynolds27k wrote:

The BEDOPS suite includes bedops, a utility which performs various set operations on BED inputs. It's not clear if this would solve your specific problem, but you can specify the degree of overlap as criteria for inclusion, when using the --element-of operator with bedops.

For example, the following command would return elements from Reads.bed which overlap elements in Map.bed by 80% or more of segment length:

$ bedops --element-of -80% Reads.bed Map.bed
ADD COMMENTlink written 6.5 years ago by Alex Reynolds27k
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: 754 users visited in the last hour