use strict; use LWP::Simple; my @fasta; my $o; # Note: @accession defined earlier has 603 elements for (my $i=0; $i<@accession; $i++) { if ($accession[$i] =~ /^(\w+)\_(\w+)\_(\w+)\;/) { #print "$accession[$i]\n"; if ($1 eq $3) { # USE $1 AS THE ACCESSION NUMBER my $seq = get "http://us.expasy.org/cgi-bin/get-sprot-fasta?$1"; #print "$seq\n"; push (@fasta, $seq); push @fasta, "> COULDN'T FIND IT" unless defined $seq; } else { # # USE $3 AS THE ACCESSION NUMBER my $seq = "http://us.expasy.org/cgi-bin/get-sprot-fasta?$3"; push (@fasta, $seq); push @fasta, "> COULDN'T FIND IT" unless defined $seq; } } } #print "$o\n"; # fasta only has 537 / 603 elements print "@fasta\n";