Codebase list cd-hit / a68a11ac-1c8a-46c5-891d-ba2259aab301/main clstr_select_rep.pl
a68a11ac-1c8a-46c5-891d-ba2259aab301/main

Tree @a68a11ac-1c8a-46c5-891d-ba2259aab301/main (Download .tar.gz)

clstr_select_rep.pl @a68a11ac-1c8a-46c5-891d-ba2259aab301/main — raw · history · blame 

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
#!/usr/bin/perl

#my $by = shift;
my $min;
my $max;
#if ($by eq "size") {
  $min = shift;
  $max = shift;
#}

$rep = "";
$no = 0;

while($ll=<>){
  if ($ll =~ /^>/) {
    if (($no >= $min) and ($no <= $max)) {
      print "$rep\n";
    }
    $rep = "";
    $no = 0;
  }
  else {
    chop($ll);
    if ($ll =~ /\*$/) {
      $rep = "";
      if ($ll =~ /\d+(aa|nt), >(.+)\.\.\./) {
        $rep = $2;
      }
      else {
        die "format error $ll";
      }
    }
    $no++;
  }
}
    if (($no >= $min) and ($no <= $max)) {
      print "$rep\n";
    }