0
|
1
|
|
2 <!--
|
|
3 Copyright (C) 2011-2014 CRS4.
|
|
4
|
|
5 This file is part of Seal.
|
|
6
|
|
7 Seal is free software: you can redistribute it and/or modify it
|
|
8 under the terms of the GNU General Public License as published by the Free
|
|
9 Software Foundation, either version 3 of the License, or (at your option)
|
|
10 any later version.
|
|
11
|
|
12 Seal is distributed in the hope that it will be useful, but
|
|
13 WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY
|
|
14 or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
|
15 for more details.
|
|
16
|
|
17 You should have received a copy of the GNU General Public License along
|
|
18 with Seal. If not, see <http://www.gnu.org/licenses/>.
|
|
19 -->
|
|
20
|
|
21
|
|
22 <tool id="seal_seqal" name="Seqal" version="13986416aa79561bd0102cb7ccc1e0668ac9f0a4">
|
|
23 <description>Map reads on Hadoop</description>
|
|
24 <requirements>
|
|
25 <requirement type="package" version="13986416aa79561bd0102cb7ccc1e0668ac9f0a4">seal</requirement>
|
|
26 <requirement type="package" version="0.11">pydoop</requirement>
|
|
27 <requirement type="package" version="0.1.3">hadoop-galaxy</requirement>
|
|
28 </requirements>
|
|
29
|
|
30 <command>
|
|
31 hadoop_galaxy
|
|
32 --input $input_data
|
|
33 --output $output1
|
|
34 --executable seal
|
|
35 seqal
|
|
36 #if $align_only.value:
|
|
37 --align-only --num-reducers 0
|
|
38 #else
|
|
39 --num-reducers $align_only.num_reducers
|
|
40 #end if
|
|
41 --trimq $trimq
|
|
42 ${reference.fields.path}
|
|
43 </command>
|
|
44
|
|
45 <inputs>
|
|
46 <param name="input_data" type="data" format="pathset" label="Input data" />
|
|
47
|
|
48 <param name="reference" type="select" label="Select a built-in reference index archive">
|
|
49 <options from_data_table="seqal_indexes">
|
|
50 </options>
|
|
51 </param>
|
|
52
|
|
53 <param name="trimq" type="integer" min="0" value="0" label="trim quality, like BWA’s -q argument" />
|
|
54
|
|
55 <conditional name="align_only">
|
|
56 <param name="value" type="boolean" default="false" label="Align only (don't identify duplicates)" />
|
|
57 <when value="false">
|
|
58 <param name="num_reducers"
|
|
59 label="Number of reduce tasks"
|
|
60 type="integer"
|
|
61 value="90"
|
|
62 min="1"
|
|
63 />
|
|
64 </when>
|
|
65 </conditional>
|
|
66 </inputs>
|
|
67
|
|
68 <outputs>
|
|
69 <data name="output1" format="pathset" />
|
|
70 </outputs>
|
|
71
|
|
72 <stdio>
|
|
73 <exit_code range="1:" level="fatal" />
|
|
74 </stdio>
|
|
75
|
|
76 <help>
|
|
77 Seqal is a distributed short read mapping and duplicate removal tool. It
|
|
78 implements a distributed version of the BWA aligner, and adds a duplicate
|
|
79 read identification feature using the same criteria as the Picard
|
|
80 MarkDuplicates command. For a full description see the `manual
|
|
81 <http://biodoop-seal.sourceforge.net/seqal_index.html>`_.
|
|
82 </help>
|
|
83 </tool>
|