求助：取最长转录本脚本

tiger4567

崇拜回复 8# huang6894


   

tiger4567

回复 8# huang6894

1. #!/usr/bin/perl -w
   
2. open (DATA,"DATA.txt") or die $!;
   
3. my(%h1, %h2);
   
4. $/="\>";
   
5. <DATA>;
   
6. while(<DATA>){
   
7.         chomp;
   
8.         my($g_t, @fastq) = (split /\n/);
   
9.         my($gene, $tran) = (split /\|/,$g_t);
   
10.         push @{$h2{$tran}},@fastq;
    
11.        $h1{$gene} ||= $tran ;
    
12.         $h1{$gene} = $tran  if length(join("",@{$h2{$h1{$gene}}})) < length(join("",@fastq));        
    
13. }
    
14. $/="\n";
    
15.      map{print ">$_|$h1{$_}\n". join("\n",@{$h2{$h1{$_}}}).$/;}(keys %h1);

复制代码

可以这样子。。。