Khader Shameer

# Perl code to report the common region between chromosome regions
# Khader Shameer
   new();
    sub new{
            my @regions = qw(35388351-55114603 37310186-52396773 38511024-48638500 32729632-48638502 17446360-51119526);
            #my @regions = qw(24328000-29166946 25388351-27114603 22310186-25239677 28511024-29638159 23729632-26385029);
            my(@starts, @ends);
            #my @regions = qw(1-20 1-80 10-20 5-25);
                    foreach my $region(@regions)
                                    {
                                    my($start, $end) = split('-', $region);
                                    push @starts, $start;
                                    push @ends, $end;
                                    }
                                             my ($left_max) = reverse sort{$a <=> $b} @starts;
                                             my ($right_min) = sort{$a <=> $b} @ends;
                                                    foreach my $each_region(@regions)
                                                            {
                                                            my ($left, $right) = split('-', $each_region);
                                                            next if $left != $left_max || $right != $right_min;
                                                            print "$left -  $right \n";
                                                            }
            }