sub versatility{ my(%genome1,%genome2)=@_; for my $key (keys %genome1) { my @genome1=$genome1{$key}=~ /(\d+)/g; my @genome2=$genome2{$key}=~ /(\d+)/g; my $m=($genome1[1]-$genome2[1]); my $n=($genome1[0]-$genome2[0]); if( $n==0){ print "number of occurance for $key in two genome are same\n";} else{ my $slope=$m/$n; print "domain versatility= $slope\n";} } }