[ensembl-dev] 5 prime UTR coordinates
Nick Fankhauser
lists at nyk.ch
Tue Mar 13 11:12:31 GMT 2012
Hi Andy
thanks for your helpful reply!
But when I try this example script, it only works on the first two
transcripts:
5' UTR 248790491-248790430 | 3' UTR 248789478-248789420
5' UTR 249120832-249119135 | 3' UTR 249106098-249104650
5' UTR 0-0 | 3' UTR 0-0
For the last one, it doesn't seem to work. Is this method maybe not
compatible with the current version of the API?
Yours,
Nick
On 09/03/12 16:55, Andy Yates wrote:
> Hi Nick,
>
> The method five_prime_utr() returns a Bio::Seq object & is meant to be a way of accessing the sequence data not the coordinates. Take a look at this script hosted at GitHub's gist:
>
> https://gist.github.com/1237394
>
> It should give you all the information you need to extract the coordinates you need
>
> All the best,
>
> Andy
>
> Andrew Yates Ensembl Core Software Project Leader
> EMBL-EBI Tel: +44-(0)1223-492538
> Wellcome Trust Genome Campus Fax: +44-(0)1223-494468
> Cambridge CB10 1SD, UK http://www.ensembl.org/
>
> On 9 Mar 2012, at 10:51, Nick Fankhauser wrote:
>
>> Hi!
>>
>> I'm trying to get the coordinates of the 5' UTRs of all transcripts of a
>> given gene. But the "start" method can't be located for the UTR object.
>> Is there a different approach that makes it possible to retrieve the
>> UTR's coordinates?
>>
>> This is what I tried:
>>
>> foreach my $transcript ( @{ $gene->get_all_Transcripts() } ) {
>> my $utr5 = $transcript->five_prime_utr();
>> next if not defined $utr5;
>> my $start5 = $utr5->start();
>> my $end5 = $utr5->end();
>> }
>>
>> Thanks!
>>
>>
>> Nick Fankhauser
>> ---
>> http://www.ccspmd.ethz.ch/people/fniklaus
>>
>> _______________________________________________
>> Dev mailing list Dev at ensembl.org
>> List admin (including subscribe/unsubscribe): http://lists.ensembl.org/mailman/listinfo/dev
>> Ensembl Blog: http://www.ensembl.info/
>
>
> _______________________________________________
> Dev mailing list Dev at ensembl.org
> List admin (including subscribe/unsubscribe): http://lists.ensembl.org/mailman/listinfo/dev
> Ensembl Blog: http://www.ensembl.info/
More information about the Dev
mailing list