求指教！遍历序列的问题！！

445236708

我把测试数据换成我的数据后，就输出氨基酸序列了{:2_169:} ，没按测试的情况运行回复 7# pitonas


   

pitonas

你数据:
奇数行序列的名字
偶数行序列

name1
GCGCAAGC
name2
GGCCAAG
name3
CCCCAAA

你贴个 前面 6行 8行的数据, 就行了

445236708

__DATA__
>scaffold16 36.5
AGCTGTTCAACGTAGCTGGTCAGGCTTAACGGTTACAGTACAGTACATGACTAGAATTGG
>scaffold35 33.2
ACGGTTACAGTACAGAGCTTCGATCGTTAGCTTGCTACGTCGTAGCTTTACGTACGTAGC
>scaffold19 33.2
ACGATCGATGCTAGATGCTACGATGATCGTGTGCTGTCGCATGCTAGCTACGTAGCTATA
>scaffold28 31.8
ACGGTTTTTATCGATTCTATGATCTCGTGATGCTCGTCGTAGCTGCATCGATGTCTGTGC
>scaffold31 33.5
ACGGGTAGCTAGTCGTTGCTGACTGTCGATTAGCTACGTACCGTAGTCGGTAGACTGTAC
>scaffold34 34.2
AGTAGTGTATGTCGTAGCCTCGATCAGCTGCTAGCTGATGGTCGTGTGGTCTCTCTTTCC回复 12# pitonas


   

pitonas

{:2_171:}
你发表于 1 楼
那个源文件中有82行，
偶数行是序列的名字，
奇数行是序列

1. #!/usr/bin/perl
   
2. 
   
3. @ARGV = '/tmp/plum_0630.scafSeq.FG';  # filename
   
4. my $length = 250;
   
5. my %data;
   
6. 
   
7. while (<>) {
   
8.     my ($name) = /^>(\w+)/;
   
9.     chomp( my $seq = <> );
   
10.     my @GC = map { tr/GC// } $seq =~ /(.{1,$length})/g;
    
11.     $data{$name} = [@GC];
    
12. }
    
13. 
    
14. while ( my ( $name, $gc ) = each %data ) {
    
15.     print "$name\t@$gc\n";
    
16. }

复制代码

回复 13# 445236708


   

445236708

额，不好意思，理解上的差错 ，不过确实是感谢你这么耐心的解答回复 14# pitonas


   

Hugo801122

我的hash table啊，到现在我都没搞明白！