Commit 5b0e0ad2 authored by James Allen's avatar James Allen
Browse files

Specify 5' or 3' UTR in the GTF file

parent 78c9983c
......@@ -404,8 +404,8 @@ sub print_feature {
push @$utrs, @{$transcript->get_all_three_prime_UTRs()};
foreach my $utr (@{$utrs}) {
my $strand = $strand_conversion{$utr->strand()};
print $fh sprintf(qq{%s\t%s\tUTR\t%d\t%d\t.\t%s\t.\t},
$idstr, $transcript->source, ($utr->start()+$sliceoffset), ($utr->end+$sliceoffset), $strand);
print $fh sprintf(qq{%s\t%s\t%s\t%d\t%d\t.\t%s\t.\t},
$idstr, $transcript->source, $utr->type, ($utr->start()+$sliceoffset), ($utr->end+$sliceoffset), $strand);
$self->_print_attribs($gene, $biotype_display, $transcript, $transcript_biotype, 0, 'UTR', undef, undef, $has_selenocysteine);
print $fh "\n";
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