Question: Intersecting unbalanced BED files, normalizing for one
gravatar for LacquerHed
6 months ago by
LacquerHed0 wrote:

Have two bed files I am intersecting, one of approximately 500 promoters +/- 2kb of TSS and the other some locations of histone marks within those promoters. There are many more intervals in the latter dataset, I would like to get the number of intersects that are unique to the promoter interval file, meaning that even if 25 intervals from the histone mark fall within one promoter I would like it to report 1 intersect. Maybe there is some way of merging across those 25 and making one unified interval that stretches across them. Not really sure of a good way to go about this, thanks.

bedops bedtools • 159 views
ADD COMMENTlink written 6 months ago by LacquerHed0

you can use bedtools intersect with -wa option to report the intersected promoter for each interval, then use sort -u to get the unique promoter entry

ADD REPLYlink written 6 months ago by yztxwd380
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: 2443 users visited in the last hour