<div dir="ltr">Any chance you might consider implementing this?<div>I think it would be very useful to be able to retrieve the underlying sequence in previous assemblies since the transform method can already give us the coordinates of any given feature on a older assembly.</div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><font style="background-color:rgb(255,255,255)" color="#999999">=========================<br>     Duarte Miguel Paulo Molha      <br></font><div><font style="background-color:rgb(255,255,255)" color="#999999">         <a href="http://about.me/duarte" target="_blank">http://about.me/duarte</a>         <br>=========================</font></div></div>
<br><div class="gmail_quote">On Tue, Oct 7, 2014 at 8:53 AM, Andy Yates <span dir="ltr"><<a href="mailto:ayates@ebi.ac.uk" target="_blank">ayates@ebi.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Duarte,<br>
<br>
You can access GRCh37 from <a href="http://ensembldb.ensembl.org" target="_blank">ensembldb.ensembl.org</a> port number 3337. Ensembl databases currently only hold the contigs and assembly for a single assembly. That's why when you try to get sequence for GRCh37 in the GRCh38 database you get N's back<br>
<br>
Hope this helps,<br>
<br>
Andy<br>
<br>
------------<br>
Andrew Yates - Ensembl Support Coordinator<br>
European Molecular Biology Laboratory<br>
European Bioinformatics Institute<br>
Wellcome Trust Genome Campus<br>
Hinxton, Cambridge<br>
CB10 1SD, United Kingdom<br>
Tel: <a href="tel:%2B44-%280%291223-492538" value="+441223492538">+44-(0)1223-492538</a><br>
Fax: <a href="tel:%2B44-%280%291223-494468" value="+441223494468">+44-(0)1223-494468</a><br>
Skype: andrewyatz<br>
<a href="http://www.ensembl.org/" target="_blank">http://www.ensembl.org/</a><br>
<div><div class="h5"><br>
On 6 Oct 2014, at 17:23, Duarte Molha <<a href="mailto:duartemolha@gmail.com">duartemolha@gmail.com</a>> wrote:<br>
<br>
> Dear developers<br>
><br>
> I have the latest API downloaded and I would want to create a scritp that could retrieve sequence information from a specified assembly<br>
><br>
> so I have made a script that tries to accomplish this:<br>
><br>
> so ...assuming this coordinates:<br>
><br>
> chr: 5<br>
> from: 112043202<br>
> to: 112046226<br>
> strand = 1<br>
> assembly = GRCh38<br>
><br>
>     my $slice = $slice_adaptor->fetch_by_region( 'chromosome', $chrom, $from, $to, $strand, $assembly );<br>
><br>
><br>
>         $seq        = $slice->seq();<br>
><br>
> I can retrieve the dna sequence:<br>
><br>
> >CHR5-112043202-112046226       chr5:112043202-112046226<br>
> AGTATATAATCACAT..............CTAAAAGCAAACA<br>
><br>
><br>
> However, if I give it the variables:<br>
><br>
> chr: 5<br>
> from: 112043202<br>
> to: 112046226<br>
> strand = 1<br>
> assembly = GRCh37 or NCBI36<br>
><br>
> i get :<br>
><br>
> >CHR5-112043202-112046226       chr5:112043202-112046226<br>
> NNNNNNNNN.........NNNNNNNNNNNNNNNN<br>
><br>
><br>
> How can I get the correct underlying sequence?<br>
><br>
> Best regards<br>
><br>
> Duarte<br>
</div></div>> _______________________________________________<br>
> Dev mailing list    <a href="mailto:Dev@ensembl.org">Dev@ensembl.org</a><br>
> Posting guidelines and subscribe/unsubscribe info: <a href="http://lists.ensembl.org/mailman/listinfo/dev" target="_blank">http://lists.ensembl.org/mailman/listinfo/dev</a><br>
> Ensembl Blog: <a href="http://www.ensembl.info/" target="_blank">http://www.ensembl.info/</a><br>
<br>
<br>
_______________________________________________<br>
Dev mailing list    <a href="mailto:Dev@ensembl.org">Dev@ensembl.org</a><br>
Posting guidelines and subscribe/unsubscribe info: <a href="http://lists.ensembl.org/mailman/listinfo/dev" target="_blank">http://lists.ensembl.org/mailman/listinfo/dev</a><br>
Ensembl Blog: <a href="http://www.ensembl.info/" target="_blank">http://www.ensembl.info/</a><br>
</blockquote></div><br></div>