--- a/rg_rnaStar.xml	Wed Feb 14 09:03:31 2024 +0000
+++ b/rg_rnaStar.xml	Tue Aug 27 14:11:16 2024 +0000
@@ -47,9 +47,11 @@
 
         ## Two pass mode
         --twopassMode ${twopass.twopassMode} ${twopass.twopass_read_subset}
-        #for $sj_input in $twopass.sj_precalculated:
-            '$sj_input'
-        #end for
+        #if str($twopass.sj_precalculated).strip():
+            #for $sj_input in $twopass.sj_precalculated:
+                '$sj_input'
+            #end for
+        #end if
         #if str($twopass.twopassMode) != 'None':
             #if str($refGenomeSource.GTFconditional.GTFselect) == 'with-gtf':
                 ## need to check first if its a cached index or from history
@@ -215,7 +217,7 @@
             #end if
 
             ## Limits
-                @LIMITS@
+            @LIMITS@
         #else:
             ## Go with STAR's default algorithmic settings,
             ## but we need to provide a reasonable default
@@ -232,12 +234,16 @@
                 #end if
             #end if
         #end if
-
         --outBAMsortingThreadN \${GALAXY_SLOTS:-4}
         --outBAMsortingBinsN $perf.outBAMsortingBinsN
         --winAnchorMultimapNmax $perf.winAnchorMultimapNmax
         --limitBAMsortRAM \$((\${GALAXY_MEMORY_MB:-0}*1000000))
 
+        #if $oformat.wasp_conditional.waspOutputMode == "wasp_mode":
+            --waspOutputMode SAMtag
+            --varVCFfile '$oformat.wasp_conditional.varVCFfile'
+        #end if
+
         ## Handle chimeric options and output
         #if str($chimOutType):
             --chimOutType $chimOutType
@@ -408,6 +414,7 @@
             primary?"/> -->
             <param name="outSAMprimaryFlag" type="hidden" value="OneBestScore" />
             <expand macro="outSAMmapqUnique"/>
+            <expand macro="wasp"/>
         </section>
         <section name="filter" title="Output filter criteria" expanded="true">
             <param name="basic_filters" type="select" display="checkboxes" multiple="true" optional="true"
@@ -565,20 +572,24 @@
         </data>
         <expand macro="outWigOutputs"/>
     </outputs>
-
     <tests>
         <test expect_num_outputs="3">
             <conditional name="singlePaired">
-                <param name="sPaired" value="single" />
-                <param name="input1" value="tophat_in2.fastqsanger" ftype="fastqsanger" />
+                <param name="sPaired" value="paired" />
+                <param name="input1" value="pbmc_1k_v2_L001.R1.10k.fastq.gz" ftype="fastqsanger.gz" />
+                <param name="input2" value="pbmc_1k_v2_L001.R2.10k.fastq.gz" ftype="fastqsanger.gz" />
             </conditional>
             <conditional name="refGenomeSource">
                 <param name="geneSource" value="history" />
-                <param name="genomeFastaFiles" value="tophat_test.fa.gz" />
+                <param name="genomeFastaFiles" value="filtered3.Homo_sapiens.GRCh38.dna.chromosome.21.fa.gz" />
                 <param name="genomeSAindexNbases" value="5" />
             </conditional>
             <section name="oformat">
                 <param name="outSAMattributes" value="NH,HI,AS,nM,NM,MD,jM,jI,MC,ch" />
+                <conditional name="wasp_conditional">
+                    <param name="waspOutputMode" value="wasp_mode"/>
+                    <param name="varVCFfile" value="filtered3.vcf" ftype="vcf" />
+                </conditional>
             </section>
             <section name="algo">
                 <conditional name="params">
@@ -586,8 +597,7 @@
                 </conditional>
             </section>
             <output name="output_log" file="rnastar_test.log" compare="re_match_multiline" />
-            <output name="splice_junctions" file="rnastar_test_splicejunctions.bed"/>
-            <output name="mapped_reads" file="rnastar_test_mapped_reads.bam" compare="sim_size" delta="634" />
+            <output name="splice_junctions" file="rnastar_test_splicejunctions_wasp.bed"/>
         </test>
         <!-- test with cached genome index -->
         <test expect_num_outputs="3">