Question: Homer annotatePeaks.pl custom GTF
0
gravatar for Lucy
16 months ago by
Lucy40
Lucy40 wrote:

Hi,

I am running Homer annotatePeaks.pl to assign my ATAC-seq peaks to genes. For the annotation, I am using a custom Ensembl GTF file and a custom Ensembl TSS file.

annotatePeaks.pl peaks.txt hg38 -gtf Ensembl_hg38.gtf -cTSS Ensembl_hg38.tss > annotated_peaks.txt

When using a custom GTF file, the Homer guidance states that "It will also use the GTF file's definition of TSS/TTS/exons/Introns for Basic Genome Annotation."

I was wondering what column of an Ensembl GTF file would be used for this? Does Homer work out which regions are intronic, intergenic etc. based on the regions that are assigned as exons or start_codons in field 3 of the GTF file, as I can't find any regions that are explicitly assigned as introns or intergenic regions?

Many thanks,

Lucy

atac-seq homer gtf • 993 views
ADD COMMENTlink modified 16 months ago • written 16 months ago by Lucy40
2

GTF never explicitely contain intergenic and intronic annotations. Intronic is the complement of all explicitely-annotated ranges of a gene within the gene border and intergenic is the complement of all gene ranges with the genome.

ADD REPLYlink written 16 months ago by ATpoint35k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 905 users visited in the last hour