[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