[ensembl-dev] Core API quits with no error messages
ian Longden
ianl at ebi.ac.uk
Fri Aug 26 15:16:03 BST 2011
What version of the API are you using?
ENSG00000256461 has been archived and is not in release 63 so that
could be part of it.
also add -verbose => 1 to the load_registry_form_db to see which
databases it is connecting to.
So i modified your script to use our server as a test
#######################################################
use Bio::EnsEMBL::Registry;
my $reg = 'Bio::EnsEMBL::Registry';
$reg->load_registry_from_db(
-host => 'ensembldb.ensembl.org',
-user => 'anonymous',
-db_version => 63,
-verbose => 1
);
my $GA = $reg->get_adaptor('Human', 'Core', 'Gene');
my $gene = $GA->fetch_by_stable_id('ENSG00000256461');
if(defined($gene)){
my @exons = @{ $gene->get_all_Exons() };
foreach my $e (@exons) {
print ">" . $e->stable_id() . "\n";
print $e->seq->seq . "\n\n";
}
}
else{
print "Gene could not be found\n";
}
exit;
##################################################################
which gave me the list of databases and also the message:-
Gene could not be found
-Ian Longden
EnsEMBL Developer.
On Fri, Aug 26, 2011 at 2:36 PM, Damian <dfermin at umich.edu> wrote:
> Hello.
>
>
> I don't know how to debug this problem with the Core API in perl.
>
> Here is the code:
> ####################################################
> use Bio::EnsEMBL::Registry;
>
> my $reg = 'Bio::EnsEMBL::Registry';
> $reg->load_registry_from_db(
> -host => '127.0.0.1',
> -user => 'nesviLab'
> );
>
>
> my $GA = $reg->get_adaptor('Human', 'Core', 'Gene');
>
> my $gene = $GA->fetch_by_stable_id('ENSG00000256461');
>
> my @exons = @{ $gene->get_all_Exons() };
>
> foreach my $e (@exons) {
> print ">" . $e->stable_id() . "\n";
> print $e->seq->seq . "\n\n";
> }
>
> exit;
>
> ####################################################
>
>
> When I run this code, the output is the first exon of the array and that's
> it:
>>ENSE00002226834
>
> The sequence of the exon is never written out.
> The script does not throw an error and the program does exit.
> Running 'top' on my computer shows that the script is not even using up CPU
> cycles.
>
> This only happens with genes on chromosome 1. The above code runs to
> completion with any genes from chromosomes 2..Y
> Any idea what's going on? Any and all help is greatly appreciated.
>
> I'm running the Ensembl Core version 62 on RHEL 5.4 (64bit).
> Perl version: 5.8.8
> MySQL version: 5.0.77
>
> Damian
>
> --
> ------------------------------------------------------------------------------
> Damian Fermin, Ph.D.
> Pathology Department
> University of Michigan
> 4237 Medical Science I
> 1301 Catherine
> Ann Arbor, MI 48109-0602
> 734.615.0302
>
> "There is No Gene for the Human Spirit"
> -- GATTACA
> ------------------------------------------------------------------------------
>
>
> _______________________________________________
> Dev mailing list Dev at ensembl.org
> List admin (including subscribe/unsubscribe):
> http://lists.ensembl.org/mailman/listinfo/dev
> Ensembl Blog: http://www.ensembl.info/
>
More information about the Dev
mailing list