# List of headers # BLAST file # Output file (BLAST format) use strict; # Specify I/O my $blastfile = $ARGV[0]; my $headerfile = $ARGV[1]; my $outfile = $ARGV[2]; # open blast file open (BLASTFILE, $blastfile) or die "Cannot open BLAST file\n"; # open header file open (HEADERFILE, $headerfile) or die "Cannot open header file\n"; # open output file open(OUTPUT, ">$outfile"); my @headerline = ; my $headerline = join ( '', @headerline); while (my $blastline = ) { my $blastcut = substr $blastline, 7; if ($headerline =~ /\Q$blastcut\E/) { print $blastcut; } } exit;