open (IN1,"<","$ARGV[0]");#chr_001 521     829     H2A.Z-NHA_peak_1        21      .       2.05858 3.0958  2.13255 chr_001 255     1182    TTHERM_00161861 .       +
open (IN2,"<","$ARGV[1]");#genebody.bed
open (OUT1,">","$ARGV[2]");#Plus
open (OUT2,">","$ARGV[3]");#Minus
while (defined ($_=<IN1>)){
	if ($_=~/\s+(TTHERM_\d+)\s+/){
		$hash{$1}=1;
		}}
while (defined ($_=<IN2>)){
	if ($_=~/(TTHERM_\d+)/){
		if (exists $hash{$1}){
		print OUT1 "$_"}
		else{print OUT2 "$_"}}}