Skip to content
Snippets Groups Projects
Commit bcec7c4f authored by Stephen Searle's avatar Stephen Searle
Browse files

Fix for subseq of sequence on small seq_regions (< 250k) which was

failing because it was requesting a sub slice which was longer than
the seq_region in the new caching code
parent 60eaf459
No related branches found
No related tags found
No related merge requests found
......@@ -279,8 +279,13 @@ sub subseq {
my $repeats = [];
my $subslice = $seq_region_slice->sub_Slice( ($block_min << $BLOCK_PWR)+1,
(($block_max+1)<<$BLOCK_PWR ));
my $sub_start = ($block_min << $BLOCK_PWR)+1;
my $sub_end = ($block_max+1)<<$BLOCK_PWR;
if ($sub_end > $seq_region_slice->length) {
$sub_end = $seq_region_slice->length ;
}
my $subslice = $seq_region_slice->sub_Slice( $sub_start, $sub_end);
foreach my $l (@$logic_names) {
push @{$repeats}, @{$subslice->get_all_RepeatFeatures($l)};
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment