Mercurial > repos > devteam > bwa

annotate bwa-mem.xml @ 15:4d82cf59895e draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit dead2eaab58698fe02b8bb068f8d859eba551e8d
author devteam
date Sat, 30 Sep 2017 13:37:06 -0400
parents 051eba708f43
children be4e38d127ae
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"?>
15
4d82cf59895e planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit dead2eaab58698fe02b8bb068f8d859eba551e8d
devteam
parents: 14
diff changeset
2 <tool id="bwa_mem" name="Map with BWA-MEM" version="@VERSION@.0">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
3 <description>- map medium and long reads (&gt; 100 bp) against reference genome</description>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
4 <macros>
5
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
5 <import>read_group_macros.xml</import>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
6 <import>bwa_macros.xml</import>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
7 </macros>
7
d8c9597bfb09 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit ddb8bdb9d62399f086b06b3469450d0aad2113bd
devteam
parents: 6
diff changeset
8 <expand macro="requirements" />
6
09a7281d24c5 planemo upload commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents: 5
diff changeset
9 <expand macro="stdio" />
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
10 <command>
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
11 <![CDATA[
11
546ada4a9f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 610045611b00099e0a24183c8cf33aebfa9635cf
devteam
parents: 10
diff changeset
12 @set_reference_fasta_filename@
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
13
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
14 ## Begin BWA-MEM command line
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
15
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
16 bwa mem
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
17 -t "\${GALAXY_SLOTS:-1}"
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
18 -v 1 ## Verbosity is set to 1 (errors only)
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
19
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
20 #if str( $fastq_input.fastq_input_selector ) == "paired_iv": ## For interleaved fastq files set -p option
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
21 -p
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
22 #if str( $fastq_input.iset_stats ): ## check that insert statistics is used
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
23 -I "${fastq_input.iset_stats}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
24 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
25 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
26
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
27 #if str( $analysis_type.analysis_type_selector ) == "pacbio":
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
28 -x pacbio
12
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
29 #elif str( $analysis_type.analysis_type_selector ) == "ont2d":
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
30 -x ont2d
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
31 #elif str( $analysis_type.analysis_type_selector ) == "intractg":
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
32 -x intractg
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
33 #elif str( $analysis_type.analysis_type_selector ) == "full":
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
34 #if str( $analysis_type.algorithmic_options.algorithmic_options_selector ) == "set": ## Algorithmic options
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
35 -k "${analysis_type.algorithmic_options.k}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
36 -w "${analysis_type.algorithmic_options.w}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
37 -d "${analysis_type.algorithmic_options.d}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
38 -r "${analysis_type.algorithmic_options.r}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
39 -y "${analysis_type.algorithmic_options.y}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
40 -c "${analysis_type.algorithmic_options.c}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
41 -D "${analysis_type.algorithmic_options.D}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
42 -W "${analysis_type.algorithmic_options.W}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
43 -m "${analysis_type.algorithmic_options.m}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
44 ${analysis_type.algorithmic_options.S}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
45 ${analysis_type.algorithmic_options.P}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
46 ${analysis_type.algorithmic_options.e}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
47 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
48
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
49 #if str( $analysis_type.scoring_options.scoring_options_selector ) == "set": ## Scoring options
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
50 -A "${analysis_type.scoring_options.A}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
51 -B "${analysis_type.scoring_options.B}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
52 -O "${analysis_type.scoring_options.O}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
53 -E "${analysis_type.scoring_options.E}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
54 -L "${analysis_type.scoring_options.L}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
55 -U "${analysis_type.scoring_options.U}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
56 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
57
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
58 #if str( $analysis_type.io_options.io_options_selector ) == "set": ## IO options
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
59 -T "${analysis_type.io_options.T}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
60 -h "${analysis_type.io_options.h}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
61 ${analysis_type.io_options.a}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
62 ${analysis_type.io_options.C}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
63 ${analysis_type.io_options.V}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
64 ${analysis_type.io_options.Y}
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
65 ${analysis_type.io_options.M}
15
4d82cf59895e planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit dead2eaab58698fe02b8bb068f8d859eba551e8d
devteam
parents: 14
diff changeset
66 ${analysis_type.io_options.five}
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
67 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
68
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
69 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
70
5
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
71 ## Handle read group options...
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
72 @define_read_group_helpers@
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
73 #if str( $fastq_input.fastq_input_selector ) == "paired":
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
74 #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
75 #else:
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
76 #set $rg_auto_name = $read_group_name_default($fastq_input.fastq_input1)
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
77 #end if
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
78 @set_use_rg_var@
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
79 @set_read_group_vars@
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
80 #if $use_rg
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
81 @set_rg_string@
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
82 -R '$rg_string'
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
83 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
84
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
85 #if str( $fastq_input.fastq_input_selector ) == "paired":
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
86 #if str( $fastq_input.iset_stats ): ## check that insert statistics is used
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
87 -I "${fastq_input.iset_stats}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
88 #end if
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
89
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
90 "${reference_fasta_filename}"
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
91 "${fastq_input.fastq_input1}" "${fastq_input.fastq_input2}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
92 #elif str( $fastq_input.fastq_input_selector ) == "paired_collection":
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
93 #if str( $fastq_input.iset_stats ): ## check that insert statistics is used
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
94 -I "${fastq_input.iset_stats}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
95 #end if
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
96
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
97 "${reference_fasta_filename}"
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
98 "${fastq_input.fastq_input1.forward}" "${fastq_input.fastq_input1.reverse}"
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
99 #else:
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
100 "${reference_fasta_filename}"
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
101 "${fastq_input.fastq_input1}"
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
102 #end if
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
103
12
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
104 | samtools sort -O bam -o '$bam_output'
10
6069ffa8b240 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 653588a3c59eb9965097f87c99e6fca1806365cd
devteam
parents: 9
diff changeset
105 ]]>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
106 </command>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
107
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
108 <inputs>
11
546ada4a9f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 610045611b00099e0a24183c8cf33aebfa9635cf
devteam
parents: 10
diff changeset
109 <expand macro="reference_source_conditional" />
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
110 <conditional name="fastq_input">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
111 <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
112 <option value="paired">Paired</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
113 <option value="single">Single</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
114 <option value="paired_collection">Paired Collection</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
115 <option value="paired_iv">Paired Interleaved</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
116 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
117 <when value="paired">
14
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
118 <param name="fastq_input1" type="data" format="fastqsanger,fastqsanger.gz,fasta" label="Select first set of reads" help="Specify dataset with forward reads"/>
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
119 <param name="fastq_input2" type="data" format="fastqsanger,fastqsanger.gz,fasta" label="Select second set of reads" help="Specify dataset with reverse reads"/>
7
d8c9597bfb09 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit ddb8bdb9d62399f086b06b3469450d0aad2113bd
devteam
parents: 6
diff changeset
120 <param name="iset_stats" type="text" optional="True" 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 &quot;250&quot; and &quot;250,25&quot; will work while &quot;250,,10&quot; will not. See below for details.">
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
121 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
122 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
123 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
124 </param>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
125 </when>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
126 <when value="single">
14
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
127 <param name="fastq_input1" type="data" format="fastqsanger,fastqsanger.gz,fasta" label="Select fastq dataset" help="Specify dataset with single reads"/>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
128 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
129 <when value="paired_collection">
14
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
130 <param name="fastq_input1" format="fastqsanger,fastqsanger.gz,fasta" type="data_collection" collection_type="paired" label="Select a paired collection" help="See help section for an explanation of dataset collections"/>
7
d8c9597bfb09 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit ddb8bdb9d62399f086b06b3469450d0aad2113bd
devteam
parents: 6
diff changeset
131 <param name="iset_stats" type="text" optional="True" 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 &quot;250&quot; and &quot;250,25&quot; will work while &quot;250,,10&quot; will not. See below for details.">
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
132 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
133 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
134 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
135 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
136 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
137 <when value="paired_iv">
14
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
138 <param name="fastq_input1" type="data" format="fastqsanger,fastqsanger.gz,fasta" label="Select fastq dataset" help="Specify dataset with interleaved reads"/>
7
d8c9597bfb09 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit ddb8bdb9d62399f086b06b3469450d0aad2113bd
devteam
parents: 6
diff changeset
139 <param name="iset_stats" type="text" optional="True" 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 &quot;250&quot; and &quot;250,25&quot; will work while &quot;250,,10&quot; will not. See below for details.">
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
140 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
141 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
142 </sanitizer>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
143 </param>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
144 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
145 </conditional>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
146
5
fbf460831036 planemo upload commit d0e3412c58bd3bdc1a483a1e2f7f9c2aa5c87a1f-dirty
devteam
parents: 4
diff changeset
147 <expand macro="read_group_conditional" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
148
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
149 <conditional name="analysis_type">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
150 <param name="analysis_type_selector" type="select" label="Select analysis mode">
4
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
151 <option value="illumina">1.Simple Illumina mode</option>
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
152 <option value="pacbio">2.PacBio mode (-x pacbio)</option>
12
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
153 <option value="ont2d">3.Nanopore 2D-reads mode (-x ont2d)</option>
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
154 <option value="intractg">4.Intra-species contigs mode (-x intractg)</option>
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
155 <option value="full">5.Full list of options</option>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
156 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
157 <when value="illumina">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
158 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
159 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
160 <when value="pacbio">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
161 <!-- do nothing. all magic happens within <command> tag -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
162 </when>
12
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
163 <when value="ont2d">
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
164 <!-- do nothing. all magic happens within <command> tag -->
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
165 </when>
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
166 <when value="intractg">
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
167 <!-- do nothing. all magic happens within <command> tag -->
bd3a1e0de84c planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 13b3ab10ff1a2e4be39b44f885df6bbe7c23099a
devteam
parents: 11
diff changeset
168 </when>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
169 <when value="full">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
170 <conditional name="algorithmic_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
171 <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
172 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
173 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
174 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
175 <when value="set">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
176 <param name="k" type="integer" value="19" label="Minimum seed length" help="-k; default=19"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
177 <param name="w" type="integer" value="100" label="Band width for banded alignment" help="-w; default=100"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
178 <param name="d" type="integer" value="100" label="Off-diagonal X-dropoff" help="-d; default=100"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
179 <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" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
180 <param name="y" type="integer" value="20" label="Seed occurrence for the 3rd round seeding" help="-y; default=20" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
181 <param name="c" type="integer" value="500" label="Skip seeds with more than that many occurrences" help="-c; default=500"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
182 <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"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
183 <param name="W" type="integer" value="0" label="Discard a chain if seeded bases shorter than THIS VALUE" help="-W; default=0"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
184 <param name="m" type="integer" value="50" label="Perform at most this many rounds of mate rescues for each read" help="-m; default=50"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
185 <param name="S" type="boolean" truevalue="-S" falsevalue="" label="Skip mate rescue" help="-S"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
186 <param name="P" type="boolean" truevalue="-P" falsevalue="" label="Skip pairing; mate rescue performed unless -S also in use" help="-P"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
187 <param name="e" type="boolean" truevalue="-e" falsevalue="" label="Discard full-length exact matches" help="-e"/>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
188 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
189 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
190 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
191 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
192 </conditional>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
193
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
194 <conditional name="scoring_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
195 <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
196 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
197 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
198 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
199 <when value="set">
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
200 <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"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
201 <param name="B" type="integer" value="4" label="Penalty for a mismatch" help="-B; default=4"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
202 <param name="O" type="text" value="6,6" label="Gap open penalties for deletions and insertions" help="-O; default=6,6">
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
203 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
204 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
205 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
206 </param>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
207 <param name="E" type="text" value="1,1" label="Gap extension penalties; a gap of size k cost &#39;-O + -E*k&#39;. 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
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
208 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
209 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
210 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
211 </param>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
212 <param name="L" type="text" value="5,5" label="Penalties for 5&#39;-end and 3&#39;-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
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
213 <sanitizer invalid_char="">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
214 <valid initial="string.digits"><add value=","/> </valid>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
215 </sanitizer>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
216 </param>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
217 <param name="U" type="integer" value="17" label="Penalty for an unpaired read pair" help="-U; default=17"/>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
218 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
219 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
220 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
221 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
222 </conditional>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
223
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
224 <conditional name="io_options">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
225 <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
226 <option value="set">Set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
227 <option value="do_not_set" selected="True">Do not set</option>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
228 </param>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
229 <when value="set">
15
4d82cf59895e planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit dead2eaab58698fe02b8bb068f8d859eba551e8d
devteam
parents: 14
diff changeset
230 <param name="five" argument="-5" type="boolean" truevalue="-5" falsevalue="" label="For split alignment, take alignment with smallest coordinate as primary" help="Useful for HiC data"/>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
231 <param name="T" type="integer" value="30" label="Minimum score to output" help="-T; default=30"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
232 <param name="h" type="integer" value="5" label="If there are less than THIS VALUE hits with score &gt;80% of the max score, output them all in the XA tag" help="-h; default=5" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
233 <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"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
234 <param name="C" type="boolean" truevalue="-C" falsevalue="" label="Append FASTA/FASTQ comment to BAM output" help="-C"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
235 <param name="V" type="boolean" truevalue="-V" falsevalue="" label="Output the reference FASTA header in the XR tag" help="-C"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
236 <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" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
237 <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&lt;1.96 compatibility" />
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
238 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
239 <when value="do_not_set">
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
240 <!-- do nothing -->
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
241 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
242 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
243 </when>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
244 </conditional>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
245 </inputs>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
246
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
247 <outputs>
6
09a7281d24c5 planemo upload commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents: 5
diff changeset
248 <data format="bam" name="bam_output" label="${tool.name} on ${on_string} (mapped reads in BAM format)">
09a7281d24c5 planemo upload commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents: 5
diff changeset
249 <expand macro="dbKeyActionsBwaMem" />
09a7281d24c5 planemo upload commit 00a7926c285bc4a339bd7deebf40b28f39c7d947-dirty
devteam
parents: 5
diff changeset
250 </data>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
251 </outputs>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
252
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
253 <tests>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
254 <test>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
255 <param name="reference_source_selector" value="history" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
256 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
257 <param name="fastq_input_selector" value="paired"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
258 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
259 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
260 <param name="analysis_type_selector" value="illumina"/>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
261 <output name="bam_output" ftype="bam" file="bwa-mem-test1.bam" lines_diff="2" />
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
262 </test>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
263 <test>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
264 <param name="reference_source_selector" value="history" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
265 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
14
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
266 <param name="fastq_input_selector" value="single"/>
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
267 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fasta1.fa"/>
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
268 <param name="analysis_type_selector" value="illumina"/>
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
269 <output name="bam_output" ftype="bam" file="bwa-mem-test1-fasta.bam" lines_diff="2" />
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
270 </test>
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
271 <test>
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
272 <param name="reference_source_selector" value="history" />
051eba708f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 33c17fa6a2d8c15972d93031b2ca4182cdf31439
devteam
parents: 13
diff changeset
273 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
13
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
274 <param name="fastq_input_selector" value="paired"/>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
275 <param name="fastq_input1" ftype="fastqsanger.gz" value="bwa-mem-fastq1.fq.gz"/>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
276 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
277 <param name="analysis_type_selector" value="illumina"/>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
278 <output name="bam_output" ftype="bam" file="bwa-mem-test1.bam" lines_diff="2" />
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
279 </test>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
280 <test>
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
281 <param name="reference_source_selector" value="history" />
53646aaaafef planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 3442b18860bc4398d917e3a35943351a9cf031ea
devteam
parents: 12
diff changeset
282 <param name="ref_file" ftype="fasta" value="bwa-mem-mt-genome.fa"/>
11
546ada4a9f43 planemo upload for repository https://github.com/galaxyproject/tools-devteam/tree/master/tools/bwa commit 610045611b00099e0a24183c8cf33aebfa9635cf
devteam
parents: 10
diff changeset
283 <param name="index_a" value="is"/>
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
284 <param name="fastq_input_selector" value="paired"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
285 <param name="fastq_input1" ftype="fastqsanger" value="bwa-mem-fastq1.fq"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
286 <param name="fastq_input2" ftype="fastqsanger" value="bwa-mem-fastq2.fq"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
287 <param name="rg_selector" value="set"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
288 <param name="ID" value="rg1"/>
4
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
289 <param name="PL" value="CAPILLARY"/>
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
290 <param name="LB" value="AARDVARK-1" />
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
291 <param name="analysis_type_selector" value="illumina"/>
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
292 <output name="bam_output" ftype="bam" file="bwa-mem-test2.bam" lines_diff="2" />
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
293 </test>
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
294 </tests>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
295 <help>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
296 **What is does**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
297
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
298 From http://arxiv.org/abs/1303.3997:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
299
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
300 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
301 It automatically chooses between local and end-to-end alignments, supports paired-end reads and performs chimeric alignment.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
302 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
303 For mapping 100bp sequences, BWA-MEM shows better performance than several state-of-art read aligners to date.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
304
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
305 It is best suited for mapping long (>70 nt) reads against large reference genomes.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
306
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
307 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
308
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
309 -----
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
310
4
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
311 **Indices: Selecting reference genomes for BWA**
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
312
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
313 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
314
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
315 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
316 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
317
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
318 If your genome of interest is not listed here you have two choices:
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
319
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
320 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
321 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
322
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
323 -----
ac30bfd3e2a8 planemo upload commit a50a3947aebc8a1d11bac39599f4efd8ed9a3bd5
devteam
parents: 2
diff changeset
324
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
325 **Galaxy-specific option**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
326
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
327 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
328
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
329 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
330 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
331 3. *Full list of options*: Allows access to all options through Galaxy interface.
2
e29bc5c169bc Uploaded
devteam
parents: 1
diff changeset
332
0
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
333 ------
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
334
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
335 **BWA MEM options**
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
336
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
337 Each Galaxy parameter widget corresponds to command line flags listed below:
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
338
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
339 Algorithm options::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
340
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
341 -k INT minimum seed length [19]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
342 -w INT band width for banded alignment [100]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
343 -d INT off-diagonal X-dropoff [100]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
344 -r FLOAT look for internal seeds inside a seed longer than {-k} * FLOAT [1.5]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
345 -y INT find MEMs longer than {-k} * {-r} with size less than INT [0]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
346 -c INT skip seeds with more than INT occurrences [500]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
347 -D FLOAT drop chains shorter than FLOAT fraction of the longest overlapping chain [0.50]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
348 -W INT discard a chain if seeded bases shorter than INT [0]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
349 -m INT perform at most INT rounds of mate rescues for each read [50]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
350 -S skip mate rescue
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
351 -P skip pairing; mate rescue performed unless -S also in use
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
352 -e discard full-length exact matches
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
353
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
354 Scoring options::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
355
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
356 -A INT score for a sequence match, which scales options -TdBOELU unless overridden [1]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
357 -B INT penalty for a mismatch [4]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
358 -O INT[,INT] gap open penalties for deletions and insertions [6,6]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
359 -E INT[,INT] gap extension penalty; a gap of size k cost '{-O} + {-E}*k' [1,1]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
360 -L INT[,INT] penalty for 5'- and 3'-end clipping [5,5]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
361 -U INT penalty for an unpaired read pair [17]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
362
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
363 Input/output options::
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
364
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
365 -p first query file consists of interleaved paired-end sequences
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
366 -R STR read group header line such as '@RG\tID:foo\tSM:bar' [null]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
367
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
368 -v INT verbose level: 1=error, 2=warning, 3=message, 4+=debugging [3]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
369 -T INT minimum score to output [30]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
370 -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
371 -a output all alignments for SE or unpaired PE
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
372 -C append FASTA/FASTQ comment to SAM output
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
373 -V output the reference FASTA header in the XR tag
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
374 -Y use soft clipping for supplementary alignments
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
375 -M mark shorter split hits as secondary
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
376
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
377 -I FLOAT[,FLOAT[,INT[,INT]]]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
378 specify the mean, standard deviation (10% of the mean if absent), max
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
379 (4 sigma from the mean if absent) and min of the insert size distribution.
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
380 FR orientation only. [inferred]
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
381
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
382 @dataset_collections@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
383
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
384 @RG@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
385
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
386 @info@
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
387 </help>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
388 <citations>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
389 <citation type="doi">10.1093/bioinformatics/btp324</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
390 <citation type="doi">10.1093/bioinformatics/btp698</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
391 <citation type="bibtex">@misc{1303.3997,
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
392 Author = {Heng Li},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
393 Title = {Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
394 Year = {2013},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
395 Eprint = {arXiv:1303.3997},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
396 url = {http://arxiv.org/abs/1303.3997},
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
397 }</citation>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
398 </citations>
ff1ae217ccc2 Uploaded
devteam
parents:
diff changeset
399 </tool>