# HG changeset patch
# User iuc
# Date 1571217885 14400
# Node ID c23da6257d6a066a3ad619af0055baef3dda206e
# Parent 58b278def57eff40681cce10a3e511c84ab1e7c4
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/rgrnastar commit 2082c018009fa73c4afee8313febab13bb807ea8"
diff -r 58b278def57e -r c23da6257d6a rg_rnaStarSolo.xml
--- a/rg_rnaStarSolo.xml Fri Sep 06 11:10:22 2019 -0400
+++ b/rg_rnaStarSolo.xml Wed Oct 16 05:24:45 2019 -0400
@@ -2,7 +2,7 @@
mapping, demultiplexing and gene quantification for single cell RNA-seq
macros.xml
-
+ 1
@@ -14,13 +14,52 @@
STAR
@REFGENOMEHANDLING@
- ## cDNA sequence always goes first, then barcode
+ ## Check that the input pairs are of the same type
+ ## otherwise STARsolo will run for a long time and then error out.
+ ## We consume either repeats of two inputs R1 + R2
+ ## or a collection of paired reads.
+
+ #try
+ #set $last = None
+ #for $x in $input_types.input_repeats:
+ #if str($input_types.use) == "repeat":
+ #set $r1 = $x.input1
+ #set $r2 = $x.input2
+ #elif str($input_types.use) == "list_paired":
+ #set $r1 = $x.forward
+ #set $r2 = $x.reverse
+ #else
+ Wrong Type
+ #stop
+ #end if
+
+ #assert $r1.datatype == $r2.datatype
+
+ ## Test that all pairs are of the same type
+ #if $last:
+ #assert $last.datatype == $r1.datatype
+ #end if
+ #set $last = $r1
+ #end for
+ #except AssertionError
+ Input types are not the same!
+ #stop
+ #end try
+
+ ## cDNA sequence(s) [R2] always go first, then barcode(s) [R1]
+ ## see: Section 3.1 of STAR manual for multiple inputs, and Section 13 for STARsolo inputs
+ #if str($input_types.use) == "repeat":
+ #set $reads2 = ','.join([ '%s' % $x.input2 for $i,$x in enumerate($input_types.input_repeats)])
+ #set $reads1 = ','.join([ '%s' % $x.input1 for $i,$x in enumerate($input_types.input_repeats)])
+ #else if str($input_types.use) == "list_paired"
+ #set $reads2 = ','.join([ '%s' % $x.reverse for $i,$x in enumerate($input_types.input_repeats)])
+ #set $reads1 = ','.join([ '%s' % $x.forward for $i,$x in enumerate($input_types.input_repeats)])
+ #end if
+
--readFilesIn
- #set $reads2 = ','.join([ '%s' % $x.input2 for $i,$x in enumerate($input_repeats)])
- #set $reads1 = ','.join([ '%s' % $x.input1 for $i,$x in enumerate($input_repeats)])
$reads2 $reads1
- #if $input_repeats[0].input1.is_of_type('fastq.gz', 'fastqsanger.gz'):
+ #if $last.is_of_type('fastq.gz', 'fastqsanger.gz'):
@FASTQ_GZ_OPTION@
#end if
@@ -29,8 +68,8 @@
## 1 - check length of barcode, 0 - do not check
## Good for checking custom chemistries
- --soloBarcodeReadLength 1
--soloCBwhitelist '$soloCBwhitelist'
+ --soloBarcodeReadLength '$solo.soloBarcodeReadLength'
#if str($solo.params.chemistry) == "CR2":
--soloCBstart 1
@@ -54,10 +93,21 @@
--soloUMIdedup '$solo.soloUMIdedup'
]]>
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
@@ -91,6 +141,7 @@
+
@@ -112,10 +163,13 @@
-
-
-
-
+
+
+
+
+
+
+
@@ -153,10 +207,13 @@
-
-
-
-
+
+
+
+
+
+
+
@@ -188,18 +245,78 @@
-
-
-
-
-
-
-
-
-
-
-
-
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+
+