I was dealing with this package called 'Chipseeker' . I tried annotating a chip-seq peak calls fro Mycobacterium tuberculosis. As there is no 'org.Mtu.eg.db' package for annotation & 'TxDB' can be used in Chipseeker, so I created TxDB package & annotated as follows :
library(GenomicFeatures) txdb <- makeTxDbFromGFF("./m.tub1.gff3",circ_seqs=character(), organism="Mycobacterium tuberculosis") peakAnno <- annotatePeak(peak1, tssRegion=c(-3000, 3000), TxDb=txdb, genomicAnnotationPriority = c("Promoter", "5UTR", "3UTR", "Exon", "Intron", "Downstream", "Intergenic"))
The chip-seq was done for 'Gyrase' but annotations from Chipseeker gave mostly 'promoters' . Is is correct or am I suppose to add some more steps? Please help me.
Thanking you in advance, Anupriya