annotate bwa.xml @ 13:53646aaaafef draft

planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
author devteam
date Mon, 30 Jan 2017 08:49:26 -0500
parents bd3a1e0de84c
children 051eba708f43
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
1 <?xml version="1.0"?>
11
546ada4a9f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 610045611b00099e0a24183c8cf33aebfa9635cf
devteam
parents: 10
diff changeset
2 <tool id="bwa" name="Map with BWA" version="@VERSION@.1">
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
3 <description>- map short reads (&lt; 100 bp) against reference genome</description>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
4 <macros>
5
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
5 <import>read_group_macros.xml</import>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
6 <import>bwa_macros.xml</import>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
7 <token name="@command_options@">
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
8 #if str( $analysis_type.analysis_type_selector ) == "full":
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
9 -n ${analysis_type.n}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
10 -o ${analysis_type.o}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
11 -e ${analysis_type.e}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
12 -i ${analysis_type.i}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
13 -d ${analysis_type.d}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
14 -l ${analysis_type.l}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
15 -k ${analysis_type.k}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
16 -m ${analysis_type.m}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
17 -M ${analysis_type.M}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
18 -O ${analysis_type.O}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
19 -E ${analysis_type.E}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
20 -R ${analysis_type.R}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
21 -q ${analysis_type.q}
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
22
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
23 #if str( $analysis_type.B ):
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
24 -B ${analysis_type.B}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
25 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
26
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
27 #if str( $analysis_type.L ):
4
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
28 -L ${analysis_type.L}
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
29 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
30 #end if
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
31 </token>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
32 <token name="@read_group_options@">
5
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
33 #if $use_rg:
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
34 @set_rg_string@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
35 -r '$rg_string'
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
36 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
37 </token>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
38
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
39 <xml name="advanced_pe_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
40 <param name="adv_pe_options_selector" type="select" label="Set advanced paired end options?" help="Provides additional controls">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
41 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
42 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
43 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
44 <when value="set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
45 <param name="a" type="integer" value="500" label="Maximum insert size for a read pair to be considered being mapped properly." help="sampe -a; This option is only used when there are not enough good alignment to infer the distribution of insert sizes; default=500"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
46 <param name="o" type="integer" value="100000" label="Maximum occurrences of a read for pairing. A read with more occurrences will be treated as a single-end read." help="sampe -o; Reducing this parameter helps faster pairing; default=100000"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
47 <param name="n" type="integer" value="3" label="Maximum number of alignments to output in the XA tag for reads paired properly." help="sampe -n; If a read has more than this many hits, the XA tag will not be written; default=3"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
48 <param name="N" type="integer" value="10" label="Maximum number of alignments to output in the XA tag for disconcordant read pairs (excluding singletons)." help="sampe -N; If a read has more than this many hits, the XA tag will not be written; default=10"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
49 <param name="c" type="float" value="0.00005" label="Prior of chimeric rate (lower bound)" help="sampe -c"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
50 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
51 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
52 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
53 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
54 </xml>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
55 <xml name="advanced_se_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
56 <param name="adv_se_options_selector" type="select" label="Set advanced single end options?" help="Provides additional controls">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
57 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
58 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
59 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
60 <when value="set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
61 <param name="n" type="integer" value="3" label="Maximum number of alignments to output in the XA tag." help="-n; If a read has more than this many hits, the XA tag will not be written; default=3"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
62 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
63 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
64 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
65 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
66 </xml>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
67 </macros>
7
d8c9597bfb09 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit ddb8bdb9d62399f086b06b3469450d0aad2113bd
devteam
parents: 6
diff changeset
68 <expand macro="requirements" />
6
09a7281d24c5 planemo upload commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents: 5
diff changeset
69 <expand macro="stdio" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
70 <command>
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
71 <![CDATA[
11
546ada4a9f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 610045611b00099e0a24183c8cf33aebfa9635cf
devteam
parents: 10
diff changeset
72 @set_reference_fasta_filename@
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
73
5
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
74 ## setup vars for rg handling...
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
75 @define_read_group_helpers@
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
76 #if str( $input_type.input_type_selector ) == "paired":
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
77 #set $rg_auto_name = $read_group_name_default($input_type.fastq_input1, $input_type.fastq_input2)
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
78 #elif str( $input_type.input_type_selector ) in ["single_bam", "paired_bam"]:
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
79 #set $rg_auto_name = $read_group_name_default($input_type.bam_input)
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
80 #else
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
81 #set $rg_auto_name = $read_group_name_default($input_type.fastq_input1)
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
82 #end if
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
83 @set_use_rg_var@
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
84 @set_read_group_vars@
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
85
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
86 ## Begin bwa command line
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
87
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
88 ####### Fastq paired
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
89
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
90 #if str( $input_type.input_type_selector ) == "paired" or str( $input_type.input_type_selector ) == "paired_collection":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
91 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
92 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
93
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
94 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
95
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
96 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
97
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
98 #if str( $input_type.input_type_selector ) == "paired_collection":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
99 "${input_type.fastq_input1.forward}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
100 #else
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
101 "${input_type.fastq_input1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
102 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
103
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
104 > first.sai &&
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
105
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
106 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
107 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
108
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
109 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
110
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
111 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
112
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
113 #if str( $input_type.input_type_selector ) == "paired_collection":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
114 "${input_type.fastq_input1.reverse}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
115 #else
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
116 "${input_type.fastq_input2}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
117 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
118
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
119 > second.sai &&
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
120
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
121 bwa sampe
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
122
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
123 #if str( $input_type.adv_pe_options.adv_pe_options_selector) == "True":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
124 -a ${$input_type.adv_pe_options.a}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
125 -o ${$input_type.adv_pe_options.o}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
126 -n ${$input_type.adv_pe_options.n}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
127 -N ${$input_type.adv_pe_options.N}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
128 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
129
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
130 @read_group_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
131
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
132 #if str( $input_type.input_type_selector ) == "paired_collection":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
133 "${reference_fasta_filename}" first.sai second.sai "${input_type.fastq_input1.forward}" "${input_type.fastq_input1.reverse}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
134 #else:
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
135 "${reference_fasta_filename}" first.sai second.sai "${input_type.fastq_input1}" "${input_type.fastq_input2}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
136 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
137
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
138 ####### Fastq single
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
139
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
140 #elif str( $input_type.input_type_selector ) == "single":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
141 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
142 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
143
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
144 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
145
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
146 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
147 "${input_type.fastq_input1}"
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
148 > first.sai &&
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
149
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
150 bwa samse
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
151
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
152 #if str( $input_type.adv_se_options.adv_se_options_selector) == "True":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
153 -n ${$input_type.adv_se_options.n}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
154 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
155
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
156 @read_group_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
157
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
158 "${reference_fasta_filename}" first.sai "${input_type.fastq_input1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
159
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
160 ####### BAM paired
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
161
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
162 #elif str( $input_type.input_type_selector ) == "paired_bam":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
163 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
164 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
165 -b
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
166 -1
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
167
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
168 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
169
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
170 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
171 "${input_type.bam_input}"
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
172 > first.sai &&
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
173
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
174 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
175 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
176 -b
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
177 -2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
178 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
179 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
180 "${input_type.bam_input}"
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
181 > second.sai &&
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
182
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
183 bwa sampe
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
184
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
185 #if str( $input_type.adv_bam_pe_options.adv_pe_options_selector) == "True":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
186 -a ${$input_type.adv_bam_pe_options.a}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
187 -o ${$input_type.adv_bam_pe_options.o}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
188 -n ${$input_type.adv_bam_pe_options.n}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
189 -N ${$input_type.adv_bam_pe_options.N}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
190 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
191
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
192 @read_group_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
193
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
194 "${reference_fasta_filename}" first.sai second.sai "${input_type.bam_input}" "${input_type.bam_input}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
195
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
196 ####### Fastq single ------------ to do next
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
197
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
198 #elif str( $input_type.input_type_selector ) == "single_bam":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
199 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
200 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
201 -b
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
202 -0
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
203
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
204 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
205
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
206 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
207 "${input_type.bam_input}"
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
208 > first.sai &&
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
209
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
210 bwa samse
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
211
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
212 #if str( $input_type.adv_bam_se_options.adv_se_options_selector) == "True":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
213 -n ${$input_type.adv_bam_se_options.n}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
214 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
215
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
216 @read_group_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
217
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
218 "${reference_fasta_filename}" first.sai "${input_type.bam_input}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
219 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
220
12
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
221 | samtools sort -O bam -o '$bam_output'
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
222 ]]>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
223 </command>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
224
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
225 <inputs>
11
546ada4a9f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 610045611b00099e0a24183c8cf33aebfa9635cf
devteam
parents: 10
diff changeset
226 <expand macro="reference_source_conditional" />
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
227 <conditional name="input_type">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
228 <param name="input_type_selector" type="select" label="Select input type" help="Select between fastq and bam datasets and between paired and single end data">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
229 <option value="paired">Paired fastq</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
230 <option value="paired_collection">Paired fastq collection</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
231 <option value="single">Single fastq</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
232 <option value="paired_bam">Paired BAM</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
233 <option value="single_bam">Single BAM</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
234 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
235 <when value="paired">
13
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
236 <param name="fastq_input1" type="data" format="fastqsanger,fastqsanger.gz" label="Select first set of reads" help="Specify dataset with forward reads"/>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
237 <param name="fastq_input2" type="data" format="fastqsanger,fastqsanger.gz" label="Select second set of reads" help="Specify dataset with reverse reads"/>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
238 <conditional name="adv_pe_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
239
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
240 <expand macro="advanced_pe_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
241
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
242 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
243 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
244
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
245 <when value="paired_collection">
13
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
246 <param name="fastq_input1" format="fastqsanger,fastqsanger.gz" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
247 <conditional name="adv_pe_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
248
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
249 <expand macro="advanced_pe_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
250
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
251 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
252 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
253
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
254 <when value="single">
13
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
255 <param name="fastq_input1" type="data" format="fastqsanger,fastqsanger.gz" label="Select fastq dataset" help="Specify dataset with single reads"/>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
256 <conditional name="adv_se_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
257
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
258 <expand macro="advanced_se_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
259
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
260 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
261 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
262
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
263 <!-- the difference between single and paired bams is in the <command> tag portion and realated to -0, -1, and -2 options -->
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
264
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
265 <when value="paired_bam">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
266 <param name="bam_input" type="data" format="bam" label="Select BAM dataset" help="Specify BAM dataset with paired reads"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
267 <conditional name="adv_bam_pe_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
268
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
269 <expand macro="advanced_pe_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
270
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
271 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
272 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
273
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
274 <when value="single_bam">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
275 <param name="bam_input" type="data" format="bam" label="Select BAM dataset" help="Specify BAM dataset with single reads"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
276 <conditional name="adv_bam_se_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
277
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
278 <expand macro="advanced_se_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
279
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
280 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
281 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
282
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
283 </conditional>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
284
5
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
285 <expand macro="read_group_conditional" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
286
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
287 <conditional name="analysis_type">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
288 <param name="analysis_type_selector" type="select" label="Select analysis mode">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
289 <option value="illumina">1.Simple Illumina mode</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
290 <option value="full">2.Full list of options</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
291 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
292 <when value="illumina">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
293 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
294 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
295 <when value="full">
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
296 <param name="n" type="text" value="0.04" label="maximum edit distance if the value is integer, or the fraction of missing alignments given 2% uniform base error rate if float. In the latter case, the maximum edit distance is automatically chosen for different read lengths." help="aln -n; default=0.04"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
297 <param name="o" type="integer" value="1" label="maximum number or gap openings" help="aln -o; default=1"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
298 <param name="e" type="integer" value="-1" label="maximum number of gap extensions" help="aln -e; -1 disables long gaps and invokes k-difference mode; default=-1"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
299 <param name="i" type="integer" value="5" label="do not put an indel within this many bp towards the ends" help="aln -i; default=5"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
300 <param name="d" type="integer" value="10" label="maximum occurrences for extending a long deletion" help="aln -d; default=10"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
301 <param name="l" type="integer" value="32" label="seed length" help="aln -l; default=32"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
302 <param name="k" type="integer" value="2" label="maximum differences in the seed" help="aln -k; default=2"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
303 <param name="m" type="integer" value="2000000" label="maximum entries in the queue" help="aln -m; default=2000000"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
304 <param name="M" type="integer" value="3" label="mismatch penalty" help="aln -M; default=3"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
305 <param name="O" type="integer" value="11" label="gap open penalty" help="aln -O; default=11"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
306 <param name="E" type="integer" value="4" label="gap extension penalty" help="aln -E; default=4"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
307 <param name="R" type="integer" value="30" label="stop searching when there are more than this value of equally best hits" help="aln -R; default=30"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
308 <param name="q" type="integer" value="0" label="quality threshold for read trimming down to 35bp" help="aln -q; default=0"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
309 <param name="B" type="integer" optional="True" label="length of barcode" help="aln -B; optional parameter"/>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
310 <param name="L" type="float" optional="True" label="log-scaled gap penalty for long deletions" help="aln -L; optional parameter"/>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
311 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
312 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
313 </inputs>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
314
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
315 <outputs>
6
09a7281d24c5 planemo upload commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents: 5
diff changeset
316 <data format="bam" name="bam_output" label="${tool.name} on ${on_string} (mapped reads in BAM format)">
09a7281d24c5 planemo upload commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents: 5
diff changeset
317 <expand macro="dbKeyActionsBwa" />
09a7281d24c5 planemo upload commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents: 5
diff changeset
318 </data>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
319 </outputs>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
320
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
321 <tests>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
322 <test>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
323 <param name="reference_source_selector" value="history" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
324 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
325 <param name="input_type_selector" value="paired"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
326 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
327 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
328 <param name="analysis_type_selector" value="illumina"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
329 <output name="bam_output" ftype="bam" file="bwa-aln-test1.bam" lines_diff="2" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
330 </test>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
331 <test>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
332 <param name="reference_source_selector" value="history" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
333 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
13
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
334 <param name="input_type_selector" value="paired"/>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
335 <param name="fastq_input1" ftype="fastqsanger.gz" value="bwa-mem-fastq1.fq.gz"/>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
336 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
337 <param name="analysis_type_selector" value="illumina"/>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
338 <output name="bam_output" ftype="bam" file="bwa-aln-test1.bam" lines_diff="2" />
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
339 </test>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
340 <test>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
341 <param name="reference_source_selector" value="history" />
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
342 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
343 <param name="input_type_selector" value="paired_bam"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
344 <param name="bam_input" ftype="bam" value="bwa-aln-bam-input.bam"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
345 <param name="analysis_type_selector" value="illumina"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
346 <output name="bam_output" ftype="bam" file="bwa-aln-test2.bam" lines_diff="2" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
347 </test>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
348 <test>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
349 <param name="reference_source_selector" value="history" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
350 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
351 <param name="input_type_selector" value="paired"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
352 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
353 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
354 <param name="rg_selector" value="set"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
355 <param name="ID" value="rg1"/>
4
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
356 <param name="PL" value="CAPILLARY"/>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
357 <param name="analysis_type_selector" value="illumina"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
358 <output name="bam_output" ftype="bam" file="bwa-aln-test3.bam" lines_diff="2" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
359 </test>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
360 </tests>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
361 <help>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
362 **What is does**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
363
4
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
364 BWA is a software package for mapping low-divergent sequences against a large reference genome, such as the human genome. The bwa-aln algorithm is designed for Illumina sequence reads up to 100bp. For longer reads use BWA-MEM algorithm distributed as a separate Galaxy tool.
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
365
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
366 This Galaxy tool wraps bwa-aln, bwa-samse and -sampe modules of bwa read mapping tool:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
367
4
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
368 - **bwa aln** - actual mapper placing reads onto the reference sequence
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
369 - **bwa samse** - post-processor converting suffix array coordinates into genome coordinates in SAM format for single reads
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
370 - **bam sampe** - post-processor for paired reads
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
371
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
372 Galaxy implementation takes fastq or BAM (unaligned BAM) datasets as input and produces output in BAM (not SAM; in reality SAM produced by the bwa is converted to BAM on the fly by samtools view command) format, which can be further processed using various BAM utilities exiting in Galaxy (BAMTools, SAMTools, Picard).
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
373
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
374 -----
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
375
4
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
376 **Indices: Selecting reference genomes for BWA**
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
377
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
378 Galaxy wrapper for BWA allows you select between precomputed and user-defined indices for reference genomes using **Will you select a reference genome from your history or use a built-in index?** flag. This flag has two options:
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
379
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
380 1. **Use a built-in genome index** - when selected (this is default), Galaxy provides the user with **Select reference genome index** dropdown. Genomes listed in this dropdown have been pre-indexed with bwa index utility and are ready to be mapped against.
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
381 2. **Use a genome from the history and build index** - when selected, Galaxy provides the user with **Select reference genome sequence** dropdown. This dropdown is populated by all FASTA formatted files listed in your current history. If your genome of interest is uploaded into history it will be shown there. Selecting a genome from this dropdown will cause Galaxy to first transparently index it using `bwa index` command, and then run mapping with `bwa aln`.
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
382
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
383 If your genome of interest is not listed here you have two choices:
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
384
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
385 1. Contact galaxy team using **Help->Support** link at the top of the interface and let us know that an index needs to be added
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
386 2. Upload your genome of interest as a FASTA file to Galaxy history and selected **Use a genome from the history and build index** option.
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
387
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
388 -----
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 3
diff changeset
389
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
390 **Galaxy-specific option**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
391
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
392 Galaxy allows three levels of control over bwa-mem options provided by **Select analysis mode** menu option. These are:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
393
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
394 1. *Simple Illumina mode*: The simplest possible bwa mem application in which it alignes single or paired-end data to reference using default parameters. It is equivalent to the following command: bwa mem &lt;reference index&gt; &lt;fastq dataset1&gt; [fastq dataset2]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
395 2. *Full list of options*: Allows access to all options through Galaxy interface.
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
396
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
397 ------
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
398
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
399 **bwa-aln options**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
400
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
401 Each Galaxy parameter widget corresponds to command line flags listed below::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
402
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
403 -n NUM max #diff (int) or missing prob under 0.02 err rate (float) [0.04]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
404 -o INT maximum number or fraction of gap opens [1]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
405 -e INT maximum number of gap extensions, -1 for disabling long gaps [-1]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
406 -i INT do not put an indel within INT bp towards the ends [5]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
407 -d INT maximum occurrences for extending a long deletion [10]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
408 -l INT seed length [32]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
409 -k INT maximum differences in the seed [2]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
410 -m INT maximum entries in the queue [2000000]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
411 -M INT mismatch penalty [3]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
412 -O INT gap open penalty [11]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
413 -E INT gap extension penalty [4]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
414 -R INT stop searching when there are >INT equally best hits [30]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
415 -q INT quality threshold for read trimming down to 35bp [0]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
416 -B INT length of barcode
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
417 -L log-scaled gap penalty for long deletions
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
418 -N non-iterative mode: search for all n-difference hits (slooow)
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
419 -I the input is in the Illumina 1.3+ FASTQ-like format
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
420 -b the input read file is in the BAM format
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
421 -0 use single-end reads only (effective with -b)
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
422 -1 use the 1st read in a pair (effective with -b)
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
423 -2 use the 2nd read in a pair (effective with -b)
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
424
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
425 **bwa-samse options**::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
426
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
427 -a INT maximum insert size [500]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
428 -o INT maximum occurrences for one end [100000]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
429 -n INT maximum hits to output for paired reads [3]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
430 -N INT maximum hits to output for discordant pairs [10]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
431 -c FLOAT prior of chimeric rate (lower bound) [1.0e-05]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
432 -r STR read group header line [null]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
433
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
434 **bwa-sampe options**::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
435
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
436 -n INT maximum hits to output for paired reads [3]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
437 -r STR read group header line [null]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
438
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
439 @dataset_collections@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
440
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
441 @RG@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
442
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
443 @info@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
444 </help>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
445 <citations>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
446 <citation type="doi">10.1093/bioinformatics/btp324</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
447 <citation type="doi">10.1093/bioinformatics/btp698</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
448 </citations>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
449 </tool>