--- a/dada2_mergePairs.xml	Thu May 23 08:54:18 2024 +0000
+++ b/dada2_mergePairs.xml	Fri May 24 11:45:13 2024 +0000
@@ -20,10 +20,12 @@
 merged <- mergePairs(dadaF, '$derepF', dadaR, '$derepR', minOverlap = $minOverlap, maxMismatch = $maxMismatch,
 ## currently not accessible by the wrapper    returnRejects = FALSE, propagateCol = character(0),
     justConcatenate = $justConcatenate, trimOverhang = $trimOverhang)
-
 saveRDS(merged, file='$merged')
-##For display in Galaxy this could be transposed (..but the user could also do it herself)
-##write.table(merged, "$merged", quote=F, sep="\t", row.names = F, col.names = T)
+#if $output_details
+    ##For display in Galaxy this could be transposed (..but the user could also do it herself)
+    merged <- apply(merged, 2, as.character)
+    write.table(merged, "$merged_details", quote=F, sep="\t", row.names = F, col.names = T)
+#end if
     ]]></configfile>
     </configfiles>
     <inputs>
@@ -35,12 +37,16 @@
         <param argument="maxMismatch" type="integer" value="0" min="0" label="Maximum number of mismatches"/>
         <param argument="justConcatenate" type="boolean" checked="false" truevalue="TRUE" falsevalue="FALSE" label="Concatenated rather than merge" help="The forward and reverse-complemented reverse read are concatenated rather than merged, with a 10 Ns spacer inserted between them"/>
         <param argument="trimOverhang" type="boolean" checked="false" truevalue="TRUE" falsevalue="FALSE" label="Trim overhangs" help="Overhangs are when the reverse read extends past the start of the forward read, and vice-versa, as can happen when reads are longer than the amplicon and read into the other-direction primer region"/>
+        <param name="output_details" type="boolean" label="Output detailed table" help="Contains detailed information on the merged sequences, e.g. number of matched/mismatches"/>
     </inputs>
     <outputs>
         <data name="merged" format="dada2_mergepairs" label="${tool.name} on ${on_string}"/>
+        <data name="merged_details" format="tabular" label="${tool.name} on ${on_string}: details">
+            <filter>output_details</filter>
+        </data>
     </outputs>
     <tests>
-        <test>
+        <test expect_num_outputs="1">
             <param name="dadaF" ftype="dada2_dada" value="dada_F3D0_S188_L001_R1.Rdata"/>
             <param name="derepF" ftype="fastq.gz" value="filterAndTrim_F3D0_R1.fq.gz"/>
             <param name="dadaR" ftype="dada2_dada" value="dada_F3D0_S188_L001_R2.Rdata"/>
@@ -48,7 +54,7 @@
             <output name="merged" value="mergePairs_F3D0_S188_L001.Rdata" ftype="dada2_mergepairs" compare="sim_size" delta="400" />
         </test>
         <!-- test non-default options -->
-        <test>
+        <test expect_num_outputs="2">
             <param name="dadaF" ftype="dada2_dada" value="dada_F3D0_S188_L001_R1.Rdata"/>
             <param name="derepF" ftype="fastq.gz" value="filterAndTrim_F3D0_R1.fq.gz"/>
             <param name="dadaR" ftype="dada2_dada" value="dada_F3D0_S188_L001_R2.Rdata"/>
@@ -58,6 +64,13 @@
             <param name="maxMismatch" value="1"/>
             <param name="justConcatenate" value="TRUE" />
             <param name="trimOverhang" value="TRUE" />
+            <param name="output_details" value="true"/>
+            <output name="merged_details">
+                <assert_contents>
+                    <has_n_lines n="50"/>
+                    <has_n_columns n="9"/>
+                </assert_contents>
+            </output>
         </test>
     </tests>
     <help><![CDATA[