#!/bin/perl
use warnings;
use strict;

my $inputfile=$ARGV[0];
#my $inputfile2=$ARGV[1];

open (file2,">clusters0h_over5b.bed") or die "$!\n";

my %matrix1;my %matrix2;my %matrix3;my %matrix4;

my $temp2;
my $temp3;
my $temp4;
my $temp1;

my $var2=0;
my $var1=0;

open (INPUT_FILE,$inputfile);

while(my $line1=<INPUT_FILE>)  
 {   
    
 @{$matrix1{$temp1}}=split(/\s/,$line1);
    $var1=${$matrix1{$temp1}}[2]-${$matrix1{$temp1}}[1];
   # $var2=${$matrix1{$temp1}}[3]/2;
  if ( $var1 >= 20 )
   { 
    if ( ${$matrix1{$temp1}}[3] >=2)
     {
print file2 "${$matrix1{$temp1}}[0]\t${$matrix1{$temp1}}[1]\t${$matrix1{$temp1}}[2]\t${$matrix1{$temp1}}[3]\t${$matrix1{$temp1}}[4]\t${$matrix1{$temp1}}[5]\n";
     }
   }
   
  $temp1++;

  }

close (file2);