Question: Beedtools intersect: more peaks intersected than intersected portions of those peaks?
0
gravatar for salamandra
9 months ago by
salamandra180
salamandra180 wrote:

When finding the number of peaks/genomic regions in A that intersects with B and C (note even the non-overlaping part of each peak is reported):

 bedtools intersect -f 0.5 -r -wa -a $A -b $B -sorted | sort -k1,1 -k2,2n | uniq > Apeaks.intersectingwith.B 
 bedtools intersect -f 0.5 -r -wa -a Apeaks.intersectingwith.B -b $C -sorted | sort -k1,1 -k2,2n | uniq > Apeaks.intersectingwith.B.and.C
 wc -l Apeaks.intersectingwith.B.and.C

we get 7 peaks/regions

But, when finding the number of regions that overlap in A, B and C (only the overlaping part of each peak is reported):

bedtools intersect -f 0.5 -r -a A.area.intersecting.B -b $C -sorted | sort -k1,1 -k2,2n | uniq > A.area.intersecting.B.and.C 
wc -l A.area.intersecting.B.and.C

we get 6 regions

How is it possible to have more peaks that intersect (in A), than intersecting areas of those peaks?? I don't get it..

ADD COMMENTlink modified 9 months ago • written 9 months ago by salamandra180
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: 928 users visited in the last hour