作者chhuang (我是瓶男我很难懂)
看板perl
标题Re: [请益] 有关写入档案
时间Sat Jun 10 02:04:19 2006
use Bio::SeqIO;
use Bio::Tools::SeqWords;
my $in = Bio::SeqIO->new( -file => "prot.seq", -format => 'Fasta' );
while ( my $seq = $in->next_seq() ) {
print $seq->id . "\n";
my $hash_ref = Bio::Tools::SeqWords->count_words( $seq, 1 );
my %hash = %$hash_ref;
foreach my $key ( sort keys %hash ) {
print "$key $hash{$key}\n";
}
}
OUTPUT:
randomseq1
A 10
C 9
D 14
E 13
F 9
G 6
H 10
I 10
K 10
L 10
M 11
N 11
P 4
Q 8
R 13
S 10
T 12
V 12
W 9
Y 9
--
我是瓶男~我很难懂!
http://blog.yam.com/chhuang
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 61.62.84.89
※ 编辑: chhuang 来自: 61.62.84.89 (06/10 02:05)