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

added hack to not return AlleleFeatures if StrainSlice is the fake reference one

parent 1e4772ec
No related branches found
No related tags found
No related merge requests found
......@@ -341,6 +341,16 @@ sub get_all_AlleleFeatures_Slice{
my $self = shift;
my $with_coverage = shift;
my $variation_db = $self->adaptor->db->get_db_adaptor('variation');
unless($variation_db) {
warning("Variation database must be attached to core database to " .
"retrieve variation information" );
return '';
}
my $indAdaptor = $variation_db->get_IndividualAdaptor();
my $ref_name = $indAdaptor->get_reference_strain_name;
return [] if ($self->strain_name eq $ref_name);
$with_coverage ||= 0; #by default, get all AlleleFeatures
if ($with_coverage == 1){
my $new_allele_features = $self->_filter_af_by_coverage($self->{'alleleFeatures'});
......
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