Annotation bug using HOMER

Dear Users,

I am getting the below mentioned issue while executing annotatePeaks.pl using HOMER.

readline() on closed filehandle IN at /data/ngs/projects/Abhi/HOMER/bin/parseGTF.pl line 77.
readline() on closed filehandle IN at /data/ngs/projects/Abhi/HOMER/bin/parseGTF.pl line 77.
    Prioritizing Annotations:
    Annotating:
        Annotation    Number of peaks    Total size (bp)    Log2 Enrichment
    Annotating:
        Annotation    Number of peaks    Total size (bp)    Log2 Enrichment
    Counting Tags in Peaks from each directory...
    Organism: unknown
    Outputing Annotation File...
PeakID (cmd=annotatePeaks.pl peaks.txt reference.fa -gtf reference.gtf)    Chr    Start    End    Strand    Peak Score    Focus Ratio/Region Size    Annotation    Detailed Annotation    Distance to TSS    Nearest PromoterID    Entrez ID    Nearest Unigene    Nearest Refseq    Nearest Ensembl    Gene Name    Gene Alias    Gene Description    Gene Type

All the annotation columns generates NA values.

Can anyone help in rectifying this bug.

Thanks!

chip-seq

• 3,587 views • link

updated 4.0 years ago by Ram

• written 10.6 years ago by pinky_pinkpinky •

What file format are you using as the input to annotatePeaks.pl? Also what is the command used that generates the error?

• 0 views • link

written 10.6 years ago by Sam

0 answers

No answers yet.

More posts like this

Exon Annotation failed using homer annotatePeaks.pl
written 2.2 years ago by Alessia •

Hi, I am analyzing RNA-seq data for alternative splicing analysis and used for this reason RefSeq human genome as reference (hg38.p14). I have a list …
How to filter out the promoter-TSS peaks from the `.narrowPeak` file?
written 4.0 years ago by Dan •

I annotated the peaks using `annotatePeaks.pl macs2/${FILE}"_peaks.narrowPeak" mm10`, The annotation file is : ``` PeakID Chr Start End Strand Peak Score Focus Ratio/Region Size Annotation …
No annotation peaks when running homer
written 5.3 years ago by Lila M

Hi there, I'm trying to use [homer][1] to annotate my peaks, what I'm doing is perl annotatePeaks.pl peaks_cntr1_peaks.narrowPeak hg38 > annotated_peaks_ctr1 and I've got this …
annotationPeaks.pl algorithm for defining promoter-Tss
written 6.0 years ago by yahelsal •

Hi, I'm using annotationPeaks.pl for extracting peaks annotations from some bed file. I assume that Homer defines a peak to be 'promoter-TSS' if the peak …
HOMER error -bash: ./ctcf.motif: Permission denied
written 6.3 years ago by gerardo_za_94 •

Hello I'm running HOMER on windows through cygwin64. installation was successful, Can run .pl commands. regardless, it seems that the commands cannot access the motif …
Empty output when use HOMER annotatePeak.pl
written 6.7 years ago by deb0612 •

when I use the following command USAGE: perl bin/annotatePeaks.pl 9A26_MarkDuplicates_summits.bed hg38 I got massage below ===================================================== Peak file = RCC/ATAC-seq/9A26_MarkDuplicates_summits.bed Genome = hg38 Organism = …
Peak annotation with HOMER. Statistics with the annStats parameter
written 7.0 years ago by nanoide

Hi all. So I've been annotating some peaks with HOMER annotatePeaks. I have followed the instructions in the [webpage][1] for changing the TSS definition: *Although …
homer-annotation of the minus strand genes
written 8.6 years ago by c2c12 •

Hi, I'm trying to use Homer to annotate peaks from the analysis in the Diffbind . I noticed recently that even when I'm assigning peaks …
Annotating significant interactions of Hi-C data of Barley (Hordeum vulgare)
written 9.5 years ago by s.singh •

Hi, I am using HOMER to analyse the Hi-C data of barley. I have found all the significant interactions. I am facing problem in annotating …
Annotate Chip-Peaks With Gene Symbol/Name
written 13.3 years ago by GeneInfo •

<p>I want to use annotatePeaks.pl from HOMER to annotate my ChIP-seq peaks. It works but the column for <code>gene name.etc</code> is missing.</p> <p>Here is the …