lparsons / kas-seq-workflow

KAS-Seq Analysis Workflow

Geek Repo:Geek Repo

Github PK Tool:Github PK Tool

Annotated peak files have data mixed with header

lparsons opened this issue · comments

The ChIPSeeker annotatePeak function is not properly reading the broadpeak files which have no header. Thus, the first peaks data is getting mixed in with the headers. Adding a header line seems to address this before running annotatePeak.

chr>start>  end>name>   score>  strand> fold_enrichment>-log10(pvalue)> -log10(qvalue)