|
2
|
1 <tool id="usearch_fastq_mergepairs" name="usearch fastq_mergepairs" version="0.0.1">
|
|
|
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'
|
|
16
|
9 #if $minovlen == "0"
|
|
13
|
10 -fastq_minovlen 1
|
|
|
11 #else
|
|
|
12 -fastq_minovlen $minovlen
|
|
7
|
13 #end if
|
|
|
14 -fastq_qmax $qmax
|
|
|
15 -fastqout '$output'
|
|
2
|
16 </command>
|
|
|
17 <inputs>
|
|
7
|
18 <param name="input_forward" type="data" format="fastq,fastqsanger,fastqcssanger" label="File with forward reads" />
|
|
|
19 <param name="input_reverse" type="data" format="fastq,fastqsanger,fastqcssanger" label="File with reverse reads" />
|
|
11
|
20 <param name="minovlen" type="integer" value="0" label="Minimum length of the overlap" help="'0' means no minimum." />
|
|
10
|
21 <param name="qmax" type="integer" value="41" label="Maximum Q score for input files" />
|
|
2
|
22 </inputs>
|
|
|
23 <outputs>
|
|
|
24 <data name='output' format='fastq' label="Merge result" />
|
|
|
25 </outputs>
|
|
|
26 <tests>
|
|
|
27 <test>
|
|
|
28 <param name="input_forward" value="fastq_mergepairs_input1.fq" ftype="fastqsanger" />
|
|
|
29 <param name="input_reverse" value="fastq_mergepairs_input2.fq" ftype="fastqsanger" />
|
|
5
|
30 <param name="qmax" value="65" />
|
|
2
|
31 <output name="output" file="fastq_mergepairs_output.fq" />
|
|
|
32 </test>
|
|
|
33 </tests>
|
|
|
34 <help>
|
|
|
35 **What it Does**
|
|
|
36
|
|
|
37 Performs merging of paired reads.
|
|
|
38
|
|
|
39 The FASTQ filename for the forward reads is specified by the -fastq_mergepairs option, and the reverse read filename is specified by the -reverse option. Output files are specified by -fastqout (for FASTQ) and / or -fastaout (for FASTA).
|
|
|
40
|
|
|
41 Forward and reverse must be in 1:1 correspondence and must appear in the same order in both files. The labels for the forward and reverse read in a given pair must be identical except for a single position where a '1' appears in the forward read label and a '2' appears in the reverse read label.
|
|
|
42
|
|
|
43 -----
|
|
|
44
|
|
|
45 **Input formats**
|
|
|
46
|
|
|
47 Forward read::
|
|
|
48
|
|
|
49 @IRIS:7:1:29:952#0/1
|
|
|
50 TGAGAAGCAAGAAGAAGGTTGGTTAGTGTTTTGGAG
|
|
|
51 +IRIS:7:1:29:952#0/1
|
|
|
52 aaabaaaaaaaaaaa`aaY`aa^aaa^a_a_`aa``
|
|
|
53
|
|
|
54 Reverse read::
|
|
|
55
|
|
|
56 @IRIS:7:1:29:952#0/2
|
|
|
57 GACTCCAAAACACTAACCAACCTTCTTCTTGCTTCT
|
|
|
58 +IRIS:7:1:29:952#0/2
|
|
|
59 aaaabaaaabaaaabbaaaa````__`__^__``__
|
|
|
60
|
|
|
61 -----
|
|
|
62
|
|
|
63 **Output**
|
|
|
64
|
|
|
65 A multiple-fastq file, for example::
|
|
|
66
|
|
|
67 @IRIS:7:1:29:952#0/1
|
|
|
68 TGAGAAGCAAGAAGAAGGTTGGTTAGTGTTTTGGAGTC
|
|
|
69 +
|
|
|
70 aaJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJJaa
|
|
|
71
|
|
|
72 ------
|
|
|
73
|
|
|
74 **Author**
|
|
|
75
|
|
|
76 Robert C. Edgar (robert@drive5.com)
|
|
|
77
|
|
|
78 **Manual**
|
|
|
79
|
|
|
80 http://drive5.com/usearch/manual/fastq_mergepairs.html
|
|
|
81
|
|
|
82 **Citation**
|
|
|
83
|
|
|
84 Please cite one of these papers if you use USEARCH in published work.
|
|
|
85
|
|
|
86 Edgar,RC (2010) Search and clustering orders of magnitude faster than BLAST, Bioinformatics 26(19), 2460-2461.
|
|
|
87 doi: 10.1093/bioinformatics/btq461
|
|
|
88 </help>
|
|
|
89 </tool>
|