隔行的数据处理

gooderpan

有一个文件file，如下>Gene.1::NM_198493::g.1::m.1 type:5prime_partial len:287 gc:universal NM_198493:2-862(+)
LRRCYGNRGGRGPGVVFLELMESEGPPESESSEFFSQQEEENEEEEAQEPEETGPKNPLLQPALTGDVEGLQKIFEDPEN
>Gene.2::NM_198493::g.2::m.2 type:complete len:153 gc:universal NM_198493:1865-2323(+)
MYEQFSFNFLKQKNTAITFSHLFGRPSGCLWSSTKKTCVCSFFFFTLL
>Gene.3::NM_198493::g.4::m.4 type:complete len:118 gc:universal NM_198493:1954-2307(+)
MEQYKEDLCLLLFFFYPPSPFSLFPSLPLLLPLPSSPFPSLPLPSPPFLSLPLPSPPFLSPPLPSSPLPSPSSSLSSSSSLSNSFFFLSYLSLSLHPRNLLDRVCLSCPGWNAVE
>Gene.4::NM_173852::g.6::m.6 type:complete len:163 gc:universal NM_173852:27-515(+)
MRIANRTRFSSPFLARGAGWTHGRGMMVVGTGTSLALSSLLSLLLFAGMQMYSRQLASTEWLTIQGGLLGSGLFVFSLTAFNNLENLVFGKGFQAKIFPEILLCLA
>Gene.9::NM_013945::g.9::m.9 type:complete len:519 gc:universal NM_013945:599-2155(+)
MAALPGTVPRMMRPAPGQNYPRTGFPLEVSTPLGQGRVNQLGGVFINGRP
>Gene.13::NM_013945::g.13::m.13 type:complete len:102 gc:universal NM_013945:1821-2126(+)
MAAWTRPPPSQPAAASGPTPSSQETACPPPRPTAHPPTAPPATAWTPWPAISTASTARVSAWCPGRPPSPFLLPPPGPPACLWRATRWCQGGECPFHRWKN

该文件内容是每两行为一组，根据“>”后面的内容，选出相同NM编号下“len”最大的那一组，想要得到的结果如下：
>Gene.1::NM_198493::g.1::m.1 type:5prime_partial len:287 gc:universal NM_198493:2-862(+)
LRRCYGNRGGRGPGVVFLELMESEGPPESESSEFFSQQEEENEEEEAQEPEETGPKNPLLQPALTGDVEGLQKIFEDPEN
>Gene.4::NM_173852::g.6::m.6 type:complete len:163 gc:universal NM_173852:27-515(+)
MRIANRTRFSSPFLARGAGWTHGRGMMVVGTGTSLALSSLLSLLLFAGMQMYSRQLASTEWLTIQGGLLGSGLFVFSLTAFNNLENLVFGKGFQAKIFPEILLCLA
>Gene.9::NM_013945::g.9::m.9 type:complete len:519 gc:universal NM_013945:599-2155(+)
MAALPGTVPRMMRPAPGQNYPRTGFPLEVSTPLGQGRVNQLGGVFINGRP

请教大神们该怎么办呀？？？？？？？谢谢啦！！！

fcl21721

$/=undef;
$content = <DATA>;
my %hash;
while($content =~ /(Gene.+\n.+)/g) {
        #push @array, $1;
        my $content = $1;
        ($key, $length) = $content =~ /(NM_\d+).+len\d+)/;
        push @array, [$content, $key, $length];
       
        $hash{$key} //= $length;
        if($hash{$key} < $length) {
                $hash{$key} = $length;
        }
}
@target = map {$_->[0]} grep { if($hash{$_->[1]} == $_->[2]){1} else{0} } @array;
print join"\n", @target;
__DATA__
>Gene.1::NM_198493::g.1::m.1 type:5prime_partial len:287 gc:universal NM_198493:2-862(+)
LRRCYGNRGGRGPGVVFLELMESEGPPESESSEFFSQQEEENEEEEAQEPEETGPKNPLLQPALTGDVEGLQKIFEDPEN
>Gene.2::NM_198493::g.2::m.2 type:complete len:153 gc:universal NM_198493:1865-2323(+)
MYEQFSFNFLKQKNTAITFSHLFGRPSGCLWSSTKKTCVCSFFFFTLL
>Gene.3::NM_198493::g.4::m.4 type:complete len:118 gc:universal NM_198493:1954-2307(+)
MEQYKEDLCLLLFFFYPPSPFSLFPSLPLLLPLPSSPFPSLPLPSPPFLSLPLPSPPFLSPPLPSSPLPSPSSSLSSSSSLSNSFFFLSYLSLSLHPRNLLDRVCLSCPGWNAVE
>Gene.4::NM_173852::g.6::m.6 type:complete len:163 gc:universal NM_173852:27-515(+)
MRIANRTRFSSPFLARGAGWTHGRGMMVVGTGTSLALSSLLSLLLFAGMQMYSRQLASTEWLTIQGGLLGSGLFVFSLTAFNNLENLVFGKGFQAKIFPEILLCLA
>Gene.9::NM_013945::g.9::m.9 type:complete len:519 gc:universal NM_013945:599-2155(+)
MAALPGTVPRMMRPAPGQNYPRTGFPLEVSTPLGQGRVNQLGGVFINGRP
>Gene.13::NM_013945::g.13::m.13 type:complete len:102 gc:universal NM_013945:1821-2126(+)
MAAWTRPPPSQPAAASGPTPSSQETACPPPRPTAHPPTAPPATAWTPWPAISTASTARVSAWCPGRPPSPFLLPPPGPPACLWRATRWCQGGECPFHRWKN

sunzhiguolu

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. local $/ = "\n>";
   
6. my ($leng, $str) = 0;
   
7. while (<>){
   
8.     chomp;
   
9.     s/\A>// if (!$leng);
   
10.     if ($. % 2 and /\blen:(\d+)/){
    
11.         ($leng, $str) = ($1, $_);
    
12.         next;
    
13.     }
    
14.     if (/\blen:(\d+)/){
    
15.         print $leng > $1 ? ">$str\n" : ">$_\n";
    
16.     }
    
17. }
    

复制代码

华小飞_Perl

换个角度来写一个~

1. #!/usr/bin/perl
   
2. 
   
3. use warnings;
   
4. use strict;
   
5. use List::Util qw(max);
   
6. 
   
7. my (%gene_info, %id_max);
   
8. my ($id, $val);
   
9. while (<>) {
   
10.         chomp;
    
11.         if (/>.*?(NM_\d+).*?len:(\d+)/) {               
    
12.                 ($id, $val) = ($1, $2);
    
13.                 $gene_info{$id}{$val} = $_;
    
14.                 next;
    
15.         } elsif (/^[^>]+/m) {
    
16.                 $id_max{$id}{$val} = $_;
    
17.         }
    
18. }
    
19. 
    
20. for my $name ( keys %id_max ) {
    
21.         my $max = max keys %{ $id_max{$name} };
    
22.         print $gene_info{$name}{$max}, "\n";
    
23.         print $id_max{$name}{$max}, "\n";       
    
24. }

复制代码