ChIPSeeker Help for Annotating Peaks

Hello,

I am trying to annotate peaks obtained from ATAC-seq experiments using ChIPSeeker. (Doing this on usegalaxy.org)

I am getting an error but not sure why. I have attached screenshots of the error and my input files (narrow peak file, GTF from UCSC).

error

narrowpeak file

UCSC GTF for mm9

Any help or thoughts would be appreciated, thanks!

chip-seq atac-seq software error

• 2,297 views • link

updated 5.3 years ago by gglim •

• written 5.6 years ago by saifdab99 •

1 answer

The error message

subscript contains invalid names

I guess...there's something wrong with your peakname i.e. 4th column of narrowPeak file.

• 0 views • link

written 5.3 years ago by gglim •

More posts like this

Fixed peak width in ENCODE TFs ChIP-seq
written 1.8 years ago by DavideBrex •

Dear all, I have a question about the TFs ChIP-seq peak width of the ENCODE project. It is probably something I am missing. If I …
Overestimate of ATAC-seq differential peaks
written 1.9 years ago by Darryl •

Hello I am comparing 2 conditions of ATAC-seq. Peaks were called using MACS2, broad setting. The number of peaks in condition A is larger than …
MACS2: broad and broad no-model output peaks are exactly the same. Is something wrong?
written 4.1 years ago by piyushjo

Hi, I am trying peak calling for H3K27ac ChIP-seq for enhancer and super enhancer analysis. Initially, I was planning to use default narrowPeak mode in …
How to perform custom annotation using garfield_annotate_uk10k.sh script in GARFIELD tool?
written 5.4 years ago by bk11 •

Hi, I am trying to use the GARFIELD V2 tool ([https://www.ebi.ac.uk/birney-srv/GARFIELD/][1]) and create custom annotation data using narrow peak bed files from Roadmap [https://egg2.wustl.edu/roadmap/data/byFileType/peaks/consolidatedImputed/narrowPeak/][2]. GARFIELD …
Genrich analyse ATAC-seq data problem
written 5.8 years ago by MatthewP

I am trying to analyse ATAC-seq data with [Genrich](https://github.com/jsh58/Genrich) with 2 replicates. I want to ask a few questions. 1. Output file(narrowPeak) by `-o` is …
Invalid contig error when running CLIPper with dm3
written 6.0 years ago by jkkbuddika

Hello, I am trying to use CLIPper (https://github.com/YeoLab/clipper) to identify peaks on eCLIP data generated from Drosophila tissues. When I am running the program using …
Roadmap epigenomics narrowPeak and broadPeak contain scores > 1000
written 7.4 years ago by F Lekschas •

How should the score of consolidated narrowPeak and broadPeak Roadmap Epigenomics files be interpreted? The bed format dictates that the score has to be between …
kmer distribution tiny peak at low coverage
written 8.4 years ago by deepti1rao •

[jellyfish.histo file][1]<br> [fastqc file][2] I generated a kmer count file using jellyfish and subsequently a histogram, which when plotted in R gave the attached graph. …
Read counts based on genes from samtools, feature count and htseq-count
written 10.2 years ago by bioinforesearchquestions

Hi All, I am currently trying to find number of reads (sequenced from DNAseq) that are mapped to human genes. So I tried three different …
MACS2 error in Galaxy
written 10.6 years ago by haithem.barbour •

Hi, I was trying to call peak using MACS2_callpeak on a ChIP-Seq data set. The alignment went fine but the peak calling failed even though …