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

added condition to check for regulatory features only if there is a FuncGen database

parent 9fe37c56
No related branches found
No related tags found
No related merge requests found
......@@ -255,9 +255,8 @@ sub type_variation {
}
##to find if a SNP is overlapping with the transcript of the regulatory region, also the SNP should be within 5kb on both side of the transcript, then check whether they overlapping
if ($tr and ref($tr) and $tr->isa('Bio::EnsEMBL::Transcript')) {
my $dbFunc = $tr->adaptor->db->get_db_adaptor("funcgen");
my $dbFunc = $tr->adaptor->db->get_db_adaptor("funcgen");
if ($tr and ref($tr) and $tr->isa('Bio::EnsEMBL::Transcript') and (defined $dbFunc)) {
my $efa = $dbFunc->get_ExternalFeatureAdaptor();
my $rfa = $dbFunc->get_RegulatoryFeatureAdaptor();
my @rf;
......
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