Mercurial > repos > plus91-technologies-pvt-ltd > softsearch
view 2.4/script/reduce_redundancy.pl @ 14:e8620ee91ca1 draft
Uploaded
author | plus91-technologies-pvt-ltd |
---|---|
date | Fri, 30 May 2014 11:00:16 -0400 |
parents | e3609c8714fb |
children |
line wrap: on
line source
open(BUFF,"$ARGV[0]") or die "no input file found\n"; $range="$ARGV[1]"; my %hash; my %store; $prev_chr=""; $next=0; while(<BUFF>) { chomp($_); #print "$.\n"; if($_ !~ m/^#/) { @array=split("\t",$_); $chr=$array[0]; $pos=$array[1]; $value=$array[@array-1]; if($prev_chr ne $chr ) { if($prev_chr ne "") { foreach $key (sort {$hash{$b} <=> $hash{$a} } keys %hash) { print "$store{$key}\n"; last; } } $next = $pos+$range; undef(%hash); undef(%store); } if($next< $pos) { foreach $key (sort {$hash{$b} <=> $hash{$a} } keys %hash) { print "$store{$key}\n"; last; } $next = $pos+$range; undef(%hash); undef(%store); } if($value eq "NA") { $hash{$chr." ".$pos." ".$.}=0; } else { $hash{$chr." ".$pos." ".$.}=$value; } $store{$chr." ".$pos." ".$.}=$_; } else { print $_."\n"; } $prev_chr = $chr; } foreach $key (sort {$hash{$b} <=> $hash{$a} } keys %hash) { print "$store{$key}\n"; last; }