3
|
1 #!/usr/bin/perl
|
|
2
|
|
3 #-------------------------------------------------------------------------------
|
|
4 # Add 'Parent' option in mRNA lines
|
|
5 #
|
|
6 # Pablo Cingolani
|
|
7 #-------------------------------------------------------------------------------
|
|
8
|
|
9 while( $l = <STDIN> ) {
|
|
10 chomp $l;
|
|
11 @t = split /\t/,$l;
|
|
12
|
|
13 $type = $t[2];
|
|
14
|
|
15 # Is it an mRNA?
|
|
16 if( $type eq 'mRNA' ) {
|
|
17 # Parse ID in options
|
|
18 $opts = $t[8];
|
|
19 if( $opts =~ /ID=(.*?)-\d+;/ ) {
|
|
20 $pid = $1;
|
|
21 $pid =~ tr/t/g/;
|
|
22 # Add 'Parent' option
|
|
23 $l .= ";Parent=$pid";
|
|
24 }
|
|
25 }
|
|
26
|
|
27 print "$l\n";
|
|
28 }
|