comparison dada2_mergePairs.xml @ 0:805daed21ef3 draft

"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/dada2 commit f8b6b6e72914ad6bcca8423dfa03f59bde80992e"
author iuc
date Fri, 08 Nov 2019 18:51:19 -0500
parents
children b1e55c248f17
comparison
equal deleted inserted replaced
-1:000000000000 0:805daed21ef3
1 <tool id="dada2_mergePairs" name="dada2: mergePairs" version="@DADA2_VERSION@+galaxy@WRAPPER_VERSION@" profile="19.09">
2 <description>Merge denoised forward and reverse reads</description>
3 <macros>
4 <import>macros.xml</import>
5 </macros>
6 <expand macro="requirements"/>
7 <expand macro="stdio"/>
8 <expand macro="version_command"/>
9 <command detect_errors="exit_code"><![CDATA[
10 Rscript '$dada2_script'
11 ]]></command>
12 <configfiles>
13 <configfile name="dada2_script"><![CDATA[
14 library(dada2, quietly=T)
15
16 dadaF <- readRDS('$dadaF')
17 dadaR <- readRDS('$dadaR')
18
19 merged <- mergePairs(dadaF, '$derepF', dadaR, '$derepR', minOverlap = $minOverlap, maxMismatch = $maxMismatch,
20 ## currently not accessible by the wrapper returnRejects = FALSE, propagateCol = character(0),
21 justConcatenate = $justConcatenate, trimOverhang = $trimOverhang)
22
23 saveRDS(merged, file='$merged')
24 ##For display in Galaxy this could be transposed (..but the user could also do it herself)
25 ##write.table(merged, "$merged", quote=F, sep="\t", row.names = F, col.names = T)
26 ]]></configfile>
27 </configfiles>
28 <inputs>
29 <param argument="dadaF" type="data" format="dada2_dada" label="Dada results for forward reads"/>
30 <param argument="derepF" type="data" format="fastq,fastq.gz" label="Forward reads" help="despite the parameter name the sequences don't need to be dereplicated "/>
31 <param argument="dadaR" type="data" format="dada2_dada" label="Dada results for reverse reads"/>
32 <param argument="derepR" type="data" format="fastq,fastq.gz" label="Reverse reads" help="despite the parameter name the sequences don't need to be dereplicated "/>
33 <param argument="minOverlap" type="integer" value="12" min="0" label="Minimum length of the overlap"/>
34 <param argument="maxMismatch" type="integer" value="0" min="0" label="Maximum number of mismatches"/>
35 <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"/>
36 <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"/>
37 </inputs>
38 <outputs>
39 <data name="merged" format="dada2_mergepairs" label="${tool.name} on ${on_string}"/>
40 </outputs>
41 <tests>
42 <test>
43 <param name="dadaF" ftype="dada2_dada" value="dada_F3D0_R1.Rdata"/>
44 <param name="derepF" ftype="fastq.gz" value="filterAndTrim_F3D0_R1.fq.gz"/>
45 <param name="dadaR" ftype="dada2_dada" value="dada_F3D0_R2.Rdata"/>
46 <param name="derepR" ftype="fastq.gz" value="filterAndTrim_F3D0_R2.fq.gz"/>
47 <output name="merged" value="mergePairs_F3D0.Rdata" ftype="dada2_mergepairs" />
48 </test>
49 <!-- test non-default options -->
50 <test>
51 <param name="dadaF" ftype="dada2_dada" value="dada_F3D0_R1.Rdata"/>
52 <param name="derepF" ftype="fastq.gz" value="filterAndTrim_F3D0_R1.fq.gz"/>
53 <param name="dadaR" ftype="dada2_dada" value="dada_F3D0_R2.Rdata"/>
54 <param name="derepR" ftype="fastq.gz" value="filterAndTrim_F3D0_R2.fq.gz"/>
55 <output name="merged" value="mergePairs_F3D0_nondefault.Rdata" ftype="dada2_mergepairs" />
56 <param name="minOverlap" value="8" />
57 <param name="maxMismatch" value="1"/>
58 <param name="justConcatenate" value="TRUE" />
59 <param name="trimOverhang" value="TRUE" />
60 </test>
61 </tests>
62 <help><![CDATA[
63 Description
64 ...........
65
66 This function attempts to merge each denoised pair of forward and reverse reads, rejecting any pairs which do not sufficiently overlap (at least 12bp by default) or which contain too many (>0 by default) mismatches in the overlap region.
67
68 Usage
69 .....
70
71 **Input**
72
73 - reads and learned error rates of the forward reads
74 - reads and learned error rates of the reverse reads
75
76 **Output**
77
78 - a data set of type dada2_mergepairs (which is a RData file containing the output of dada2's mergePairs function).
79
80 Details
81 .......
82
83 @HELP_OVERVIEW@
84
85
86 ]]></help>
87 <expand macro="citations"/>
88 </tool>