[Bioperl-l] Blasting multiple protein sequences.

Chet Seligman cseligman at earthlink.net
Mon Apr 26 20:56:01 EDT 2010


I am a student and new to BioPerl.
I have a fasta formatted file with 17 protein sequences in it and I would
like to blast them all.
My problem is that only the first one gets blasted.
Can anyone make a suggestion?
Thanks in advance - Chet
Here's my code:
#!/usr/bin/perl
use Bio::Perl;
use strict;

my $fastafilename = 'protein.fa';
# reads an array of sequences
my @seq_object_array = read_all_sequences($fastafilename,'fasta');
 
# BLAST a sequence (assumes an internet connection)

my  $blast_report = blast_sequence(@seq_object_array);

write_blast(">blast.out",$blast_report);
  
exit;  


More information about the Bioperl-l mailing list