Question: Chip-Seq peaks annotation
0
gravatar for tarek.mohamed
24 months ago by
tarek.mohamed250
tarek.mohamed250 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="org.Hs.eg.db")

Thanks Tarek

chipseeker chip-seq • 1.2k views
ADD COMMENTlink written 24 months ago by tarek.mohamed250
Please log in to add an answer.

Help
Access

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