[ensembl-dev] Havana Gene in Compara
Nick Fankhauser
lists at nyk.ch
Tue May 10 16:13:20 BST 2011
I'm trying to get a member adaptor for Havana Genes (ncRNAs), for
example OTTHUMG00000002858.
The following script works for Ensembl Genes, but how do I have to
formulate it for Havana Genes?
#!/usr/bin/perl
use lib "/opt/ensembl/modules";
use lib "/opt/ensembl-compara/modules";
use strict;
use Bio::EnsEMBL::Registry;
Bio::EnsEMBL::Registry->load_registry_from_db(-host =>
'ensembldb.ensembl.org', -user => 'anonymous');
my $member_adaptor = Bio::EnsEMBL::Registry->get_adaptor('Multi',
'compara', 'Member');
throw("Cannot connect to Compara") if (!$member_adaptor);
# works:
my $member =
$member_adaptor->fetch_by_source_stable_id('ENSEMBLGENE','ENSG00000004059');
die("Gene not found!") if (!$member);
print $member->description,"\n";
# does not work:
my $member =
$member_adaptor->fetch_by_source_stable_id('HAVANAGENE','OTTHUMG00000002858');
die("Gene not found") if (!$member);
print $member->description,"\n";
The program dies on the second "fetch_by_source_stable_id" for
OTTHUMG00000002858. How can I fetch an adaptor for it?
More information about the Dev
mailing list