[Bioperl-l] standaloneblast large seq retrieving

[Guojun Yang]

Hi,all,
I was trying to use the following sub to get seq after a standaloneblast. It worked with DB with short entries (~200kb), but it failed to work with a DB with much longer entries (up to ~30 Mb an entry). Can anybody give me a hint?

 sub getseq {
my $name=$_[0];
my $file_name = $_[1];
my $inx=Bio::Index::Fasta->new (-filename => $file_name.".idx",
                                -write_flag => 1);
$inx->id_parser(\&get_id);
$inx->make_index($file_name);
$seq = $inx->fetch($name);  
return $seq;
	    }

Thanks,
Yang