Skip to content
Snippets Groups Projects
Commit 5aa3583b authored by Yuan Chen's avatar Yuan Chen
Browse files

add upstream/downstream flanking 5kb around transcript when consider a...

add upstream/downstream flanking 5kb around transcript when consider a variation is involved in the transcript
parent 4ac75038
No related branches found
No related tags found
No related merge requests found
......@@ -172,10 +172,19 @@ sub calculate_same_codon{
sub type_variation {
my $tr = shift;
my $var = shift;
my $UPSTREAM = 5000;
my $DOWNSTREAM = 5000;
if (!$var->isa('Bio::EnsEMBL::Variation::ConsequenceType')){
throw("Not possible to calculate the consequence type for ",ref($var)," : Bio::EnsEMBL::Variation::ConsequenceType object expected");
}
if (($var->start < $tr->start - $UPSTREAM) || ($var->end > $tr->end + $DOWNSTREAM)){
#since the variation is more than UPSTREAM and DOWNSTREAM of the transcript, there is no effect in the transcript
return [];
}
my $tm = $tr->get_TranscriptMapper();
my @coords = $tm->genomic2cdna($var->start,
$var->end,
......
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