[Bioperl-l] Needed help

Raman Deep Singh deep.raman at gmail.com
Tue Feb 7 20:16:48 UTC 2006


Hi all
     I have a huge task of retrieving a number of sequences from the
swiss prot databases on some fixed criteria. FOr that i want to index
the swiss prot database on my local disk. I have downloaded the whole
swiss prot database on my local disc  (the january 2006 release).

  I am currently using the bioperl on linux machine . I am using the
code listed below


=======================

    use Bio::Index::Swissprot;

    my $Index_File_Name = shift;
    my $inx = Bio::Index::Swissprot->new('-filename' => $Index_File_Name,
 '-write_flag' => 'WRITE');
    $inx->make_index(@ARGV);
-----------------------------------------
    # Print out several sequences present in the index
    # in gcg format
    use Bio::Index::Swissprot;
    use Bio::SeqIO;

    my $out = Bio::SeqIO->new( '-format' => 'gcg', '-fh' => \*STDOUT );
    my $Index_File_Name = shift;
    my $inx = Bio::Index::Swissprot->new('-filename' => $Index_File_Name);

    foreach my $id (@ARGV) {
        my $seq = $inx->fetch($id); # Returns Bio::Seq object
        $out->write_seq($seq);
    }

    # alternatively

    my $seq1 = $inx->get_Seq_by_id($id);
    my $seq2 = $inx->get_Seq_by_acc($acc);


-- -------------------------------
i am running teh script as

 perl getseqfromid.pl sample.dat

from the shell

and i am getting this error repeatedly

------------- EXCEPTION  -------------
MSG: Can't open 'DB_File' dbm file 'swiss100.dat' : No such file or directory
STACK Bio::Index::Abstract::open_dbm
/usr/lib/perl5/site_perl/5.8.5/Bio/Index/Abstract.pm:389
STACK Bio::Index::Abstract::new
/usr/lib/perl5/site_perl/5.8.5/Bio/Index/Abstract.pm:150
STACK Bio::Index::AbstractSeq::new
/usr/lib/perl5/site_perl/5.8.5/Bio/Index/AbstractSeq.pm:91
STACK toplevel i.pl:6


--------------------------
At some place online, i also found some document that some variables
need to be exported. I also did the same but still got teh same errors

kindly  help




Ramandeep Singh




More information about the Bioperl-l mailing list