use Bio::DB::Fasta; my $fastaFile = shift; my $queryFile = shift; my $db = Bio::DB::Fasta->new( $fastaFile ); open (IN, $queryFile); while (){ chomp; $seq = $_; my $seqstr = $db->get_Seq_by_id($seq)->seq; if (!defined( $seqstr )) { die "Sequence $seq not found. \n" } print ">$seq\n", "$seqstr\n"; }