[Bio] / FigKernelPackages / PropagateGBMetadata.pm Repository:
ViewVC logotype

Diff of /FigKernelPackages/PropagateGBMetadata.pm

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.2, Wed Nov 18 22:20:02 2015 UTC revision 1.3, Mon Feb 29 17:37:45 2016 UTC
# Line 111  Line 111 
111          next unless $f->{type} =~ /rna/i;          next unless $f->{type} =~ /rna/i;
112          next unless exists $f->{genbank_feature};          next unless exists $f->{genbank_feature};
113    
114            #
115            # From Set::IntervalTree doc:
116            # All intervals are half-open, i.e. [1,3), [2,6), etc.
117            # Thus we bump the right endpoint.
118            #
119            $max++;
120    
121          my $tree = $trees{$ctg};          my $tree = $trees{$ctg};
122          $tree->insert($f->{id}, $min, $max);          $tree->insert($f->{id}, $min, $max);
123      }      }
# Line 127  Line 134 
134    
135          my($ctg, $min, $max, $dir, $len) = GenomeTypeObject::bounds($f);          my($ctg, $min, $max, $dir, $len) = GenomeTypeObject::bounds($f);
136          my $tree = $trees{$ctg};          my $tree = $trees{$ctg};
137          my $overlap = $tree->fetch($min, $max);          my $overlap = $tree->fetch($min, $max+1);
138          my @overlap = grep { $_ ne $fid } @$overlap;          my @overlap = grep { $_ ne $fid } @$overlap;
139          #          #
140          # @overlap contains the genbank feature we are mapping from.          # @overlap contains the genbank feature we are mapping from.

Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3

MCS Webmaster
ViewVC Help
Powered by ViewVC 1.0.3