This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Homer annotatePeaks.pl custom GTF

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

homer atac-seq gtf

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.

0 answers

No answers yet.

Log in to answer this question.