rsem v1.1.21

[rsem.git] / rsem-plot-transcript-wiggles

diff --git a/rsem-plot-transcript-wiggles b/rsem-plot-transcript-wiggles
index 5054f5bc0d4fa3256bc4ac4899162c4b4e570847..25e8b355dae5992cfd003bb3a659a0d28662f29a 100755 (executable)
--- a/rsem-plot-transcript-wiggles
+++ b/rsem-plot-transcript-wiggles
@@ -19,8 +19,12 @@ pod2usage(-msg => "Invalid number of arguments!", -exitval => 2, -verbose => 2)
 my ($fn, $dir, $suf) = fileparse($0);
 my $command = "";
 
+unless (-e "$ARGV[0].transcript.sorted.bam") {
+    $command = $dir."sam/samtools sort $ARGV[0].transcript.bam $ARGV[0].transcript.sorted";
+    &runCommand($command);
+}
 unless (-e "$ARGV[0].transcript.readdepth") {
-    $command = $dir."rsem-bam2readdepth $ARGV[0].transcript.sorted.bam > $ARGV[0].transcript.readdepth";
+    $command = $dir."rsem-bam2readdepth $ARGV[0].transcript.sorted.bam $ARGV[0].transcript.readdepth";
     &runCommand($command);
 }
 
@@ -34,7 +38,7 @@ if ($show_unique) {
        &runCommand($command);
     }
     unless (-e "$ARGV[0].uniq.transcript.readdepth") {
-       $command = $dir."rsem-bam2readdepth $ARGV[0].uniq.transcript.sorted.bam > $ARGV[0].uniq.transcript.readdepth";
+       $command = $dir."rsem-bam2readdepth $ARGV[0].uniq.transcript.sorted.bam $ARGV[0].uniq.transcript.readdepth";
        &runCommand($command);
     }
 }