Homer annotatePeaks.pl output
0
0
Entering edit mode
25 days ago
Rakesh • 0

Hi All,

I am using homer annotatePeaks.pl (version 4.11) by providing the bed file. Surprisingly, in the output file all peaks start position is increasing by +1, while the end remains same as in the input bed file. I read through homer annotatePeaks.pl help page but not finding any reason specific to this context.

"Command used":

annotatePeaks.pl  Tumor_vs_Normal_peak_annotation_input_file.bed  hg38  -annStats  Tumor_vs_Normal_peaks_stats.txt  -CpG   -go GO_analysis_results   -genomeOntology   genome_ontology_results  >Tumor_vs_Normal_homer_annotation_output.txt"

Input file: tail Tumor_vs_Normal_peak_annotation_input_file.bed

chrM    10857   11082   peak_677899

chrM    13955   14612   peak_677900

As we seen in output file chr start site is. 13955 and 10857.

"Output file (only first few columns are shown):" tail Tumor_vs_Normal_homer_annotation_output.txt

peak_677899 chrM    10858   11082   +   0    NA      Intergenic ncRNA   0.013393    0.417778

peak_677900 chrM    13956   14612   +   0    NA  Intergenic ncRNA   0.013720    0.429224

As we seen in output file chr start site is. 13956 and 10858.

I do understand that annotation tools can be 0 and 1-based. But is there a way to retain peaks regions similar to input bed file? or I need to subtract 1 from all entries in that column?

Thanks

Rakesh

start annotatePeaks.pl chr homer position • 232 views
ADD COMMENT

Login before adding your answer.

Traffic: 6036 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6