How to find overlapping coordinates of a gtf file
I am trying to list the entries of a gtf file (gencode.vM25.annotation.gtf) if they have overlap. I have been looking around for some tool to perform it - is there such a tool to do what I want? Thanks
• 3,804 views
•
link
updated
by
Alex Reynolds
•
written
by
Apex92 •
•
written
by
Apex92 •
2 answers
I assume you mean features in a GFF file that overlap other features in that same GFF file?
If so, have a look at AGAT , it certainly has some sub-programs that can do this.
here is more info: AGAT - Another Gff/Gtf Analysis Toolkit
• 0 views
•
link
written
by
lieven.sterck
You could use BEDOPS gtf2bed and bedmap to map entries to themselves, filtering out any that are disjoint with awk and cut:
$ gtf2bed < annotations.gtf \
| bedmap --count --echo --echo-map - \
| awk -v FS="\t" -v OFS="\t" '($1 > 1)' \
| cut -f2- \
> answer.bed
References:
• 0 views
•
link
written
by
Alex Reynolds
Log in to answer this question.
More posts like this
-
Extracting human exome from GTF file
written by Eugenio •
-
Genome annotation: convert gtf to gff3 file
written by shinyjj •
-
How robust is rMATS when swapping samples?
written by Apex92 •
-
What is a suitable Gene Ontology Enrichment Analaysis tool?
written by Apex92 •
-
RefSeq obtained from genome.ucsc.edu
written by Apex92 •
-
Which GTF file to use for FeatureCounts and STAR: low percent of assigned reads with primary assemb…
written by science_lizard •
-
Simple genome annotation using transcripts data
written by liorglic
-
Recurrent CNV analysis
written by enho •
-
Do you guys know how to use net file format?
written by thejustpark
-
List Of All Transcription Factors In Human
written by dfernan
You need to start with BedTools suite. It has many features to deal with coordinates, e.g cluster, merge, etc and the documentation is really easy to follow.
This is definitely helpful - thank you
There are some tools to do interval overlap, bedtools and R/Bioconductor (rtracklayer+GenomicRanges) comes to mind. If you want to do it on the commandline use the former, or inside R use Bioconductor.