#!/usr/bin/perl

use strict;
use warnings;
use Bio::SeqIO;
use Data::Dumper;
my %seq_id;
my $fasta_id;
	    open HEADER , "<FASTA.headers" or die $!;
		while (<HEADER>){
			chomp $_;
			#print $_."\n";
			$fasta_id = $_; 
			$fasta_id =~ s/_.*//g ;
			#print $fasta_id."\n";				
			%seq_id = ("$fasta_id" => "$_");
#			print Dumper (\%seq_id);


            my $infile=$ARGV[0] || die ("Please give me an input fasta file\n");
            my $inseq = new Bio::SeqIO(-format => 'fasta',
                        -file => $infile);
            while (my $seq_obj = $inseq->next_seq ) {
 		    my $id = $seq_obj->id ;
			chomp $id;
#			print $id."\n";
			my $seq = $seq_obj->seq ;
				if (exists ($seq_id{$id})) {
				print ">";
				print $seq_id{$fasta_id};
				print "\n".$seq."\n";
				} 
				

			}
}