#!/usr/bin/perl

open(IN,$ARGV[0]);
while(<IN>) {
/(\S+)/;
$get{$1} = 1;
}
close(IN);

open(IN,"get_indels.out.hd");

while(<IN>) {
/^(\S+)\s+(.*)/;
$gene = $1;
@fields = split(/\s+/,$2);
$sum = 0;
	foreach $f (@fields) {
	$sum = $sum +  $f;
		if($get{$gene} == 1) {
		push(@deletions,$f);
		}
	}
$sum_del{$gene} = $sum;
}
close(IN);

open(IN,"get_indels.out.hi");

while(<IN>) {
/^(\S+)\s+(.*)/;
$gene = $1;
@fields = split(/\s+/,$2);
$sum = 0;
	foreach $f (@fields) {
	$sum = $sum +  $f;
		if($get{$gene} == 1) {
		push(@insertions,$f);
		}
	}
$sum_ins{$gene} = $sum;
}
close(IN);

foreach $g (keys(%get)) {
#print  "$g $sum_ins{$g} $sum_del{$g}\n";
}

print join(" ",@deletions);
#print join(" ",@insertions);
