[Bioperl-l] Swissprot query - Help!

SG Edwards s0460205 at sms.ed.ac.uk
Tue Mar 15 13:37:36 EST 2005


Hi, sorry for the obvious question but I'm really new to Perl/BioPerl!!

I want to run a script that sends a query to swissprot and returns the list of
sequences as Seq objects. I have tried the following code which throws an
exception 'MSG: Must speciy a value for uids to query'.

use Bio::DB::SwissProt;

$query = "Arabidopsis[ORGN] AND topoisomerase[TITL]";

$sp_obj = Bio::DB::SwissProt->new;

$stream_obj = $sp_obj->get_Stream_by_query($query);

while ($seq_obj = $stream_obj->next_seq) {
#print out the id
    print $seq_obj->display_id, "\n";
}

exit;


Any help is greatly appreciated!


More information about the Bioperl-l mailing list