This is a test version of Biostars. For the public version, visit https://www.biostars.org.
checking overlap between two genomic coordinates

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
• 3,694 views • link
updated by Tm Tm • written by kanwarjag kanwarjag

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 views • link
written by michael.j.apostolides • michael.j.apostolides

Hi, an answer with bedtools intersect already exists. If you have more to add, please add it as a comment to that answer so people going down that path can find all relevant information in the same thread.

• 0 views • link
written by Ram Ram

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 views • link
written by Tm Tm

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

• 0 views • link
written by GenoMax GenoMax

3 answers

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

• 0 views • link
written by Ram Ram

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

• 0 views • link
written by jordi.planells jordi.planells

Have a look at BedSect and Galaxy online platforms.

• 0 views • link
written by Shalu Jhanwar Shalu Jhanwar

Log in to answer this question.

More posts like this