Chip-Seq peaks annotation
0
0
Entering edit mode
7.0 years ago
tarek.mohamed ▴ 360

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

ChIP-Seq chipseeker • 2.9k views
ADD COMMENT

Login before adding your answer.

Traffic: 1957 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