#! /usr/bin/perl my @row1=(); my @row2=(); my $ini1=0; my $ini2=0; my $id1=""; my $id2=""; my $fam1=""; my $fam2=""; my $c1=0; my $c2=0; while () { chomp($_); if (($c1==0)&&($c2==0)){ @row1=split (/\t/,$_); } if((!$row1) && ($_=~/Seleno/)){ @row1=split (/\t/,$_); @fam1=split (/\(/,$_); @fam1=split (/\)/,$fam1[1]); $fam1=$fam1[0]; $c1++; } if(($_!~/Seleno/)&&(!$row2)){ @row2=split(/\t/,$_); @fam2=split(/\./,$_); $fam2=$fam2[0]; $c2++; } if ((@row1)&&(@row2)){ $chr1=$row1[2]; $chr2=$row2[2]; $ini1=$row1[3]; $ini2=$row2[3]; $id2=$row2[0]; $id1=$row1[0]; if((($chr1 == $chr2) || ($chr1 eq $chr2))&&(($ini1>=$ini2-75000)&&($ini1<=$ini2+75000))){ print "MATCH:\t$id1 $id2 $fam1 $fam2 $chr1 $chr2 $ini1 $ini2\n"; }else{ print "NO:\t$id1 $id2 $fam1 $fam2 $chr1 $chr2 $ini1 $ini2\n"; } $row1=(); $row2=(); } }