rsem v1.1.12 extract-transcript-to-gene-map-from-trinity changed

diff --git a/extract-transcript-to-gene-map-from-trinity b/extract-transcript-to-gene-map-from-trinity
index e04e99305aa67f8ed6b3076b9e435de9f45f033a..260cd834972f663f9df143f10a74e35518074028 100755 (executable)
--- a/extract-transcript-to-gene-map-from-trinity
+++ b/extract-transcript-to-gene-map-from-trinity
@@ -20,9 +20,9 @@ while (substr($tag, 0, 1) eq ">") {
     }
     if ($cnt == 0) { print "Warning: Fasta entry $tag has an empty sequence, it is omitted.\n"; }
     else {
-       my ($tid, $gid) = split(/ /, $tag);
-       my ($comp, $c, @tmp) = split(/_/, $gid);
-       $gid = join("_", $comp, $c);
+       my ($tid, @tmp) = split(/ /, $tag);
+       my ($comp, $c, @tmp) = split(/_/, $tag);
+       my $gid = join("_", $comp, $c);
        print OUTPUT "$gid\t$tid\n";
     }
     $tag = $line; chomp($tag);