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

Added calculation on the flyu of new consequence type: SARA (Same As Reference Allele)

parent 1d0ccb0c
No related branches found
No related tags found
No related merge requests found
......@@ -98,10 +98,15 @@ sub get_all_ConsequenceType {
}
my $consequence_type = Bio::EnsEMBL::Variation::ConsequenceType->new($transcript->dbID(),'',$allele->start, $allele->end, $transcript->strand, [$string]);
if ($allele->start == 16165){
print "hello";
#calculate the consequence type of the Allele if different from the reference Allele
if ($allele->ref_allele_string eq $string){
#same allele as reference, there is no consequence, called SARA
$consequence_type->type('SARA');
push @out, $consequence_type;
next;
}
#calculate the consequence type of the Allele
my $ref_consequences = type_variation($transcript,"",$consequence_type);
if ($allele->start != $allele->end){
#do not calculate for indels effects of 2 or more in same codon
......
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