Question: Chip-Seq peaks annotation
gravatar for tarek.mohamed
3.8 years ago by
tarek.mohamed270 wrote:

Hi All, I am analyzing a chip seq data trying to identify DNA binding sites to a specific protein.I did the alignment using bwa, and peaks calling using macs2. Now I am using R package "chipseeker" to annotate these peaks. My question is how can I excluded micro-RNA encoding genes from the annotation. I would guess that I have to remove them from the txdb I used for annotating the peaks. Is there an easy way to do this?

library(ChIPseeker) library("TxDb.Hsapiens.UCSC.hg38.knownGene") txdb<-TxDb.Hsapiens.UCSC.hg38.knownGene peakAnnotation <- annotatePeak("narrowpeak_bowtie_q0.05_IgG.txt", tssRegion=c(-3000, 3000), TxDb=txdb, annoDb="")

Thanks Tarek

chipseeker chip-seq • 2.1k views
ADD COMMENTlink written 3.8 years ago by tarek.mohamed270
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1529 users visited in the last hour