小菜求助：获取指定序列

sunzhiguolu

回复 10# ahy114

经过你对第二个文件的修改, 发现我对第一个文件与第二个文件之间的 大区间 与 小区间之间的关系有点糊涂了; 最后得出的结果与区间的数值之间有什么关系请再举例说明一下.

ahy114

回复 11# sunzhiguolu

我用下面这个图说明了下

sunzhiguolu

回复 12# ahy114
试下,

1. #!/usr/bin/perl
   
2. use strict;
   
3. use warnings;
   
4. 
   
5. sub load_data{
   
6.     local @ARGV = shift;
   
7.     my $data = shift;
   
8.     my $merge = defined $data;
   
9.     my %hCnt;
   
10.     while (<>){
    
11.         next if (/\A\D/);
    
12.         my ($group, $start, $stop, $type) = /\w+/g;
    
13.         if ($merge){
    
14.             if (exists $data->{$group}{$type}){
    
15.                 $hCnt{$group,$type}++ ? push (@{$data->{$group}{$type}}, [$start, $stop]) : ($data->{$group}{$type}->[0][-1] = $stop);
    
16.             }
    
17.             next;
    
18.         }
    
19.         push (@{$data->{$group}{$type}}, [$start, $stop]);
    
20.     }
    
21.     $data;
    
22. }
    
23. 
    
24. my ($data, $grpID);
    
25. $data = load_data (shift, $data);
    
26. $data = load_data (shift, $data);
    
27. 
    
28. while (<>){
    
29.     if (/\A>(\d)/){
    
30.         $grpID = $1;
    
31.         next;
    
32.     }
    
33.     next if (!exists $data->{$grpID});
    
34.     my $line = $_;
    
35.     foreach my $intervals (values %{$data->{$grpID}}){
    
36.         my $str = "";
    
37.         foreach (@$intervals){
    
38.             my ($start, $stop) = @$_;
    
39.             $start -= 1;
    
40.             $str .= substr ($line, $start, ($stop - $start));
    
41.         }
    
42.         print "$str\n";
    
43.     }
    
44. }
    
45. 
    

复制代码

perl abc.pl a b c
-----------------------------------
1223abc
33aab

ahy114

完美，感谢大神的帮助！！！