Skip to content
Snippets Groups Projects
Commit 03e0a891 authored by Ian Longden's avatar Ian Longden
Browse files

Fix for multiple assemblies breaks haplotype projection so removed

parent 2a740284
No related branches found
No related tags found
No related merge requests found
......@@ -223,10 +223,17 @@ sub map_coordinates{
my $self_coord = $pair->{$from};
my $target_coord = $pair->{$to};
if($self_coord->{'start'} < $start){
$start = $orig_start;
$rank++;
}
#
# But not the case for haplotypes!! need to test for this case???
# so removing this till a better solution is found
#
#
# if($self_coord->{'start'} < $start){
# $start = $orig_start;
# $rank++;
# }
if(defined($last_target_coord) and $target_coord->{'id'} ne $last_target_coord){
......
......@@ -808,10 +808,11 @@ sub project {
my $coord_end = $coord->end();
my $length = $coord_end - $coord_start + 1;
if( $last_rank != $coord->rank){
$current_start = 1;
}
$last_rank = $coord->rank;
# if( $last_rank != $coord->rank){
# $current_start = 1;
# print "LAST rank has changed to ".$coord->rank."from $last_rank \n";
# }
# $last_rank = $coord->rank;
#skip gaps
if($coord->isa('Bio::EnsEMBL::Mapper::Coordinate')) {
......
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