use strict;
use warnings;
use DB_File;


my $myhashfile = "hash.$$";
tie my %hash1, "DB_File", $myhashfile, O_RDWR|O_CREAT, 0666, $DB_HASH
    or die "cannot open file $myhashfile: $!";


open(OUT,">output.out")or die "can not open";
open(my $fh1, "<", "file1.db") or die "file1.db: $!";


foreach (<$fh1>){
    chomp;
    my @dbinfo = split(/\s+/);
    $hash1{"$dbinfo[1]#$dbinfo[2]"} = "$dbinfo[4]##$dbinfo[0]";
}
close($fh1);


my $c=0;
open(my $fh2, "<", "file2.query") or die "file2.query: $!";
my @snplist;
foreach (<$fh2>) {
chomp($_);@snplist=();

    my @queryinfo = split(/[\s>]+/);
    my $values = $hash1{"$queryinfo[0]#$queryinfo[1]"};
    my ($variant,$rs)=split("##",$values);

    my $flag_left=0;my $flag_right=0;
    @snplist=split("/",$variant);

    if(defined($variant)){ 
	foreach my $lis(@snplist){
		if($queryinfo[2] eq $lis){$flag_left=1;}
		elsif($queryinfo[3] eq $lis){$flag_right=1;}
	}

        if(($flag_left == 1) && ($flag_right == 1)){
		print OUT "$_\t$rs\n";
		$c++;
		$flag_left=0;$flag_right=0;
    	}
  }

  else{
	print OUT "$_\t-\n";
  }
}


print OUT "\nTotal number of HITS: $c\n";
close($fh2);

untie %hash1;
unlink($myhashfile);
$myhashfile=();