[Bioperl-l] (no subject)
Hermann Norpois
hnorpois at googlemail.com
Tue Jun 5 21:10:13 UTC 2012
Hello,
actually I wish to get ensembl transcript IDs for some gene names. For
addressing this I tried some syntax ( to become more familiar with bioperl
and ensembl) but it doesnt work even for easier things, for instance for
getting DBEntries. So I guess there is a principal mistake in my script.
Can anybody help me?
Thanks,
Hermann Norpois
#!/bin/perl -w
use strict;
use Bio::EnsEMBL::Registry;
my $registry = 'Bio::EnsEMBL::Registry';
$registry->load_registry_from_db(
-host => 'ensembldb.ensembl.org', # alternatively '
useastdb.ensembl.org'
-user => 'anonymous'
);
#my $name = "Bdnf";
my $gene_adaptor = $registry -> get_adaptor( 'Mouse', 'Core', 'Gene');
my $gene = $gene_adaptor->fetch_by_stable_id('ENSG00000184129');
#my $gene = $gene_adaptor -> fetch_all_by_external_name ($name);
my @dbentries = @{ $gene->get_all_DBEntries() };
my $info = join (":", @dbentries);
print "$info\n";
# print $gene -> seq(), "\n";
More information about the Bioperl-l
mailing list