annotate snap_caller.xml @ 23:5db0545b9004 draft

update to v0.1.7.3
author wolma
date Thu, 21 Jul 2016 03:55:49 -0400
parents c46406466625
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
1 <tool id="read_alignment" name="SNAP Read Alignment" version="0.1.7.3">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
2 <description>Map sequence reads to a reference genome using SNAP</description>
2
72d20758ba2c final upload
wolma
parents: 1
diff changeset
3 <macros>
72d20758ba2c final upload
wolma
parents: 1
diff changeset
4 <import>toolshed_macros.xml</import>
72d20758ba2c final upload
wolma
parents: 1
diff changeset
5 </macros>
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
6 <expand macro="requirements" />
21
c46406466625 Uploaded
wolma
parents: 9
diff changeset
7 <version_command>python3 -m MiModD version -q</version_command>
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
8 <command>
21
c46406466625 Uploaded
wolma
parents: 9
diff changeset
9 python3 -m MiModD snap-batch -s
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
10 ## SNAP calls (considering different cases)
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
11
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
12 #for $i in $datasets
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
13 "snap ${i.mode_choose.mode} '$ref_genome'
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
14 #if $str($i.mode_choose.mode) == "paired" and $str($i.mode_choose.input.iformat) in ("fastq", "gz"):
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
15 '${i.mode_choose.input.ifile1}' '${i.mode_choose.input.ifile2}'
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
16 #else:
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
17 '${i.mode_choose.input.ifile}'
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
18 #end if
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
19 --ofile '$outputfile' --iformat ${i.mode_choose.input.iformat} --oformat $oformat
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
20 --idx-seedsize '$set.seedsize'
6
85214e4428fd upgrade to mimodd version 0.1.7.0
Wolfgang Maier
parents: 5
diff changeset
21 --idx-slack '$set.slack' --maxseeds '$set.maxseeds' --maxhits '$set.maxhits' --clipping $set.clipping --maxdist '$set.maxdist' --confdiff '$set.confdiff' --confadapt '$set.confadpt'
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
22 #if $i.mode_choose.input.header:
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
23 --header '${i.mode_choose.input.header}'
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
24 #end if
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
25 #if $str($i.mode_choose.mode) == "paired":
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
26 --spacing '$set.sp_min' '$set.sp_max'
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
27 #end if
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
28 #if $str($set.selectivity) != "off":
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
29 --selectivity '$set.selectivity'
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
30 #end if
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
31 #if $str($set.filter_output) != "off":
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
32 --filter-output $set.filter_output
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
33 #end if
4
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
34 #if $str($set.sort) == "off":
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
35 --no-sort
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
36 #end if
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
37 #if $str($set.mmatch_notation) != "general":
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
38 -X
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
39 #end if
4
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
40 #if $set.discard_overlapping_mates:
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
41 --discard-overlapping-mates
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
42 ## remove ',' (and possibly adjacent whitespace) and replace with ' '
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
43 '#echo ("' '".join($set.discard_overlapping_mates.replace(" ", "").split(',')))#'
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
44 #end if
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
45 --verbose
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
46 "
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
47 #end for
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
48 </command>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
49
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
50 <inputs>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
51 ## mandatory arguments (and mode-conditionals)
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
52
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
53 <param format="fasta" help="The fasta reference genome that SNAP should align reads against." label="reference genome" name="ref_genome" type="data" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
54
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
55 <repeat default="1" min="1" name="datasets" title="datasets">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
56 <conditional name="mode_choose">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
57 <param help="Reads obtained from single-end sequencing runs should be aligned in 'single' mode, paired-end reads in 'paired' mode. **WARNING**: if the read input file is in SAM/BAM format, the current version of this tool will **not** verify the mode and may produce erroneous alignments with wrong settings!" label="choose mode" name="mode" type="select">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
58 <option value="single">single-end</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
59 <option value="paired">paired-end</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
60 </param>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
61
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
62 <when value="single">
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
63 <conditional name="input">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
64 <param label="input file format" name="iformat" type="select">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
65 <option value="bam">BAM</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
66 <option value="sam">SAM</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
67 <option value="gz">gz</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
68 <option value="fastq">fastq</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
69 </param>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
70 <when value="bam">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
71 <param format="bam" label="input file" name="ifile" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
72 <param format="sam" label="custom header file" name="header" optional="true" type="data" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
73 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
74 <when value="sam">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
75 <param format="sam" label="input file" name="ifile" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
76 <param format="sam" label="custom header file" name="header" optional="true" type="data" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
77 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
78 <when value="gz">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
79 <param label="input file" name="ifile" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
80 <param format="sam" label="header file" name="header" type="data" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
81 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
82 <when value="fastq">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
83 <param format="fastq" label="input file" name="ifile" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
84 <param format="sam" label="header file" name="header" type="data" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
85 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
86 </conditional>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
87 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
88 <when value="paired">
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
89 <conditional name="input">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
90 <param label="input file format" name="iformat" type="select">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
91 <option value="bam">BAM</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
92 <option value="sam">SAM</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
93 <option value="gz">gz</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
94 <option value="fastq">fastq</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
95 </param>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
96 <when value="bam">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
97 <param format="bam" label="input file" name="ifile" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
98 <param format="sam" label="custom header file" name="header" optional="true" type="data" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
99 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
100 <when value="sam">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
101 <param format="sam" label="input file" name="ifile" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
102 <param format="sam" label="custom header file" name="header" optional="true" type="data" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
103 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
104 <when value="fastq">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
105 <param format="fastq" label="inputfile with the first set of reads of paired-end data" name="ifile1" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
106 <param format="fastq" label="inputfile with the second set of reads of paired-end data" name="ifile2" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
107 <param format="sam" help="required" label="header file" name="header" type="data" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
108 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
109 <when value="gz">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
110 <param label="inputfile with the first set of reads of paired-end data" name="ifile1" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
111 <param label="inputfile with the second set of reads of paired-end data" name="ifile2" type="data" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
112 <param format="sam" help="required" label="header file" name="header" type="data" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
113 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
114 </conditional>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
115 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
116 </conditional>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
117 </repeat>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
118
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
119 <param label="output file format" name="oformat" type="select">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
120 <option value="bam">BAM</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
121 <option value="sam">SAM</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
122 </param>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
123
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
124 ## optional arguments
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
125
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
126 <conditional name="set">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
127 <param help="This section lets you specify the detailed parameter settings for the SNAP aligner. Only change them if you know what you are doing, i.e., read the documentation first." label="further parameter settings" name="settings_mode" type="select">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
128 <option value="default">default settings</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
129 <option value="change">change settings</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
130 </param>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
131
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
132 ## default settings
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
133
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
134 <when value="default">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
135 <param name="seedsize" type="hidden" value="20" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
136 <param name="slack" type="hidden" value="0.3" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
137 <param name="sp_min" type="hidden" value="100" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
138 <param name="sp_max" type="hidden" value="10000" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
139 <param name="maxdist" type="hidden" value="8" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
140 <param name="confdiff" type="hidden" value="2" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
141 <param name="confadpt" type="hidden" value="7" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
142
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
143 <param name="maxseeds" type="hidden" value="25" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
144 <param name="maxhits" type="hidden" value="250" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
145 <param name="clipping" type="hidden" value="++" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
146
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
147 <param name="selectivity" type="hidden" value="off" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
148 <param name="filter_output" type="hidden" value="off" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
149 <param name="sort" type="hidden" value="0" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
150 <param name="mmatch_notation" type="hidden" value="general" />
4
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
151 <param name="discard_overlapping_mates" type="hidden" value="" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
152 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
153
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
154 ## change settings
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
155
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
156 <when value="change">
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
157 <param help="Length of the seeds used in the reference genome hash table (SNAP index option -s)." label="seed size (default: 20)" name="seedsize" type="integer" value="20" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
158 <param help="Corresponds to the -h option of SNAP index." label="hash table slack size (default: 0.3)" name="slack" type="float" value="0.3" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
159
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
160 ## paired-end specific options
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
161 <param help="Corresponds to the first value of the SNAP option -s. Affects paired-end data only." label="minimum spacing to allow between paired ends (default: 100)" name="sp_min" type="integer" value="100" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
162 <param help="Corresponds to the second value of the SNAP option -s. Affects paired-end data only." label="maximum spacing to allow between paired ends (default: 10000)" name="sp_max" type="integer" value="10000" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
163 <param display="checkboxes" help="Consider overlapping mate pairs of the given orientation type(s) anomalous and discard them; allowed values: RF, FR, FF, RR; multiple types may be specified as a comma-separated list and ALL can be used as a shortcut for discarding all overlapping mate pairs; leave blank to retain all overlapping pairs. Affects paired-end data only." label="discard overlapping read pairs of type" multiple="true" name="discard_overlapping_mates" type="text" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
164 <param help="maximum edit distance allowed per read or pair (SNAP option -d); higher values allow more divergent alignments to be found, but increase the rate of misalignments." label="edit distance (default: 8)" name="maxdist" type="integer" value="8" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
165 <param help="Maximum hits to consider per seed (SNAP option -h); don't use a seed region in the alignment process if it matches more than maxhits regions in the reference genome. Higher values reduce the rate of misalignments, but reduce performance." label="maximum hits per seed (default: 250)" name="maxhits" type="integer" value="250" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
166 <param help="Confidence threshold (SNAP option -c); the minimum edit distance difference between two alternate alignments required to reject the poorer alignment as suboptimal; higher values increase the rate of ambiguously aligned reads." label="confidence threshold (default: 2)" name="confdiff" type="integer" value="2" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
167 <param help="Specifies how many seeds of a read may be ignored (based on the maximum hits value above) before the confidence threshold above gets increased by one for that read; helps fine-tuning alignment accuracy in repetitive regions of the genome." label="adaptive confdiff behaviour (default: 7)" name="confadpt" type="integer" value="7" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
168 <param help="Number of seeds to use per read (SNAP option -n) when trying to match it to the reference genome; higher numbers will increase the rate of aligned reads and reduce the rate of misalignments, but will reduce performance." label="maximum seeds per read (default: 25)" name="maxseeds" type="integer" value="25" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
169 <param help="Specifies from which end of a read low-quality bases should be clipped (SNAP option -Cxx)" label="read clipping (default: from back and front)" name="clipping" type="select">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
170 <option value="++">from back and front</option>
6
85214e4428fd upgrade to mimodd version 0.1.7.0
Wolfgang Maier
parents: 5
diff changeset
171 <option value="x+">from back only</option>
85214e4428fd upgrade to mimodd version 0.1.7.0
Wolfgang Maier
parents: 5
diff changeset
172 <option value="+x">from front only</option>
85214e4428fd upgrade to mimodd version 0.1.7.0
Wolfgang Maier
parents: 5
diff changeset
173 <option value="xx">no clipping</option>
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
174 </param>
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
175 <param help="randomly choose 1/selectivity of the reads to score (SNAP option -S). The tool uses the default of 1 (or a 0 setting) to indicate that all reads should be worked with." label="selectivity (default: 1)" name="selectivity" type="integer" value="1" />
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
176 <param help="filter output (SNAP option -F for certain classes of reads." label="filter output (default: no filtering)" name="filter_output" type="select">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
177 <option value="off">no filtering</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
178 <option value="a">aligned only</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
179 <option value="s">single-aligned only</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
180 <option value="u">unaligned only</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
181 </param>
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
182 <param help="Sort the output file by alignment location (SNAP option --so)." label="output sorting (default: sort by read coordinates)" name="sort" type="select">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
183 <option value="0">sort by read coordinates</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
184 <option value="off">no sorting</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
185 </param>
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
186 <param help="Indicates whether CIGAR strings in the generated SAM/BAM file should use M (alignment match) rather than = and X (sequence (mis-)match). Warning: Downstream variant calling based on samtools currently relies on the old-style M notation!!" label="CIGAR symbols for alignment matches/mismatches (default: M notation)" name="mmatch_notation" type="select">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
187 <option value="general">use M for both matches and mismatches</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
188 <option value="differentiate">use = for matches, X for mismatches</option>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
189 </param>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
190 </when>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
191 </conditional>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
192 </inputs>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
193
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
194 <outputs>
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
195 <data format="bam" label="Aligned reads from MiModd ${tool.name} on ${on_string}" name="outputfile">
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
196 <change_format>
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
197 <when format="sam" input="oformat" value="sam" />
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
198 </change_format>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
199 </data>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
200 </outputs>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
201
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
202 <help>
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
203 .. class:: infomark
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
204
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
205 **What it does**
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
206
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
207 The tool aligns the sequenced reads in an arbitrary number of input datasets against a common reference genome and stores the results in a single, possibly multi-sample output file. It supports a variety of different sequenced reads input formats, i.e., SAM, BAM, fastq and gzipped fastq, and both single-end and paired-end data.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
208
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
209 Internally, the tool uses the ultrafast, hashtable-based aligner SNAP (http://snap.cs.berkeley.edu), hence its name.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
210
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
211 **Notes:**
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
212
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
213 1) In its standard configuration Galaxy will decompress any .gz files during their upload, so the option to align gzipped fastq input is useful only with customized Galaxy instances or by using linked files as explained in our `recipe for using gzipped fastq files in Galaxy`_ from the `MiModD user guide`_.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
214
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
215 2) To use paired-end fastq data with the tool the read mate information needs to be split over two fastq files in corresponding order.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
216
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
217 **TIP:** If your paired-end data is arranged differently, you may look into the *fastq splitter* and *fastq de-interlacer* tools for Galaxy from the `Fastq Manipulation category`_ of the Galaxy Tool Shed to see if they can convert your files to the expected format.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
218
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
219 3) The tool supports the alignment of reads from the same sequencing run, but distributed across several input files.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
220
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
221 Generally, it expects the reads from each input dataset to belong to one read-group and will abort with an error message if any input dataset declares more than one read group or sample names in its header. Different datasets, however, are allowed to contain reads from the same read-group (as indicated by matching read-group IDs and sample names in their headers), in which case the reads will be combined into one group in the output.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
222
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
223 4) Read-group information is required for every input dataset!
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
224
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
225 We generally recommend to store NGS datasets in SAM/BAM format with run metadata stored in the file header. You can use the *NGS Run Annotation* and *Convert* tools to convert data in fastq format to SAM/BAM with added run information.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
226
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
227 While it is not our recommended approach, you can, if you prefer it, align reads from fastq files or SAM/BAM files without header read-group information. To do so, you **must** specify a SAM file that provides the missing information in its header along with the input dataset. You can generate a SAM header file with the *NGS Run Annotation* tool.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
228
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
229 Optionally, a SAM header file can also be used to replace existing read-group information in a headered SAM/BAM input file. This can be used to resolve read-group ID conflicts between multiple input files at tool runtime.
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
230
4
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
231 5) The options available under *further parameter settings* can have **big** effects on the alignment quality. You are strongly encouraged to consult the `tool documentation`_ for detailed explanations of the available options.
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
232
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
233 6) Currently, you cannot configure aligner-specific options separately for specific input files from within this Galaxy tool. If you need this advanced level of control, you should use the command line tool ``mimodd snap-batch``.
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
234
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
235 .. _Fastq Manipulation category: https://toolshed.g2.bx.psu.edu/repository/browse_repositories_in_category?id=310ff67d4caf6531
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
236 .. _recipe for using gzipped fastq files in Galaxy: http://mimodd.readthedocs.org/en/latest/recipes.html#use-gzipped-fastq-files-in-galaxy
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
237 .. _MiModD user guide: http://mimodd.readthedocs.org/en/latest
4
ffee8534a5c4 upgrade to mimodd version 0.1.6
Wolfgang Maier
parents: 2
diff changeset
238 .. _tool documentation: http://mimodd.readthedocs.org/en/latest/tool_doc.html#snap
0
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
239
6231ae8f87b8 Uploaded
wolma
parents:
diff changeset
240 </help>
23
5db0545b9004 update to v0.1.7.3
wolma
parents: 21
diff changeset
241 </tool>