|
32
|
1 <tool id="usearch_fastq_mergepairs" name="usearch fastq_mergepairs" version="0.0.2">
|
|
2
|
2 <description>merging of paired reads</description>
|
|
|
3 <version_command>usearch -version</version_command>
|
|
12
|
4 <command interpreter='bash'>usearch_wrapper.sh
|
|
|
5 usearch
|
|
|
6 -quiet
|
|
7
|
7 -fastq_mergepairs '$input_forward'
|
|
|
8 -reverse '$input_reverse'
|
|
17
|
9 #if $minovlen.value != 0
|
|
13
|
10 -fastq_minovlen $minovlen
|
|
7
|
11 #end if
|
|
22
|
12 #if $minmergelen.value != 0
|
|
|
13 -fastq_minmergelen $minmergelen
|
|
|
14 #end if
|
|
24
|
15 #if $maxmergelen.value != 0
|
|
|
16 -fastq_maxmergelen $maxmergelen
|
|
|
17 #end if
|
|
|
18 #if $maxdiffs.value != 0
|
|
|
19 -fastq_maxdiffs $maxdiffs
|
|
|
20 #end if
|
|
|
21 #if $truncqual.value != 0
|
|
|
22 -fastq_truncqual $truncqual
|
|
|
23 #end if
|
|
|
24 #if $minlen.value != 0
|
|
|
25 -fastq_minlen $minlen
|
|
|
26 #end if
|
|
25
|
27 $allowmergestagger
|
|
|
28 -fastq_ascii $ascii
|
|
|
29 -fastq_qmin $qmin
|
|
7
|
30 -fastq_qmax $qmax
|
|
25
|
31 -fastq_qmaxout $qmaxout
|
|
34
|
32 #if $out_format.value == "fastq"
|
|
26
|
33 -fastqout $output
|
|
|
34 #else
|
|
|
35 -fastaout $output
|
|
|
36 #end if
|
|
2
|
37 </command>
|
|
|
38 <inputs>
|
|
33
|
39 <param name="input_forward" type="data" format="fastq,fastqsanger,fastqcssanger" label="1. File with forward reads" />
|
|
|
40 <param name="input_reverse" type="data" format="fastq,fastqsanger,fastqcssanger" label="2. File with reverse reads" />
|
|
|
41 <param name="minovlen" type="integer" value="0" label="3. Minimum length of the overlap" help="'0' means no minimum." />
|
|
|
42 <param name="minmergelen" type="integer" value="0" label="4. Minimum length of the merged read" help="'0' means no minimum." />
|
|
|
43 <param name="maxmergelen" type="integer" value="0" label="5. Maximum length of the merged read" help="'0' means no maximum." />
|
|
|
44 <param name="maxdiffs" type="integer" value="0" label="6. Maximum number of mismatches allowed in the overlap region" help="'0' means any number of mismatches allowed." />
|
|
|
45 <param name="truncqual" type="integer" value="0" label="7. Truncate the forward and reverse reads at the first Q that is equal or less than this value, if present"
|
|
|
46 help="'0' means no quality truncation. This truncation is performed before aligning the pair. With Illumina paired reads, it is recommended to set this to 2 or higher, as low-quality tails will otherwise often cause alignment of the pair to fail." />
|
|
|
47 <param name="minlen" type="integer" value="0" label="8. Minimum length of the forward and reverse read, after truncating per option 7 if applicable" help="'0' means no minimum." />
|
|
27
|
48 <param name="allowmergestagger" type="boolean" truevalue="-fastq_allowmergestagger" falsevalue="" checked="false" label="9. Allow merge of a pair where the alignment is staggered" help="By default, pairs with staggered alignments are discarded." />
|
|
33
|
49 <param name="ascii" type="integer" value="33" label="10. ASCII_BASE constant" help="See http://drive5.com/usearch/manual/fastq_params.html" />
|
|
|
50 <param name="qmin" type="integer" value="0" label="11. Minimum Q score" />
|
|
|
51 <param name="qmax" type="integer" value="41" label="12. Maximum Q score for input files" />
|
|
|
52 <param name="qmaxout" type="integer" value="41" label="13. Maximum Q score for output files" />
|
|
|
53 <param name="out_format" type="select" label="Output format">
|
|
|
54 <option value="fastq">FASTQ</option>
|
|
|
55 <option value="fasta">FASTA</option>
|
|
|
56 </param>
|
|
17
|
57 </inputs>
|
|
2
|
58 <outputs>
|
|
33
|
59 <data format="fastq" name="output" label="Merge output">
|
|
|
60 <change_format>
|
|
|
61 <when input="out_format" value="fasta" format="fasta" />
|
|
|
62 </change_format>
|
|
|
63 </data>
|
|
2
|
64 </outputs>
|
|
|
65 <tests>
|
|
|
66 <test>
|
|
33
|
67 <param name="input_forward" value="fastq_mergepairs_input1.fq" ftype="fastq" />
|
|
|
68 <param name="input_reverse" value="fastq_mergepairs_input2.fq" ftype="fastq" />
|
|
5
|
69 <param name="qmax" value="65" />
|
|
33
|
70 <output name="output" file="fastq_mergepairs_output1.fastq" />
|
|
2
|
71 </test>
|
|
33
|
72 <test>
|
|
|
73 <param name="input_forward" value="fastq_mergepairs_input1.fq" ftype="fastq" />
|
|
|
74 <param name="input_reverse" value="fastq_mergepairs_input2.fq" ftype="fastq" />
|
|
|
75 <param name="minovlen" value="30" />
|
|
|
76 <param name="qmax" value="65" />
|
|
|
77 <param name="out_format" value="fasta" />
|
|
|
78 <output name="output" file="fastq_mergepairs_output2.fasta" />
|
|
|
79 </test>
|
|
|
80 </tests>
|
|
2
|
81 <help>
|
|
|
82 **What it Does**
|
|
|
83
|
|
|
84 Performs merging of paired reads.
|
|
|
85
|
|
31
|
86 Forward and reverse must be in 1:1 correspondence and must appear in the same
|
|
|
87 order in both files. The labels for the forward and reverse read in a given
|
|
|
88 pair must be identical except for a single position where a '1' appears in the
|
|
|
89 forward read label and a '2' appears in the reverse read label.
|
|
2
|
90
|
|
|
91 -----
|
|
|
92
|
|
|
93 **Input formats**
|
|
|
94
|
|
33
|
95 Forward reads::
|
|
2
|
96
|
|
|
97 @IRIS:7:1:29:952#0/1
|
|
|
98 TGAGAAGCAAGAAGAAGGTTGGTTAGTGTTTTGGAG
|
|
|
99 +IRIS:7:1:29:952#0/1
|
|
|
100 aaabaaaaaaaaaaa`aaY`aa^aaa^a_a_`aa``
|
|
|
101
|
|
33
|
102 Reverse reads::
|
|
2
|
103
|
|
|
104 @IRIS:7:1:29:952#0/2
|
|
|
105 GACTCCAAAACACTAACCAACCTTCTTCTTGCTTCT
|
|
|
106 +IRIS:7:1:29:952#0/2
|
|
|
107 aaaabaaaabaaaabbaaaa````__`__^__``__
|
|
|
108
|
|
|
109 -----
|
|
|
110
|
|
|
111 **Output**
|
|
|
112
|
|
|
113 A multiple-fastq file, for example::
|
|
|
114
|
|
|
115 @IRIS:7:1:29:952#0/1
|
|
|
116 TGAGAAGCAAGAAGAAGGTTGGTTAGTGTTTTGGAGTC
|
|
|
117 +
|
|
|
118 aaJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJaa
|
|
|
119
|
|
|
120 ------
|
|
|
121
|
|
|
122 **Manual**
|
|
|
123
|
|
33
|
124 * USEARCH fastq_mergepairs options: http://drive5.com/usearch/manual/fastq_mergepairs.html
|
|
30
|
125 * FASTQ format options: http://drive5.com/usearch/manual/fastq_params.html
|
|
2
|
126
|
|
|
127 **Citation**
|
|
|
128
|
|
|
129 Please cite one of these papers if you use USEARCH in published work.
|
|
|
130
|
|
|
131 Edgar,RC (2010) Search and clustering orders of magnitude faster than BLAST, Bioinformatics 26(19), 2460-2461.
|
|
|
132 doi: 10.1093/bioinformatics/btq461
|
|
|
133 </help>
|
|
|
134 </tool>
|
