[BioC] BiomaRt and extracting genomic sequence

Rob Illingworth rillingw at staffmail.ed.ac.uk
Wed Jan 28 09:08:29 CET 2009


Hi Steffen

Thank you very much for your help that works perfectly.

I will keep checking to see if and when ensembl make the current build 
available for this.

kind regards

Rob Illingworth


steffen at stat.Berkeley.EDU wrote:
> Hi Rob,
>
> This is currently not possible using biomaRt in web service mode.  However
> the Ensembl team will probably make this available to us in the near
> future.
>
> If you're using the current release version of biomaRt (1.16.0) then you
> should be able to get genomic sequences (unfortunately only from Ensembl
> version 47) by:
>
>   
>> mart = useMart("ensembl_mart_47", dataset="hsapiens_gene_ensembl",
>>     
> mysql=TRUE)
>   
>> seq = getSequence(chromosome = 20, start=20000, end = 20010, mart=mart)
>> seq
>>     
>   chromosome start   end    sequence
> 1         20 20000 20010 TGAGGTCTTAC
>
>
> Note that this MySQL mode retrieval of genomic sequences will be replaced
> by the default web service mode use of biomaRt.
>
> Cheers,
> Steffen
>
>   
>> Hello to the list
>>
>> I have recently started mapping next generation sequencing data to the
>> human genome and would now like to map this to genes and other
>> annotated genomic features.  Whilst I have found that I can map and
>> sort genes with the biomaRt package I have completely failed to
>> extract genome sequence using simple chromosome coordinates. I have
>> looked through the documentation for this package and have tried to
>> use the getSequence function but to no avail. The 'seqType' values all
>> seem to relate to gene associated sequences and 'type' has no obvious
>> value for genomic sequence and as such return an error msg.  I think
>> there is probably a very straight forward solution which I am missing
>> and  I would be very grateful if anyone could help.
>>
>> Thanks for your assistance
>>
>> kind regards
>>
>> Rob Illingworth
>>
>> --
>> WTCCB
>> Institute of Cell Biology
>> The University of Edinburgh
>> Michel Swann Building
>> Mayfield Road
>> Edinburgh
>> EH9 3JR
>> TEL:+44(0)131 650 8695
>>
>>
>>
>>
>> --
>> The University of Edinburgh is a charitable body, registered in
>> Scotland, with registration number SC005336.
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at stat.math.ethz.ch
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>>
>>     

>
>
>
>   

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/bioconductor/attachments/20090128/ce52395c/attachment.pl>


More information about the Bioconductor mailing list