annotate bwa-mem.xml @ 1:c71dd035971e draft

Uploaded
author devteam
date Wed, 14 Jan 2015 13:51:07 -0500
parents ff1ae217ccc2
children e29bc5c169bc
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"?>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
2 <tool id="bwa_mem" name="BWA-MEM" version="0.1">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
3
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
4 <macros>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
5 <import>bwa_macros.xml</import>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
6 </macros>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
7
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
8 <requirements>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
9 <requirement type="package" version="0.7.10.039ea20639">bwa</requirement>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
10 <requirement type="package" version="1.1">samtools</requirement>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
11 </requirements>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
12 <description>- map medium and long reads (&gt; 100 bp) against reference genome</description>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
13 <command>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
14
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
15 #set $reference_fasta_filename = "localref.fa"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
16
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
17 #if str( $reference_source.reference_source_selector ) == "history":
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
18
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
19 ln -s "${reference_source.ref_file}" "${reference_fasta_filename}" &amp;&amp;
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
20
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
21 ## The following shell commands decide with of the BWA indexing algorithms (IS or BWTSW) will be run
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
22 ## depending ob the size of the input FASTA dataset
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
23
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
24 (
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
25 size=`stat -c %s "${reference_fasta_filename}" 2&gt;/dev/null`; ## Linux
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
26 if [ $? -eq 0 ];
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
27 then
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
28 if [ \$size -lt 2000000000 ];
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
29 then
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
30 bwa index -a is "${reference_fasta_filename}";
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
31 echo "Generating BWA index with is algorithm";
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
32 else
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
33 bwa index -a bwtsw "${reference_fasta_filename}";
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
34 echo "Generating BWA index with bwtsw algorithm";
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
35 fi;
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
36 fi;
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
37
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
38 eval \$(stat -s "${reference_fasta_filename}"); ## OSX
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
39 if [ $? -eq 0 ];
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
40 then
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
41 if [ \$st_size -lt 2000000000 ];
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
42 then
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
43 bwa index -a is "${reference_fasta_filename}";
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
44 echo "Generating BWA index with is algorithm";
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
45 else
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
46 bwa index -a bwtsw "${reference_fasta_filename}";
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
47 echo "Generating BWA index with bwtsw algorithm";
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
48 fi;
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
49 fi;
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
50 ) &amp;&amp;
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
51
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
52 #else:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
53 #set $reference_fasta_filename = str( $reference_source.ref_file.fields.path )
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
54 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
55
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
56 ## Begin BWA-MEM command line
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
57
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
58 bwa mem
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
59 -t "\${GALAXY_SLOTS:-1}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
60 -v 1 ## Verbosity is set to 1 (errors only)
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
61
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
62 #if str( $fastq_input.fastq_input_selector ) == "paired_iv": ## For interleaved fastq files set -p option
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
63 -p
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
64 #if str( $fastq_input.iset_stats ): ## check that insert statistics is used
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
65 -I "${fastq_input.iset_stats}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
66 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
67 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
68
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
69 #if str( $analysis_type.analysis_type_selector ) == "pacbio":
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
70 -x
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
71
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
72 #elif str( $analysis_type.analysis_type_selector ) == "full":
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
73
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
74 #if str( $analysis_type.algorithmic_options.algorithmic_options_selector ) == "True": ## Algorithmic options
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
75
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
76 -k "${analysis_type.algorithmic_options.k}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
77 -w "${analysis_type.algorithmic_options.w}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
78 -d "${analysis_type.algorithmic_options.d}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
79 -r "${analysis_type.algorithmic_options.r}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
80 -y "${analysis_type.algorithmic_options.y}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
81 -c "${analysis_type.algorithmic_options.c}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
82 -D "${analysis_type.algorithmic_options.D}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
83 -W "${analysis_type.algorithmic_options.W}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
84 -m "${analysis_type.algorithmic_options.m}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
85 ${analysis_type.algorithmic_options.S}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
86 ${analysis_type.algorithmic_options.P}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
87 ${analysis_type.algorithmic_options.e}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
88
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
89 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
90
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
91 #if str( $analysis_type.scoring_options.scoring_options_selector ) == "True": ## Scoring options
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
92
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
93 -A "${analysis_type.scoring_options.A}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
94 -B "${analysis_type.scoring_options.B}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
95 -O "${analysis_type.scoring_options.O}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
96 -E "${analysis_type.scoring_options.E}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
97 -L "${analysis_type.scoring_options.L}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
98 -U "${analysis_type.scoring_options.U}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
99
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
100 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
101
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
102 #if str( $analysis_type.io_options.io_options_selector ) == "True": ## IO options
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
103
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
104 -T "${analysis_type.io_options.T}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
105 -h "${analysis_type.io_options.h}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
106 ${analysis_type.io_options.a}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
107 ${analysis_type.io_options.C}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
108 ${analysis_type.io_options.V}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
109 ${analysis_type.io_options.Y}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
110 ${analysis_type.io_options.M}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
111
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
112 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
113
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
114 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
115
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
116 #if str( $rg.rg_selector ) == "True":
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
117 -R "@RG\tID:$rg.ID\tSM:$rg.SM"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
118 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
119
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
120 #if str( $fastq_input.fastq_input_selector ) == "paired":
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
121
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
122 #if str( $fastq_input.iset_stats ): ## check that insert statistics is used
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
123 -I "${fastq_input.iset_stats}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
124 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
125
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
126 "${reference_fasta_filename}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
127
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
128 "${fastq_input.fastq_input1}" "${fastq_input.fastq_input2}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
129
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
130 #elif str( $fastq_input.fastq_input_selector ) == "paired_collection":
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
131
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
132 #if str( $fastq_input.iset_stats ): ## check that insert statistics is used
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
133 -I "${fastq_input.iset_stats}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
134 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
135
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
136
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
137 "${reference_fasta_filename}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
138
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
139 "${fastq_input.fastq_input1.forward}" "${fastq_input.fastq_input1.reverse}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
140
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
141 #else:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
142
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
143
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
144 "${reference_fasta_filename}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
145
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
146 "${fastq_input.fastq_input1}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
147
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
148 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
149
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
150 | samtools view -Sb - > temporary_bam_file.bam &amp;&amp;
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
151
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
152 samtools sort -f temporary_bam_file.bam ${bam_output}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
153
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
154 </command>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
155
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
156 <inputs>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
157
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
158 <conditional name="reference_source">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
159 <param name="reference_source_selector" type="select" label="Load reference genome from">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
160 <option value="cached">Local cache</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
161 <option value="history">History</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
162 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
163 <when value="cached">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
164 <param name="ref_file" type="select" label="Using reference genome" help="Select genome from the list">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
165 <options from_data_table="bwa_mem_indexes">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
166 <filter type="sort_by" column="2" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
167 <validator type="no_options" message="No indexes are available" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
168 </options>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
169 <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
170 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
171 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
172 <when value="history">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
173 <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
174 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
175 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
176 <conditional name="fastq_input">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
177 <param name="fastq_input_selector" type="select" label="Single or Paired-end reads" help="Select between paired and single end data">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
178 <option value="paired">Paired</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
179 <option value="single">Single</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
180 <option value="paired_collection">Paired Collection</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
181 <option value="paired_iv">Paired Interleaved</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
182 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
183 <when value="paired">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
184 <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
185 <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
186 <param name="iset_stats" type="text" optional="True" size="10" label="Enter mean, standerd deviation, max, and min for insert lengths." help="-I; This parameter is only used for paired reads. Only mean is required while sd, max, and min will be inferred. Examples: both &quot;250&quot; and &quot;250,25&quot; will work while &quot;250,,10&quot; will not. See below for details.">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
187 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
188 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
189 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
190 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
191 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
192 <when value="single">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
193 <param name="fastq_input1" type="data" format="fastqsanger" label="Select fastq dataset" help="Specify dataset with single reads"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
194 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
195 <when value="paired_collection">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
196 <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
197 <param name="iset_stats" type="text" optional="True" size="10" label="Enter mean, standerd deviation, max, and min for insert lengths." help="-I; This parameter is only used for paired reads. Only mean is required while sd, max, and min will be inferred. Examples: both &quot;250&quot; and &quot;250,25&quot; will work while &quot;250,,10&quot; will not. See below for details.">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
198 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
199 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
200 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
201 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
202 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
203 <when value="paired_iv">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
204 <param name="fastq_input1" type="data" format="fastqsanger" label="Select fastq dataset" help="Specify dataset with interleaved reads"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
205 <param name="iset_stats" type="text" optional="True" size="10" label="Enter mean, standerd deviation, max, and min for insert lengths." help="-I; This parameter is only used for paired reads. Only mean is required while sd, max, and min will be inferred. Examples: both &quot;250&quot; and &quot;250,25&quot; will work while &quot;250,,10&quot; will not. See below for details.">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
206 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
207 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
208 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
209 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
210 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
211 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
212
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
213
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
214 <conditional name="rg">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
215 <param name="rg_selector" type="select" label="Set read groups information?" help="Specifying readgroup information can greatly simplify your downstream analyses by allowing combining multiple datasets. See help below for more details">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
216 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
217 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
218 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
219 <when value="set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
220 <param name="ID" type="text" value="readgroup1" size="20" label="Specify readgroup ID" help="This value must be unique among multiple samples in your experiment">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
221 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
222 <param name="SM" type="text" value="blood" size="20" label="Specify readgroup sample name (SM)" help="This value should be descriptive">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
223 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
224 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
225 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
226 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
227 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
228 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
229
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
230 <conditional name="analysis_type">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
231 <param name="analysis_type_selector" type="select" label="Select analysis mode">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
232 <option value="illumina">1.Simple Illumina mode</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
233 <option value="pacbio">2.PacBio mode</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
234 <option value="full">3.Full list of options</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
235 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
236 <when value="illumina">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
237 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
238 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
239 <when value="pacbio">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
240 <!-- do nothing. all magic happens within <command> tag -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
241 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
242 <when value="full">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
243 <conditional name="algorithmic_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
244 <param name="algorithmic_options_selector" type="select" label="Set algorithmic options?" help="Sets -k, -w, -d, -r, -y, -c, -D, -W, -m, -S, -P, and -e options.">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
245 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
246 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
247 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
248 <when value="set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
249 <param name="k" type="integer" value="19" label="minimum seed length" help="-k; default=19"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
250 <param name="w" type="integer" value="100" label="band width for banded alignment" help="-w; default=100"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
251 <param name="d" type="integer" value="100" label="off-diagonal X-dropoff" help="-d; default=100"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
252 <param name="r" type="float" value="1.5" label="look for internal seeds inside a seed longer than -k * THIS VALUE" help="-r; default=1.5"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
253 <param name="y" type="integer" value="0" label="find maximum exact matches (MEMs) longer than -k * -r with size less than THIS VALUE" help="-y; default=0"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
254 <param name="c" type="integer" value="500" label="skip seeds with more than that many occurrences" help="-c; default=500"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
255 <param name="D" type="float" value="0.5" label="drop chains shorter than this fraction of the longest overlapping chain" help="-D; default=0.5"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
256 <param name="W" type="integer" value="0" label="discard a chain if seeded bases shorter than" help="-W; default=0"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
257 <param name="m" type="integer" value="50" label="perform at most this many rounds of mate rescues for each read" help="-m; default=50"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
258 <param name="S" type="boolean" truevalue="-S" falsevalue="" label="skip mate rescue" help="-S"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
259 <param name="P" type="boolean" truevalue="-P" falsevalue="" label="skip pairing; mate rescue performed unless -S also in use" help="-P"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
260 <param name="e" type="boolean" truevalue="-e" falsevalue="" label="discard full-length exact matches" help="-e"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
261 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
262 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
263 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
264 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
265 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
266
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
267 <conditional name="scoring_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
268 <param name="scoring_options_selector" type="select" label="Set scoring options?" help="Sets -A, -B, -O, -E, -L, and -U options.">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
269 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
270 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
271 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
272 <when value="set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
273 <param name="A" type="integer" value="1" label="score for a sequence match" help="-A; scales options -T, -d, -B, -O, -E, -L, and -U; default=1"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
274 <param name="B" type="integer" value="4" label="penalty for mismatch" help="-B; default=4"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
275 <param name="O" type="text" value="6,6" label="gap open penalty for deletions and insertions" help="-O; default=6,6">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
276 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
277 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
278 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
279 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
280 <param name="E" type="text" value="1,1" label="gap extension penalty; a gap of size k cost &#39;-O + -E*k&#39; " help="-E; default=1,1">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
281 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
282 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
283 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
284 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
285 <param name="L" type="text" value="5,5" label="penalty for 5&#39;-end and 3&#39;-end clipping" help="-L; default=5,5">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
286 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
287 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
288 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
289 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
290 <param name="U" type="integer" value="17" label="penalty for an unpaired read pair" help="-U; default=17"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
291 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
292 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
293 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
294 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
295 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
296
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
297 <conditional name="io_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
298 <param name="io_options_selector" type="select" label="Set input/output options" help="Sets -T, -h, -a, -C, -V, -Y, and -M options.">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
299 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
300 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
301 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
302 <when value="set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
303 <param name="T" type="integer" value="30" label="minimum score to output" help="-T; default=30"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
304 <param name="h" type="integer" value="5" label="if there are this many hits with score >80% of the max score, output all in XA tag" help="-h; default=5"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
305 <param name="a" type="boolean" truevalue="-a" falsevalue="" label="output all alignments for single-ends or unpaired paired-ends" help="-a"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
306 <param name="C" type="boolean" truevalue="-C" falsevalue="" label="append FASTA/FASTQ comment to BAM output" help="-C"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
307 <param name="V" type="boolean" truevalue="-V" falsevalue="" label="output the reference FASTA header in the XR tag" help="-C"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
308 <param name="Y" type="boolean" truevalue="-Y" falsevalue="" label="use soft clipping for supplementary alignments" help="-Y"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
309 <param name="M" type="boolean" truevalue="-M" falsevalue="" label="mark shorter split hits as secondary" help="-M"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
310 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
311 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
312 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
313 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
314 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
315 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
316 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
317 </inputs>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
318
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
319 <outputs>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
320 <data format="bam" name="bam_output" label="${tool.name} on ${on_string} (mapped reads in BAM format)"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
321 </outputs>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
322
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
323 <tests>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
324 <test>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
325 <param name="reference_source_selector" value="history" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
326 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
327 <param name="fastq_input_selector" value="paired"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
328 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
329 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
330 <param name="analysis_type_selector" value="illumina"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
331 <output name="bam_output" ftype="bam" file="bwa-mem-test1.bam" lines_diff="2" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
332 </test>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
333 </tests>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
334 <stdio>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
335 <exit_code range="1:" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
336 </stdio>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
337 <help>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
338
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
339 **What is does**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
340
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
341 From http://arxiv.org/abs/1303.3997:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
342
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
343 BWA-MEM is a new alignment algorithm for aligning sequence reads or long query sequences against a large reference genome such as human.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
344 It automatically chooses between local and end-to-end alignments, supports paired-end reads and performs chimeric alignment.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
345 The algorithm is robust to sequencing errors and applicable to a wide range of sequence lengths from 70bp to a few megabases.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
346 For mapping 100bp sequences, BWA-MEM shows better performance than several state-of-art read aligners to date.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
347
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
348 It is best suited for mapping long (>70 nt) reads against large reference genomes.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
349
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
350 This Galaxy tool wraps bwa-mem module of bwa read mapping tool. Galaxy implementation takes fastq files as input and produces output in BAM (not SAM) format, which can be further processed using various BAM utilities exiting in Galaxy (BAMTools, SAMTools, Picard).
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
351
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
352 -----
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
353
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
354 **Galaxy-specific option**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
355
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
356 Galaxy allows four levels of control over bwa-mem options provided by **Select analysis mode** menu option. These are:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
357
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
358 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
359 2. *PacBio mode*: The mode adjusted specifically for mapping of long PacBio subreads. Equivalent to the following command: bwa mem -k17 -W40 -r10 -A1 -B1 -O1 -E1 -L0 &lt;reference index&gt; &lt;PacBio dataset in fastq format&gt;
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
360 3. *Full list of options*: Allows access to all options through Galaxy interface.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
361
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
362 ------
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
363
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
364 **BWA MEM options**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
365
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
366 Each Galaxy parameter widget corresponds to command line flags listed below:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
367
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
368 Algorithm options::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
369
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
370 -k INT minimum seed length [19]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
371 -w INT band width for banded alignment [100]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
372 -d INT off-diagonal X-dropoff [100]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
373 -r FLOAT look for internal seeds inside a seed longer than {-k} * FLOAT [1.5]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
374 -y INT find MEMs longer than {-k} * {-r} with size less than INT [0]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
375 -c INT skip seeds with more than INT occurrences [500]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
376 -D FLOAT drop chains shorter than FLOAT fraction of the longest overlapping chain [0.50]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
377 -W INT discard a chain if seeded bases shorter than INT [0]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
378 -m INT perform at most INT rounds of mate rescues for each read [50]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
379 -S skip mate rescue
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
380 -P skip pairing; mate rescue performed unless -S also in use
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
381 -e discard full-length exact matches
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
382
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
383 Scoring options::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
384
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
385 -A INT score for a sequence match, which scales options -TdBOELU unless overridden [1]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
386 -B INT penalty for a mismatch [4]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
387 -O INT[,INT] gap open penalties for deletions and insertions [6,6]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
388 -E INT[,INT] gap extension penalty; a gap of size k cost '{-O} + {-E}*k' [1,1]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
389 -L INT[,INT] penalty for 5'- and 3'-end clipping [5,5]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
390 -U INT penalty for an unpaired read pair [17]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
391
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
392 Input/output options::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
393
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
394 -p first query file consists of interleaved paired-end sequences
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
395 -R STR read group header line such as '@RG\tID:foo\tSM:bar' [null]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
396
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
397 -v INT verbose level: 1=error, 2=warning, 3=message, 4+=debugging [3]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
398 -T INT minimum score to output [30]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
399 -h INT if there are &lt;INT hits with score &gt;80% of the max score, output all in XA [5]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
400 -a output all alignments for SE or unpaired PE
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
401 -C append FASTA/FASTQ comment to SAM output
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
402 -V output the reference FASTA header in the XR tag
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
403 -Y use soft clipping for supplementary alignments
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
404 -M mark shorter split hits as secondary
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
405
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
406 -I FLOAT[,FLOAT[,INT[,INT]]]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
407 specify the mean, standard deviation (10% of the mean if absent), max
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
408 (4 sigma from the mean if absent) and min of the insert size distribution.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
409 FR orientation only. [inferred]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
410
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
411
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
412 @dataset_collections@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
413
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
414 @RG@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
415
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
416 @info@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
417
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
418
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
419
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
420 </help>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
421 <citations>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
422 <citation type="doi">10.1093/bioinformatics/btp324</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
423 <citation type="doi">10.1093/bioinformatics/btp698</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
424 <citation type="bibtex">@misc{1303.3997,
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
425 Author = {Heng Li},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
426 Title = {Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
427 Year = {2013},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
428 Eprint = {arXiv:1303.3997},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
429 url = {http://arxiv.org/abs/1303.3997},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
430 }</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
431 </citations>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
432 </tool>