I have two Chip seq bed files with genomic coordinates one is from control and second is from pull down (treated). These files are from GEO and no other information is available. I want to identify unique non-overlapping coordinates in second (treated) bed file as compared to first file. I know I can use intersect, Bed tool to find overlapping regions. How should I find out non-overlapping unique regions enriched in treated file.
Thanks
•
written
by
kanwarjag
1 answer
One option is to use the -v argument of bedtools intersect which will return peaks that have no overlapping regions with the control peaks. Your second option would be to use bedtools subtract to subtract out from the experimental peaks parts that overlap with the control peaks.
Taking a step back for a second, usually a GEO entry will have an associated SRA accession, in which case you could download the fastq files, processes them, and then use MACS2 to properly call the peaks in the experimental condition using the control.
Log in to answer this question.
-
How to find overlapping regions in a single bed file
written by AA •
-
Histone marks enrichment analysis
written by gspirito •
-
Minimum overlap when comparing BED files
written by pablo.anton.garcia •
-
Genomics regions that common to atleast 10 cell types
written by rohitsatyam102
-
Find unique non overlapping intervals between two files?
written by a.rex
-
Find consensus overlap between many bed files
written by goodez
-
Combine regions of an interaction matrix
written by bharath.mh24 •
-
Fraction of bases covered by each element in bed
written by avinash.sahu •
-
Motif And Down Stream Analysis
written by kanwarjag
-
Chip-Seq Visualization Without Using Any Peak Caller
written by kanwarjag
