Biostar Beta. Not for public use.
Homer annotatePeaks.pl custom GTF
0
Entering edit mode
20 months ago
Lucy • 0

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

ADD COMMENTlink
1
Entering edit mode

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

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.3.1