I'd like to perform a motif analysis of all promoter regions of the gene ENSMUSG00000020538 in mus musculus. To do so, I wanted to use Biomart :
library('biomaRt')
mart <- useMart("ensembl")
mart <- useDataset("mmusculus_gene_ensembl", mart = mart)
all.genes <- getSequence(id = 'ENSMUSG00000020538',
type="ensembl_gene_id",
seqType="coding_gene_flank",
upstream=100,
mart=mart)
But this provide me only the flanking region of the sepcified gene, not all its promoter regions. How could I do ?