Ubuntu中文论坛

将这个程序编写完整，这个是我的作业，不会做，就来论坛求助了，呵呵

统计Canis_familiaris.seq中各蛋白序列的长度和各氨基酸的个数。编号为gi|7387799的蛋白序列中K和C的个数分别是多少？
提示：1)用一函数count统计氨基酸个数，调用形式为count($seq, $letter)
sub count{
my ($seq, $letter) = @_;
my $num = s/$letter//gi;
return $num;
}
2)用函数getletters得到序列中氨基酸字母
sub getletters{
my ($protein) = @_;
my %letters;
my $len = length($protein);
foreach (0..$len-1){
my $letter = substr($protein,$_,1);
$letters{$letter} = 1;
}
return keys %letters;
}

原文件Canis_familiaris.seq什么样子哦。
长沙什么学校开的perl啊。