求指导 ~~~半天也没写出来 求大神指导~~

albertkingking

现有文件如下：
>comp0_c0_seq1 len=245 path=[391:0-244]
ATCATTGGAGGCATATTGAAATCTTAGTAGCAGGTATTATTAGAACTTGAGAGAA
>comp0_c0_seq1 len=381 path=[359:0-380]
CGTGGAGAGAAGCCAGCTTGAATTGAATGGTGCAAGACTCGGCTCAA
>comp3_c1_seq1 len=721 path=[717:0-720]
GTGAAGAAACTGTGTTTGGCACCTGTGTCGGCCGCAAAGCACCAATCGTGGAGAGAAGCCACAA
>comp2_c0_seq1 len=468 path=[1:0-467]
ATTAAGAGTTTGGAAGATTAGGTAGCCATTGACAGTTTTACACTGTTTAGAGAAAGACTTATCTTAATTTTGTTTTGGGTTTATTTGGAAATTCAAACTGTATCAAATAAGGGA
>comp1_c1_seq1 len=264 path=[1132:0-263]
GGGTATGAGGATATGGATTTTATTTCATGTATGATTTTAGTCAAATACAAGTATTGTATTAAAAAAAAAAAGTGAGGTG
>comp3_c0_seq1 len=258 path=[481:0-207 110:208-257]
AAGTCAAGAATAAGGTAGAAAAATAGCATACCTGTAATAAGAGTTCCACGGGATGCTGGTGA
>comp2_c0_seq1 len=441 path=[1:0-97 633:98-440]
CACTGATTGATATTTAATATTCCTAAAATTCTTTCATATCTAAGTAGTCTTCCTTCTATTTTTTCTAGGTCTTTCGAAAGGATAACAGAA
>comp3_c0_seq1 len=417 path=[607:0-416]
GTAAGAACTGTCACAAAGCCAAATTCAGGCTATCAGATGGTGCCCAG
>comp1_c1_seq1 len=721 path=[717:0-720]
GTGAAGAAACTGTGTTTGGCACCTGTGTCGGCCGCAAAGCACCAATCGTGGAGAGAAGCCACAA

目的是：按照comp后面紧跟的数字将这些序列文件分组，然后把每一组中的最长的那条序列打印出来，例如按照给出的例子的结果文件是：

>comp0_c0_seq1 len=245 path=[391:0-244]
ATCATTGGAGGCATATTGAAATCTTAGTAGCAGGTATTATTAGAACTTGAGAGAA
>comp3_c0_seq1 len=441 path=[1:0-97 633:98-440]
CACTGATTGATATTTAATATTCCTAAAATTCTTTCATATCTAAGTAGTCTTCCTTCTATTTTTTCTAGGTCTTTCGAAAGGATAACAGAA
>comp2_c0_seq1 len=468 path=[1:0-467]
ATTAAGAGTTTGGAAGATTAGGTAGCCATTGACAGTTTTACACTGTTTAGAGAAAGACTTATCTTAATTTTGTTTTGGGTTTATTTGGAAATTCAAACTGTATCAAATAAGGGA
>comp1_c1_seq1 len=264 path=[1132:0-263]
GGGTATGAGGATATGGATTTTATTTCATGTATGATTTTAGTCAAATACAAGTATTGTATTAAAAAAAAAAAGTGAGGTG

求大神指导~~~~跪谢~~~~{:3_199:}

FTDlinux

这个样子还是比较容易看懂
假设文件叫data

1. #!/usr/bin/perl
   
2. use warnings;
   
3. use strict;
   
4. use Data::Dumper;
   
5. 
   
6. $/ = ">";
   
7. my %hash;
   
8. 
   
9. open IN, "<", "data" or die "Can't open data:$!\n";
   
10. 
    
11. <IN>;
    
12. 
    
13. while (<IN>) {
    
14.         chomp;
    
15.         my @tmp = split /\n/;
    
16.         my $comp = (split /_/,$tmp[0])[0];
    
17.         if (!exists $hash{$comp}) {
    
18.                 $hash{$comp}{ID}  = $tmp[0];
    
19.                 $hash{$comp}{seq} = $tmp[1];
    
20.         } elsif (length $tmp[1] > length ($hash{$comp}{seq})) {
    
21.                 $hash{$comp}{ID}  = $tmp[0];
    
22.                 $hash{$comp}{seq} = $tmp[1];
    
23.         }
    
24. }
    
25. 
    
26. close IN;
    
27. 
    
28. print ">$hash{$_}{ID}\n$hash{$_}{seq}\n" for keys %hash;
    

复制代码

./data.pl

>comp2_c0_seq1 len=468 path=[1:0-467]
ATTAAGAGTTTGGAAGATTAGGTAGCCATTGACAGTTTTACACTGTTTAGAGAAAGACTTATCTTAATTTTGTTTTGGGTTTATTTGGAAATTCAAACTGTATCAAATAAGGGA
>comp3_c1_seq1 len=721 path=[717:0-720]
GTGAAGAAACTGTGTTTGGCACCTGTGTCGGCCGCAAAGCACCAATCGTGGAGAGAAGCCACAA
>comp0_c0_seq1 len=245 path=[391:0-244]
ATCATTGGAGGCATATTGAAATCTTAGTAGCAGGTATTATTAGAACTTGAGAGAA
>comp1_c1_seq1 len=264 path=[1132:0-263]
GGGTATGAGGATATGGATTTTATTTCATGTATGATTTTAGTCAAATACAAGTATTGTATTAAAAAAAAAAAGTGAGGTG

rubyish

如下：

1. #!/usr/bin/perl
   
2. 
   
3. my ( %h, @a );
   
4. open my $d, 'data';
   
5. while (<$d>) {
   
6.     /^.(\w+)/;
   
7.     my $n = <$d>;
   
8.     my $l = length $n;
   
9.     push @a, $1 if !$h{$1};
   
10.     $h{$1} = [ $l, $_ . $n ] if !$h{$1} or $h{$1}[0] < $l;
    
11. }
    
12. print $h{$_}[1] for @a;

复制代码

albertkingking

谢谢指导 ~~~~
我得先看看学习学习 {:3_193:} 回复 2# FTDlinux


   

albertkingking

谢谢啦 ~~~~
学习 学习了~~~{:3_193:} 回复 3# rubyish


   

albertkingking

大神 太简略了 膜拜了~~~回复 3# rubyish


   

albertkingking

太清晰了 对于我们这种菜鸟 学习起来比较舒服  膜拜 了 回复 2# FTDlinux