#!/usr/bin/perl -w
#usage: ./matchMF_HCtug.pl MF.TUG.fasta MF_matchGparental.txt > MF_matchGparental.fsa

my %SEQ;
$/ = '>';
open(FILE, "$ARGV[0]") || die("Cannot open");
while(<FILE>){
	next unless /\w/;
	my @lines = split(/\n/, $_);
	my ($ID) = $lines[0] =~ /(\S+)/;
	shift @lines;
	my $seq = join('', @lines);
	$seq =~ s/\W+//g;
	$SEQ{$ID} = $seq;
}
close(FILE);
$/ = "\n";

my %ID;
open(FILE, "grep -v '#' $ARGV[1] | awk '{print \$6}' | ") || die("Canot open");
while(<FILE>){
	my ($ID) = /(\S+)/;
	next unless (defined $ID);
	next if(exists $ID{$ID});
	$ID{$ID} = 1;
	print ">$ID\n$SEQ{$ID}\n";
}
close(FILE);