Skip to content

Commit ee024fe

Browse files
author
David Jones
committed
Change tabix->query to use tabix->query_full
1 parent f30bcef commit ee024fe

File tree

1 file changed

+7
-5
lines changed

1 file changed

+7
-5
lines changed

perl/bin/utilities/ascatSnpPanelGenerator.pl

Lines changed: 7 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -154,11 +154,13 @@ sub get_regions {
154154
for my $t_range(@top_ranges) {
155155
my ($c_chr, $c_start, $c_end) = @{$t_range};
156156
if(defined $exclude) {
157-
my $iter = $exclude->query(sprintf "%s:%d-%d", $c_chr, $c_start, $c_end);
158-
while(my $record = $iter->next){
159-
my ($start0, $end1) = (split /\t/, $record)[1..2];
160-
push @refined_regions, [sprintf("%s:%d-%d", $c_chr, $c_start, $start0), $c_chr, $c_start, $start0] unless($start0 == 0);
161-
$c_start = $end1;
157+
my $iter = $exclude->query_full($c_chr, $c_start, $c_end);
158+
if(defined($iter)){
159+
while(my $record = $iter->next){
160+
my ($start0, $end1) = (split /\t/, $record)[1..2];
161+
push @refined_regions, [sprintf("%s:%d-%d", $c_chr, $c_start, $start0), $c_chr, $c_start, $start0] unless($start0 == 0);
162+
$c_start = $end1;
163+
}
162164
}
163165
}
164166
push @refined_regions, [sprintf("%s:%d-%d", $c_chr, $c_start, $c_end), $c_chr, $c_start, $c_end];

0 commit comments

Comments
 (0)