Entering edit mode
10.4 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!
Is this true when you use biomart web instead biomaRt?
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.