HOMER annotatePeaks.pl with GFF. How to specify a feature?
Entering edit mode
19 months ago
callapus • 0

HI, Having an issue running annotatePeaks.pl with a GFF annotation file (working in a non-standard organism and this is all I have available to me). Unfortunately in my GFF file exons are actually marked as CDS so annotatePeaks is not able to annotate my called peaks. The output:

Features that will be considered: exon

0 lines contained useful feature information This seem very low - consider specifying specific features (i.e. 'exon', 'transcript', etc.) or use "-keepAll"

My problem is I can't figure out how to specify the feature 'CDS' in my GFF. It does seem like it's possible to specify features however. Ive gone back and changed CDS to exon in my GFF file and annotatePeaks runs fine but would still find this really usefull if anyone knows how to do it! Thanks in advance!! C

HOMER ChIP-Seq peakcalling annotation GFF • 543 views
Entering edit mode
10 months ago
Qiongyi ▴ 120

You may be able to use the option of "-keepAll" or "-features" as I saw the source code of parseGTF.pl offer these options. However, exon is the essential feature that GTF/GFF file should contain, which is used to determine the TSS and TES, etc. It is better that you use a proper GTF/GFF file for the annotation. Simply replacing CDS to exon is not a correct way (CDS != exon).


Login before adding your answer.

Traffic: 2161 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6