[ensembl-dev] Difficulty finding genomic slices using Ensembl Perl API.

Allan Kamau kamauallan at gmail.com
Wed Mar 9 08:04:12 GMT 2011


I am trying to obtain sections of sequences (and perhaps accompanying
annotation) of several organisms, for each such request I have the
following pieces of information.
1)species name
2)genome build name
3)chromosome, scaffold or contig name
4)start and end positions
5)orientation of the sequence (fwd or rev)

When I query for the above for some species such as
"anopheles_gambiae" I get no results found in the object
"@db_adaptors" (see below).


use Bio::EnsEMBL::Registry;

my $registry='Bio::EnsEMBL::Registry';

$registry->load_registry_from_db(
       -host=>'ensembldb.ensembl.org',
       -user=>'anonymous'
);
my @db_adaptors=@{Bio::EnsEMBL::Registry->get_all_DBAdaptors(-species=>"anopheles_gambiae")};

What could I be doing wrong.

The code

$species_name2='Anopheles gambiae';
$species_name_available=$registry->get_alias($species_name2);

returns
"anopheles_gambiae"

Meaning that "anopheles_gambiae" is known to ensembl as a species name
but the outcome of running the preceding code indicates that there are
no db_adaptors for this species.

Is there a graceful way of obtaining the information I am looking for
(which I pointed out in the beginning)?

Allan.




More information about the Dev mailing list