[Bioperl-l] remote db retrieval
   
    Mikaela Ilinca Gabrielli
     
    MILG@lundbeck.com
       
    Wed, 20 Nov 2002 14:20:30 +0100
    
    
  
Can somebody help me?
I'd like to retrieve sequences from a remote database in FASTA format. I
cannot find how to do so allthough there should be a simple way.
Here's my code:
#! /usr/bin/perl
use Bio::SeqIO;
use Bio::DB::RefSeq;
my @acc=('NP_061843.1', 'NP_110411.1', 'NP_072093.1');
foreach my $ACC(@acc){
my $fasta= GetSeqDB($ACC);
print "$fasta\n";
}
sub GetSeqDB{
        my $acc = $_[0];
        my $db = new Bio::DB::RefSeq(-format=>'fasta');
        $db->proxy(['http', 'ftp'],'http://fw.lundbeck.com:8080');
        $seq = $db->get_Seq_by_acc($acc);
	  return $seq->seq();
        }
--Thanks!
/Mikaela