This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Peaks annotation on bacterial genome

Hi all, I'm trying to annotate peaks identified in ChIP-seq experiments in M. tuberculosis using ChIPpeakAnno on R. I would like to distinguish between peaks that are located in promoter regions or in coding regions. Is it that possible? Is there any code I could use to selectively distinguish between these two regions? Can I use a custom genome that includes also promoters in addition to CDS?

Thank you in advance for your help.

chippeakanno annotation bacteria

2 answers

You can call annoPeaks with a GRanges annotation object that includes promoters, then use assignChromosomeRegion to summarize peak distributions associated with different genome elements. See the manual:

assignChromosomeRegion

Summarize peak distribution over exon, intron, enhancer, proximal promoter, 5 prime UTR and 3 prime UTR

Thank you for your kind reply! I have just a further question. Which is the difference between using annoPeaks and annotatePeaksinbatch with ChipPeakAnno?

Thank you for your help

Log in to answer this question.