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