annotate bwa.xml @ 3:607ca4b95837 draft

Uploaded
author devteam
date Fri, 20 Mar 2015 12:21:16 -0400
parents e29bc5c169bc
children ac30bfd3e2a8
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"?>
3
607ca4b95837 Uploaded
devteam
parents: 2
diff changeset
2 <tool id="bwa" name="Map with BWA" version="0.2.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>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
5 <import>bwa_macros.xml</import>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
6 <token name="@command_options@">
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
7 #if str( $analysis_type.analysis_type_selector ) == "full":
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
8 -n ${analysis_type.n}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
9 -o ${analysis_type.o}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
10 -e ${analysis_type.e}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
11 -i ${analysis_type.i}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
12 -d ${analysis_type.d}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
13 -l ${analysis_type.l}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
14 -k ${analysis_type.k}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
15 -m ${analysis_type.m}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
16 -M ${analysis_type.M}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
17 -O ${analysis_type.O}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
18 -E ${analysis_type.E}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
19 -R ${analysis_type.R}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
20 -q ${analysis_type.q}
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
21
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
22 #if str( $analysis_type.B ):
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
23 -B ${analysis_type.B}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
24 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
25
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
26 #if str( $analysis_type.L ):
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
27 -B ${analysis_type.L}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
28 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
29 #end if
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
30 </token>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
31 <token name="@read_group_options@">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
32 #if str( $rg.rg_selector ) == "set":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
33 @set_rg_string@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
34 -r '$rg_string'
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
35 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
36 </token>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
37
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
38 <xml name="advanced_pe_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
39 <param name="adv_pe_options_selector" type="select" label="Set advanced paired end options?" help="Provides additional controls">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
40 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
41 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
42 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
43 <when value="set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
44 <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
45 <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
46 <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
47 <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
48 <param name="c" type="float" value="0.00005" label="Prior of chimeric rate (lower bound)" help="sampe -c"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
49 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
50 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
51 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
52 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
53 </xml>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
54 <xml name="advanced_se_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
55 <param name="adv_se_options_selector" type="select" label="Set advanced single end options?" help="Provides additional controls">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
56 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
57 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
58 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
59 <when value="set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
60 <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
61 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
62 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
63 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
64 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
65 </xml>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
66 </macros>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
67
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
68 <requirements>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
69 <requirement type="package" version="0.7.10.039ea20639">bwa</requirement>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
70 <requirement type="package" version="1.1">samtools</requirement>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
71 </requirements>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
72 <stdio>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
73 <exit_code range="1:" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
74 </stdio>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
75 <command>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
76 #set $reference_fasta_filename = "localref.fa"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
77
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
78 #if str( $reference_source.reference_source_selector ) == "history":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
79 ln -s "${reference_source.ref_file}" "${reference_fasta_filename}" &amp;&amp;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
80
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
81 ## The following shell commands decide with of the BWA indexing algorithms (IS or BWTSW) will be run
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
82 ## depending ob the size of the input FASTA dataset
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
83 (
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
84 size=`stat -c %s "${reference_fasta_filename}" 2&gt;/dev/null`; ## Linux
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
85 if [ $? -eq 0 ];
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
86 then
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
87 if [ "\$size" -lt 2000000000 ];
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
88 then
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
89 bwa index -a is "${reference_fasta_filename}";
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
90 else
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
91 bwa index -a bwtsw "${reference_fasta_filename}";
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
92 fi;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
93 fi;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
94
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
95 eval \$(stat -s "${reference_fasta_filename}" 2&gt;/dev/null); ## OSX
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
96 if [ -n "\$st_size" ];
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
97 then
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
98 if [ "\$st_size" -lt 2000000000 ];
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
99 then
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
100 bwa index -a is "${reference_fasta_filename}";
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
101 echo "Generating BWA index with is algorithm";
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
102 else
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
103 bwa index -a bwtsw "${reference_fasta_filename}";
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
104 echo "Generating BWA index with bwtsw algorithm";
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
105 fi;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
106 fi;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
107 ) &amp;&amp;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
108 #else:
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
109 #set $reference_fasta_filename = str( $reference_source.ref_file.fields.path )
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
110 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
111
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
112 ## Begin bwa command line
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
113
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
114 ####### Fastq paired
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
115
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
116 #if str( $input_type.input_type_selector ) == "paired" or str( $input_type.input_type_selector ) == "paired_collection":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
117 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
118 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
119
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
120 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
121
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
122 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
123
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
124 #if str( $input_type.input_type_selector ) == "paired_collection":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
125 "${input_type.fastq_input1.forward}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
126 #else
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
127 "${input_type.fastq_input1}"
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 > first.sai &amp;&amp;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
131
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
132 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
133 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
134
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
135 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
136
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
137 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
138
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
139 #if str( $input_type.input_type_selector ) == "paired_collection":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
140 "${input_type.fastq_input1.reverse}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
141 #else
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
142 "${input_type.fastq_input2}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
143 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
144
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
145 > second.sai &amp;&amp;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
146
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
147 bwa sampe
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
148
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
149 #if str( $input_type.adv_pe_options.adv_pe_options_selector) == "True":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
150 -a ${$input_type.adv_pe_options.a}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
151 -o ${$input_type.adv_pe_options.o}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
152 -n ${$input_type.adv_pe_options.n}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
153 -N ${$input_type.adv_pe_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 #if str( $input_type.input_type_selector ) == "paired_collection":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
159 "${reference_fasta_filename}" first.sai second.sai "${input_type.fastq_input1.forward}" "${input_type.fastq_input1.reverse}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
160 #else:
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
161 "${reference_fasta_filename}" first.sai second.sai "${input_type.fastq_input1}" "${input_type.fastq_input2}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
162 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
163
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
164 ####### Fastq single
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
165
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
166 #elif str( $input_type.input_type_selector ) == "single":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
167 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
168 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
169
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
170 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
171
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
172 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
173 "${input_type.fastq_input1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
174 > first.sai &amp;&amp;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
175
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
176 bwa samse
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
177
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
178 #if str( $input_type.adv_se_options.adv_se_options_selector) == "True":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
179 -n ${$input_type.adv_se_options.n}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
180 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
181
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
182 @read_group_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
183
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
184 "${reference_fasta_filename}" first.sai "${input_type.fastq_input1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
185
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
186 ####### BAM paired
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
187
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
188 #elif str( $input_type.input_type_selector ) == "paired_bam":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
189 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
190 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
191 -b
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
192 -1
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
193
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
194 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
195
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
196 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
197 "${input_type.bam_input}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
198 > first.sai &amp;&amp;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
199
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
200 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
201 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
202 -b
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
203 -2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
204 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
205 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
206 "${input_type.bam_input}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
207 > second.sai &amp;&amp;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
208
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
209 bwa sampe
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
210
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
211 #if str( $input_type.adv_bam_pe_options.adv_pe_options_selector) == "True":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
212 -a ${$input_type.adv_bam_pe_options.a}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
213 -o ${$input_type.adv_bam_pe_options.o}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
214 -n ${$input_type.adv_bam_pe_options.n}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
215 -N ${$input_type.adv_bam_pe_options.N}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
216 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
217
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
218 @read_group_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
219
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
220 "${reference_fasta_filename}" first.sai second.sai "${input_type.bam_input}" "${input_type.bam_input}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
221
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
222 ####### Fastq single ------------ to do next
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
223
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
224 #elif str( $input_type.input_type_selector ) == "single_bam":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
225 bwa aln
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
226 -t "\${GALAXY_SLOTS:-1}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
227 -b
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
228 -0
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
229
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
230 @command_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
231
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
232 "${reference_fasta_filename}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
233 "${input_type.bam_input}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
234 > first.sai &amp;&amp;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
235
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
236 bwa samse
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
237
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
238 #if str( $input_type.adv_bam_se_options.adv_se_options_selector) == "True":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
239 -n ${$input_type.adv_bam_se_options.n}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
240 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
241
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
242 @read_group_options@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
243
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
244 "${reference_fasta_filename}" first.sai "${input_type.bam_input}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
245 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
246
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
247 | samtools view -Sb - > temporary_bam_file.bam &amp;&amp;
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
248
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
249 samtools sort -f temporary_bam_file.bam ${bam_output}
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
250 </command>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
251
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
252 <inputs>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
253
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
254 <conditional name="reference_source">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
255 <param name="reference_source_selector" type="select" label="Load reference genome from">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
256 <option value="cached">Local cache</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
257 <option value="history">History</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
258 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
259 <when value="cached">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
260 <param name="ref_file" type="select" label="Using reference genome" help="Select genome from the list">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
261 <options from_data_table="bwa_mem_indexes">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
262 <filter type="sort_by" column="2" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
263 <validator type="no_options" message="No indexes are available" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
264 </options>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
265 <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
266 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
267 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
268 <when value="history">
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
269 <param name="ref_file" type="data" format="fasta" label="Use the following dataset as the reference sequence" help="You can upload a FASTA sequence to the history and use it as reference" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
270 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
271 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
272 <conditional name="input_type">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
273 <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
274 <option value="paired">Paired fastq</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
275 <option value="paired_collection">Paired fastq collection</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
276 <option value="single">Single fastq</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
277 <option value="paired_bam">Paired BAM</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
278 <option value="single_bam">Single BAM</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
279 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
280 <when value="paired">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
281 <param name="fastq_input1" type="data" format="fastqsanger" label="Select first set of reads" help="Specify dataset with forward reads"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
282 <param name="fastq_input2" type="data" format="fastqsanger" label="Select second set of reads" help="Specify dataset with reverse reads"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
283 <conditional name="adv_pe_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
284
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
285 <expand macro="advanced_pe_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
286
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
287 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
288 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
289
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
290 <when value="paired_collection">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
291 <param name="fastq_input1" format="fastqsanger" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
292 <conditional name="adv_pe_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
293
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
294 <expand macro="advanced_pe_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
295
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
296 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
297 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
298
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
299 <when value="single">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
300 <param name="fastq_input1" type="data" format="fastqsanger" label="Select fastq dataset" help="Specify dataset with single reads"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
301 <conditional name="adv_se_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
302
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
303 <expand macro="advanced_se_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
304
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
305 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
306 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
307
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
308 <!-- 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
309
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
310 <when value="paired_bam">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
311 <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
312 <conditional name="adv_bam_pe_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
313
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
314 <expand macro="advanced_pe_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
315
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
316 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
317 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
318
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
319 <when value="single_bam">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
320 <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
321 <conditional name="adv_bam_se_options">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
322
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
323 <expand macro="advanced_se_options" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
324
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
325 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
326 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
327
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
328 </conditional>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
329
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
330 <expand macro="readgroup_params" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
331
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
332 <conditional name="analysis_type">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
333 <param name="analysis_type_selector" type="select" label="Select analysis mode">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
334 <option value="illumina">1.Simple Illumina mode</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
335 <option value="full">2.Full list of options</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
336 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
337 <when value="illumina">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
338 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
339 </when>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
340 <when value="full">
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
341 <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
342 <param name="o" type="integer" value="1" label="maximum number or gap openings" help="aln -o; default=1"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
343 <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
344 <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
345 <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
346 <param name="l" type="integer" value="32" label="seed length" help="aln -l; default=32"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
347 <param name="k" type="integer" value="2" label="maximum differences in the seed" help="aln -k; default=2"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
348 <param name="m" type="integer" value="2000000" label="maximum entries in the queue" help="aln -m; default=2000000"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
349 <param name="M" type="integer" value="3" label="mismatch penalty" help="aln -M; default=3"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
350 <param name="O" type="integer" value="11" label="gap open penalty" help="aln -O; default=11"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
351 <param name="E" type="integer" value="4" label="gap extension penalty" help="aln -E; default=4"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
352 <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
353 <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
354 <param name="B" type="integer" optional="True" label="length of barcode" help="aln -B; optional parameter"/>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
355 <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
356 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
357 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
358 </inputs>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
359
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
360 <outputs>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
361 <data format="bam" name="bam_output" label="${tool.name} on ${on_string} (mapped reads in BAM format)"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
362 </outputs>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
363
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
364 <tests>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
365 <test>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
366 <param name="reference_source_selector" value="history" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
367 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
368 <param name="input_type_selector" value="paired"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
369 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
370 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
371 <param name="analysis_type_selector" value="illumina"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
372 <output name="bam_output" ftype="bam" file="bwa-aln-test1.bam" lines_diff="2" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
373 </test>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
374 <test>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
375 <param name="reference_source_selector" value="history" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
376 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
377 <param name="input_type_selector" value="paired_bam"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
378 <param name="bam_input" ftype="bam" value="bwa-aln-bam-input.bam"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
379 <param name="analysis_type_selector" value="illumina"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
380 <output name="bam_output" ftype="bam" file="bwa-aln-test2.bam" lines_diff="2" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
381 </test>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
382 <test>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
383 <param name="reference_source_selector" value="history" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
384 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
385 <param name="input_type_selector" value="paired"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
386 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
387 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
388 <param name="rg_selector" value="set"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
389 <param name="ID" value="rg1"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
390 <param name="analysis_type_selector" value="illumina"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
391 <output name="bam_output" ftype="bam" file="bwa-aln-test3.bam" lines_diff="2" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
392 </test>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
393 </tests>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
394 <help>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
395 **What is does**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
396
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
397 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 separate Galaxy tool.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
398
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
399 This Galaxy tool wraps bwa-aln, bwa-samse and -sampe modules of bwa read mapping tool:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
400
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
401 - bwa aln - actual mapper placing reads onto the reference sequence
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
402 - bwa samse - post-processor converting suffix array coordinates into genome coordinates in SAM format for single reads
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
403 - bam sampe - post-processor for paired reads
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
404
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
405 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
406
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
407 -----
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
408
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
409 **Galaxy-specific option**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
410
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
411 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
412
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
413 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
414 2. *Full list of options*: Allows access to all options through Galaxy interface.
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
415
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
416 ------
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
417
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
418 **bwa-aln options**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
419
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
420 Each Galaxy parameter widget corresponds to command line flags listed below::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
421
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
422 -n NUM max #diff (int) or missing prob under 0.02 err rate (float) [0.04]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
423 -o INT maximum number or fraction of gap opens [1]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
424 -e INT maximum number of gap extensions, -1 for disabling long gaps [-1]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
425 -i INT do not put an indel within INT bp towards the ends [5]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
426 -d INT maximum occurrences for extending a long deletion [10]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
427 -l INT seed length [32]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
428 -k INT maximum differences in the seed [2]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
429 -m INT maximum entries in the queue [2000000]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
430 -M INT mismatch penalty [3]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
431 -O INT gap open penalty [11]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
432 -E INT gap extension penalty [4]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
433 -R INT stop searching when there are >INT equally best hits [30]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
434 -q INT quality threshold for read trimming down to 35bp [0]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
435 -B INT length of barcode
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
436 -L log-scaled gap penalty for long deletions
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
437 -N non-iterative mode: search for all n-difference hits (slooow)
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
438 -I the input is in the Illumina 1.3+ FASTQ-like format
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
439 -b the input read file is in the BAM format
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
440 -0 use single-end reads only (effective with -b)
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
441 -1 use the 1st read in a pair (effective with -b)
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
442 -2 use the 2nd read in a pair (effective with -b)
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
443
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
444 **bwa-samse options**::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
445
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
446 -a INT maximum insert size [500]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
447 -o INT maximum occurrences for one end [100000]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
448 -n INT maximum hits to output for paired reads [3]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
449 -N INT maximum hits to output for discordant pairs [10]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
450 -c FLOAT prior of chimeric rate (lower bound) [1.0e-05]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
451 -r STR read group header line [null]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
452
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
453 **bwa-sampe options**::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
454
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
455 -n INT maximum hits to output for paired reads [3]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
456 -r STR read group header line [null]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
457
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
458 @dataset_collections@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
459
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
460 @RG@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
461
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
462 @info@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
463 </help>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
464 <citations>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
465 <citation type="doi">10.1093/bioinformatics/btp324</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
466 <citation type="doi">10.1093/bioinformatics/btp698</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
467 </citations>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
468 </tool>