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";
}
|
