perl求助！

yang7473453

NW_010786190.1         46474        46670                 XP_010360621.1
NW_010786190.1         55641        55860                 XP_010360621.1
NW_010786190.1         66091        66240                 XP_010360621.1
NW_010786190.1         66489        66704                 XP_010360621.1
NW_010786203.1         9752        9838                 XP_010370913.1
NW_010786203.1         13535        13735                 XP_010370913.1
NW_010786203.1                     14755        14968                 XP_010370913.1
NW_010786203.1             9752        9838                 XP_010377713.1
NW_010786203.1         13535        13735                 XP_010377713.1

现有文件如上，当最右边一样时，最左边也是一样的，想要将文件变成以下格式

NW_010786190.1；XP_010360621.1；46474-46670  55641-55860  66091-66240  66489-66704
NW_010786203.1；XP_010370913.1；9752-9838  13535-13735  14755-14968
NW_010786203.1；XP_010377713.1；9752-9838  13535-13735

求perl脚本

huang6894

1. #!/usr/bin/perl -w
   
2. my %hash;
   
3. while(<DATA>){
   
4.         chomp;
   
5.         my @tmp = split;
   
6.         push @{$hash{$tmp[0].";".$tmp[-1].";"}},$tmp[1]."-".$tmp[2];
   
7. }
   
8. foreach my $keys(sort keys %hash){
   
9.         print join("\t",$keys,@{$hash{$keys}}).$/;
   
10. }
    
11. __DATA__
    
12. NW_010786190.1         46474        46670                 XP_010360621.1
    
13. NW_010786190.1         55641        55860                 XP_010360621.1
    
14. NW_010786190.1         66091        66240                 XP_010360621.1
    
15. NW_010786190.1         66489        66704                 XP_010360621.1
    
16. NW_010786203.1         9752        9838                 XP_010370913.1
    
17. NW_010786203.1         13535        13735                 XP_010370913.1
    
18. NW_010786203.1                     14755        14968                 XP_010370913.1
    
19. NW_010786203.1             9752        9838                 XP_010377713.1
    
20. NW_010786203.1         13535        13735                 XP_010377713.1

复制代码