Merge pull request #87 from drpatelh/rsem_index

Use STAR to build initial index for RSEM

software/rsem/preparereference/main.nf

@@ -24,15 +24,41 @@ process RSEM_PREPAREREFERENCE {
 
     script:
     def software = getSoftwareName(task.process)
-    """
+    def args     = options.args.tokenize()
-    mkdir rsem
+    if (args.contains('--star')) {
-    rsem-prepare-reference \\
+        args.removeIf { it.contains('--star') }
-        --gtf $gtf \\
+        def memory = task.memory ? "--limitGenomeGenerateRAM ${task.memory.toBytes() - 100000000}" : ''
-        --num-threads $task.cpus \\
+        """
-        $options.args \\
+        mkdir rsem
-        $fasta \\
+        STAR \\
-        rsem/genome
+            --runMode genomeGenerate \\
+            --genomeDir rsem/ \\
+            --genomeFastaFiles $fasta \\
+            --sjdbGTFfile $gtf \\
+            --runThreadN $task.cpus \\
+            $memory \\
+            $options.args2
+        
+        rsem-prepare-reference \\
+            --gtf $gtf \\
+            --num-threads $task.cpus \\
+            ${args.join(' ')} \\
+            $fasta \\
+            rsem/genome
 
-    rsem-calculate-expression --version | sed -e "s/Current version: RSEM v//g" > ${software}.version.txt
+        rsem-calculate-expression --version | sed -e "s/Current version: RSEM v//g" > ${software}.version.txt
-    """
+        """
+    } else {
+        """
+        mkdir rsem
+        rsem-prepare-reference \\
+            --gtf $gtf \\
+            --num-threads $task.cpus \\
+            $options.args \\
+            $fasta \\
+            rsem/genome
+
+        rsem-calculate-expression --version | sed -e "s/Current version: RSEM v//g" > ${software}.version.txt
+        """
+    }
 }