I am currently trying to perform part of my ChIPseq analysis with the ChIPseeker program.
I would like to use the ''annotatePeak'' function but for that apart for the bed file with my peaks (derived from MACS), which I have it, I need two more files, the ''MeSH.Miy.eg.db'', which I also have it (from https://www.bioconductor.org/packages/devel/data/annotation/), and a ``TxDb.Lotus japonicus.knownGene`` file that I cannot find it somewhere.
Could anyone help with that?
This is an example of the command that I should follow (from the manual of the ChIPseeker):
>peakAnno <- annotatePeak(files[], tssRegion=c(-3000, 3000), TxDb=txdb, annoDb="org.Hs.eg.db")
# txdb <- TxDb.Hsapiens.UCSC.hg19.knownGene
I would like to do the same but with Lotus japonicus.