[Bioperl-l] Out of memory errors running Bio::ASN1::EntrezGene against latest Homo_sapiens.ags file
Susan Wilson
smwilson at hpc.unm.edu
Fri Oct 12 13:45:41 EDT 2007
Hi,
I downloaded the latest ftp://ftp.ncbi.nlm.nih.gov/gene/DATA/
ASN_BINARY/Mammalia/Homo_sapiens.ags.gz and ran gene2xml on it to
generate Homo_sapiens.xml which is 5821420628 bytes. I cannot parse
this file with Bio::ASN1::EntrezGene, even on a machine with 256GB of
memory. I get a simple "Out of memory" output even with the
following code:
#!/usr/bin/perl
use strict;
use Bio::ASN1::EntrezGene;
my $parser = Bio::ASN1::EntrezGene->new('file' =>
"Homo_sapiens.xml");
while(my $result = $parser->next_seq)
{
}
Thanks.
Susan
More information about the Bioperl-l
mailing list