Skip to content
Snippets Groups Projects
Commit c1a07cb6 authored by Magali Ruffier's avatar Magali Ruffier
Browse files

using biotype group directly rather than list of biotypes belonging to said group

parent 1d969007
No related branches found
No related tags found
No related merge requests found
......@@ -172,14 +172,11 @@ sub exon {
$difference = diff(\$old_seq, \$projected_seq);
my $biotype_group = $biotype_manager->fetch_biotype($parent_gene)->biotype_group;
my $group_list = $biotype_manager->group_members($biotype_group);
foreach my $group (@$group_list) {
if ($group eq 'coding') {
$state = '!!';
last;
}
if ($biotype_group eq 'coding') {
$state = '!!';
} else {
$state = '%%';
}
$state = '%%' unless $state; # Middle badness
}
$location = sprintf('%d : %d', $projected_exon->start(), $projected_exon->end());
......@@ -240,17 +237,11 @@ sub transcript {
#Check if it was protein coding
my $biotype_group = $biotype_manager->fetch_biotype($projected_transcript)->biotype_group;
my $group_list = $biotype_manager->group_members($biotype_group);
my $is_pc = 0;
foreach my $group (@$group_list) {
if ($group eq 'coding') {
$is_pc = 1;
last;
}
}
my $is_pc = 0;
#Now check for protein sequence mis-match
if ($is_pc) {
if ($biotype_group eq 'coding') {
$is_pc = 1;
my $old_seq = $old_transcript->translate()->seq();
my $new_seq = $projected_transcript->translate()->seq();
$total_length = length($old_seq);
......
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