#!/usr/bin/perl 


use warnings;
use strict;

my @outs;
my %grupos1;
my %grupos2;
my @split; 
#my $out;	
my @element2; 
my @element1; 


open my $IN , '<', " $ARGV[0]"; 
open my $OUT, '>', "$out.out" or die "open failed: $!\n";



#########################----

while  (my $line = <$IN>){
	chomp $line;
	my @split = split /\t/, $line;
	my $id1 = join ("\t", @split);

	if (!defined $grupos1{$id1}) {
		$grupos1{$id1}= [];
	} #close defined hash %grupos1 if

	push @{$grupos1{$id1}}, [$split[0],$split[1],$split[2],$split[3],$split[4]];  

}# close while
close $IN;


#-----------------------------------
open my $IN2 , '<', "$ARGV[1]"; 

while (my $line2=<$IN2>) {
	chomp $line2;                  
	my @split2 = split /\t/, $line2;
	my $id2 = join ("\t", @split); 

	if (!defined $grupos2{$id2}) {
			$grupos2{$id2}= [];

	} #close defined hash %grupos2 if

	push @{$grupos2{$id2}}, [$split2[12],$split2[1], $split2[10], $split2[11]];  

} #close while

close $IN2;


foreach my $key1 ( sort keys %grupos1 ) {
	for my $element1 ( @{ $grupos1{$key1} } ) {    #derefrence

		foreach my $key2 ( sort keys %grupos2 ) {
			for my $element2 ( @{ $grupos2{$key2} } ) {

				if ( $element1->[0] eq $element2->[0]){

					my $new_start= $element1->[2]+$element2->[2]; 
					my $new_end= $element1->[3]+$element2->[3];  
					
					print  $OUT, "$element1->[1]\t$new_start\t$new_end\t$element2[1]\t.\t$element1[4].\t.$element2[0]\n";
					
				} ## end if
			}    #end foreach key1
		}
	}    #end foreach key2
} #end foreach key1
print "$out is done\n";
%grupos1 = (); #empty hash
%grupos2 = (); #empty hash
close $OUT;
#} #close genome