Permalink
Browse files

test script for larger systems

svn path=/bioperl-corba-client/trunk/; revision=58
  • Loading branch information...
1 parent 593b96d commit c17e7a909acad03e7433cf37a1334f0aabf273fa birney committed Feb 28, 2002
Showing with 17 additions and 4 deletions.
  1. +17 −4 scripts/test_seqdb.pl
View
21 scripts/test_seqdb.pl
@@ -17,13 +17,26 @@
print "ref is $ref\n";
my $db = Bio::CorbaClient::SeqDB->new( '-corbaref' => $ref);
-my $seq = $db->get_Seq_by_id('HUMBDNF');
#print "sequence is ",$seq->id,"\n";
-print "sequence is ",$seq->seq,"\n";
+open(F,shift) || die "Could not open file!";
+
+while( <F> ) {
+ ($acc) = split;
+ push(@accs,$acc);
+}
+
+$seqio = Bio::SeqIO->new( '-format' => 'fasta', -fh => \*STDOUT);
+
+foreach $acc ( @accs ) {
+ my $seq = $db->get_Seq_by_acc($acc);
+ #foreach $sf ( $seq->top_SeqFeatures ) {
+ #print STDERR "Got $sf with ".$sf->corbaref()."\n";
+ #print STDERR "Location ",$sf->corbaref()->location_as_string," vs ",$sf->location->to_FTstring,"\n";
+ #}
+ $seqio->write_seq($seq);
+}
-$seqio = Bio::SeqIO->new( '-format' => 'embl', -fh => \*STDOUT);
-$seqio->write_seq($seq);

0 comments on commit c17e7a9

Please sign in to comment.