checking overlap between two genomic coordinates
3
0
Entering edit mode
15 months ago
kanwarjag ★ 1.1k

I have two different Chipseq datasets files with peak coordinates. and other information like gene name etc. I want check the overlap of coordinates in two files. Is there any online tool that can help me in checking the overlap between two files. Thanks

Kanwar

ChIP-Seq • 927 views
0
Entering edit mode

In my opinion the most clean way to do this is using bedtools intersect. This is the most standard and "old school" way, and can be easily integrated into pipelines with only one easy to install dependency (bedtools)

0
Entering edit mode

0
Entering edit mode

If you have gene annotation in gff/gtf format and reference genome fasta with you, then it is very easy to annotate peaks using annotatePeaks.pl, a perl script which is a part of HOMER package.

Command goes like this:

perl annotatePeaks.pl sample_peak.bed Reference_Genome.fasta -gft Reference_genome_annotation.gft   >homer_annotation.txt

0
Entering edit mode

Original question is not about annotating peaks. It is for finding overlap between two sets.

1
Entering edit mode
15 months ago
Ram 35k

findOverlaps from GenomicRanges might be of use. See: https://stackoverflow.com/questions/19101849/overlapping-genomic-ranges

1
Entering edit mode
15 months ago

You can also try bedtools intersect. See here for more details

0
Entering edit mode
15 months ago
Shalu Jhanwar ▴ 500

Have a look at BedSect and Galaxy online platforms.