Commit 3d896a5a authored by James Allen's avatar James Allen
Browse files

For trans-spliced transcripts, augment the CDS IDs with the location, to...

For trans-spliced transcripts, augment the CDS IDs with the location, to prevent multi-row features having different seq_region names and/or strands
parent 28154de0
......@@ -206,6 +206,11 @@ sub print_feature {
$value = 'gene:' . $value;
} elsif ($feature->isa('Bio::EnsEMBL::Exon')) {
$key = 'Name';
} elsif ($feature->isa('Bio::EnsEMBL::CDS')) {
my $trans_spliced = $feature->transcript->get_all_Attributes('trans_spliced');
if (scalar(@$trans_spliced)) {
$value = $so_term . ':' . join('_', $value, $feature->seq_region_name, $feature->seq_region_strand);
}
} else {
$value = $so_term . ':' . $value;
}
......
Markdown is supported
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