[Bioperl-l] PARSE GENBANK SEQUENCE FILES
Jason Stajich
jason@cgt.mc.duke.edu
Mon, 10 Jun 2002 13:31:18 -0400 (EDT)
Whoops - didn't finish that:
> Now if you wanted to only print out sequences that match a criteria, like
> perhaps were of a certain species as you allude to below you could augment
> the while loop above like this
while( my $seq = $input->next_seq ){
next unless ( $seq->species->species =~ /rattus/i );
$output->write_seq($seq);
}
--
Jason Stajich
Duke University
jason at cgt.mc.duke.edu