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

quick fix to remove indels from strain/individual sequence

parent 5cfce118
No related branches found
No related tags found
No related merge requests found
......@@ -249,6 +249,7 @@ sub seq {
my @variation_features_ordered = sort {$b->start() <=> $a->start()} @{$self->{'alleleFeatures'}} if (defined $self->{'alleleFeatures'});
foreach my $vf (@variation_features_ordered){
next if ($vf->start != $vf->end); #quick fix to remove indels from sequence
$vf->apply_edit($reference_sequence); #change, in the reference sequence, the vf
}
#need to find coverage information if different from reference
......
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