Obtaining promoter sequences with biomaRt R package
0
0
Entering edit mode
9.8 years ago
ilyco ▴ 60

Hello,

I tried to obtain the promoter sequences of 22 genes (2000 bp upstream) using the biomaRt R package and their Ensembl Gene Ids. However, the getBM() function returns 'Sequence Unavailable' for 20 of the genes and only retrieves the sequences for 2 of the genes. Do you know why this is happening and what alternatives do I have for obtaining the sequences?

I am using an archived version of the Ensembl mart for mm9 mouse genome and the following code:

promoters <- getSequence(id = gene.names, type='ensembl_gene_id',seqType='coding_gene_flank',
                upstream=2000,mart=ensembl67)

where gene.names is a list of ensembl gene ids.

Thank you very much in advance!

sequence biomaRt promoter • 5.1k views
ADD COMMENT
0
Entering edit mode

Is this true when you use biomart web instead biomaRt?

ADD REPLY
0
Entering edit mode

Thank you!

I downloaded all the upstream sequences of all the genes from biomaRt and selected the genes I wanted afterwards using the Ensembl Ids.

ADD REPLY

Login before adding your answer.

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