clstr_reduce.pl

#!/usr/bin/env perl


$file90 = shift;
$segs = shift;
$reduce_rate = shift;

@clstr_nos = ();
@segs = split(/,/, $segs);
@segs_no = ();
@no2seg_idx = ();
for ($i=0; $i<@segs; $i++) {
  $seg = $segs[$i];
  if ($seg =~ /-/) { ($b, $e) = split(/-/, $seg); }
  else             {  $b= $e  = $seg; }
  for($j=$b; $j<=$e; $j++) { $no2seg_idx[$j]=$i; }
  $segs_no[$i] = 0;
}


open(TMP, $file90) || die "Can not open file";
$readin = 0;
$this_no = 0;
$this_clstr = "";
$ll;
while($ll=<TMP>) {
  if ($ll =~ /^>/ ) {
    if ($readin) { 
      $this_seg = $no2seg_idx[$this_no];
      if (($segs_no[$this_seg] % $reduce_rate) == 0) {
        print $this_clstr;
      }
      $segs_no[$this_seg]++;
    }
    $this_no=0;
    $this_clstr = $ll;
  }
  else {
    $this_clstr .= $ll;
    $readin = 1; $this_no++;
  }
}
close(TMP);

    if ($readin) { 
      $this_seg = $no2seg_idx[$this_no];
      if (($segs_no[$this_seg] % $reduce_rate) == 0) {
        print $this_clstr;
      }
      $segs_no[$this_seg]++;
    }