<div dir="ltr">Thanks Kieron<div><br></div><div>I understand your point of view... but I still think there is a case for a couple of methods to be implemented in the transcript object:</div><div>@{$transcript->get_coding_regions} and @{$transcript->get_noncoding_regions} </div><div><br></div><div>Both returning feature objects. Am I the only one to find these useful? I hope not :)</div><div><br></div><div>Thanks</div><div><br></div><div>Duarte</div><div><br></div><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><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></div>
<br><div class="gmail_quote">On 16 April 2015 at 16:19, Kieron Taylor <span dir="ltr"><<a href="mailto:ktaylor@ebi.ac.uk" target="_blank">ktaylor@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>
The coordinates you’re getting back are pre-splicing. The method you’re calling is from the Transcript class, hence the response is with reference to that object. If you’re after exon coordinates, you should be attempting to work with exon objects, such as fetching the exons of the transcript and asking them for coding_region_start($transcript) until numbers start appearing. Your workaround is also a valid approach.<br>
<br>
My explanation isn’t very satisfactory, but we try to avoid writing methods that need complex return types, such as the list of lists required for your usecase. More often than not, users require other attributes of the objects too, so you would still end up with a list of exons. I hope that helps.<br>
<br>
Regards,<br>
<br>
Kieron<br>
<br>
<br>
Kieron Taylor PhD.<br>
Ensembl Core senior software developer<br>
<br>
EMBL, European Bioinformatics Institute<br>
<div class="HOEnZb"><div class="h5"><br>
<br>
<br>
<br>
<br>
> On 16 Apr 2015, at 09:15, Duarte Molha <<a href="mailto:duartemolha@gmail.com">duartemolha@gmail.com</a>> wrote:<br>
><br>
> Anyone able to provide me some help on this?<br>
><br>
> I have now found away around this issue by finding the exonic regions within the reported URT, but would very much like to understand the thinking behind this.<br>
><br>
> Best regards<br>
><br>
> Duarte<br>
><br>
><br>
> =========================<br>
>      Duarte Miguel Paulo Molha<br>
>          <a href="http://about.me/duarte" target="_blank">http://about.me/duarte</a><br>
> =========================<br>
><br>
> On 14 April 2015 at 13:40, Duarte Molha <<a href="mailto:duartemolha@gmail.com">duartemolha@gmail.com</a>> wrote:<br>
> Dear Developers<br>
><br>
> Please consider the transcript :<br>
><br>
> ENST00000470357<br>
><br>
> I am trying to retrieve the coordinates of UTR regions of this transcript<br>
> To this end I have a script that basicaly starts with the transcript feature object $transcript<br>
><br>
> my $five_prime  = $transcript->five_prime_utr_Feature;<br>
><br>
> $feature_params->{start} = $five_prime->start;<br>
> $feature_params->{end}  = $five_prime->end;<br>
><br>
> However, in this case the script will output the coordinates from the start of the 1st non_coding exon to the end of the non-coding portion of the 3rd exon (chr1     7772707 7777171).<br>
> How can I change this so that the script will only output the coordinates of the non-coding exon portions?<br>
><br>
> In this case I would like to output:<br>
><br>
> chr1  7772707 7773198<br>
> chr1  7773442 7773511<br>
> chr1  7777160 7777171<br>
><br>
> This there a simple way of achieving this?<br>
><br>
> Many thanks<br>
><br>
> Duarte<br>
><br>
><br>
</div></div><div class="HOEnZb"><div class="h5">> _______________________________________________<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>
</div></div></blockquote></div><br></div>