#!/usr/bin/perl -w

my %cluster21;
open(FILE, "grep '^21' $ARGV[0] |");
while(<FILE>){
	my ($plasmid) = /^\d+\s+(\S+)/;
	$cluster21{$plasmid} = 1;
}
close(FILE);

$/ = '>';
open(FILE, "$ARGV[1]");
while(<FILE>){
        next unless (/\w/);
        my @lines = split(/\n/, $_);
        my ($ID) = $lines[0] =~ /^(\S+)/;
        shift(@lines);
        my $Seq = join('', @lines);
        $Seq =~ s/\W+//g;
	my ($plasmid) = $ID =~ /(\S+)XY/;
	next unless (exists  $cluster21{$plasmid});
	print ">$ID\n$Seq\n";
}