Skip to content
Snippets Groups Projects
Commit 0a9346d3 authored by Daniel Rios's avatar Daniel Rios
Browse files

check there are differences between strainSlice and Slice

parent d1e0c62c
No related branches found
No related tags found
No related merge requests found
......@@ -210,8 +210,11 @@ sub _get_indels{
#go throuh all the strains getting ONLY the indels (length_diff <> 0)
my @indels;
foreach my $strainSlice (@{$self->strains}){
my @results = grep {$_->length_diff != 0} @{$strainSlice->get_all_differences_Slice()};
push @indels, @results;
my $differences = $strainSlice->get_all_differences_Slice(); #need to check there are differences....
if (defined $differences){
my @results = grep {$_->length_diff != 0} @{$differences};
push @indels, @results;
}
}
#need to overlap the gaps using the RangeRegistry module
my $range_registry = Bio::EnsEMBL::Mapper::RangeRegistry->new();
......
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