[ensembl-dev] accessing the tilepath entries programatically

Duarte Molha duartemolha at gmail.com
Tue Jun 30 11:50:10 BST 2015


I used this code to get all the gebnomic coordinates of your subcontigs:


my @slices = @{ $slice_adaptor->fetch_all('clone') };
foreach my $slice (@slices){
$progress->update();
my $clone_name =  $slice->seq_region_name();
my $projection = $slice->project('toplevel');
foreach my $segment ( @{$projection} ) {
my $to_slice = $segment->to_Slice();
print join "\t", ("chr".$to_slice->seq_region_name(), $to_slice->start(),
$to_slice->end(), $clone_name."\n");
}
}

However, by doing this, the database does not fetch the original clone name

for example.. using this script I get
chr4    47567235        47733411        AC092597.1

However I would like to get :

chr4    47567235        47733411        RP11-100N21

Can someone explain what I am doing wrong?

Thanks

Duarte



=========================
     Duarte Miguel Paulo Molha
         http://about.me/duarte
=========================

On 30 June 2015 at 09:45, Duarte Molha <duartemolha at gmail.com> wrote:

> Dear devs
>
> How can I search for a specific clone id present on your tilepath
>
> for example this: RP5-892C22
>
> I would like to use the perl API if possible
>
> Many thanks
>
> Duarte
>
>
>
> =========================
>      Duarte Miguel Paulo Molha
>          http://about.me/duarte
> =========================
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mail.ensembl.org/pipermail/dev_ensembl.org/attachments/20150630/914f444c/attachment.html>


More information about the Dev mailing list