annotate variant_recalibrator.xml @ 0:340633249b3d draft

Uploaded
author bgruening
date Mon, 02 Dec 2013 06:18:36 -0500
parents
children f244b8209eb8
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
340633249b3d Uploaded
bgruening
parents:
diff changeset
1 <tool id="gatk2_variant_recalibrator" name="Variant Recalibrator" version="0.0.7">
340633249b3d Uploaded
bgruening
parents:
diff changeset
2 <description></description>
340633249b3d Uploaded
bgruening
parents:
diff changeset
3 <expand macro="requirements" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
4 <macros>
340633249b3d Uploaded
bgruening
parents:
diff changeset
5 <import>gatk2_macros.xml</import>
340633249b3d Uploaded
bgruening
parents:
diff changeset
6 </macros>
340633249b3d Uploaded
bgruening
parents:
diff changeset
7 <command interpreter="python">
340633249b3d Uploaded
bgruening
parents:
diff changeset
8 gatk2_wrapper.py
340633249b3d Uploaded
bgruening
parents:
diff changeset
9 --stdout "${output_log}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
10 #for $var_count, $variant in enumerate( $reference_source.variants ):
340633249b3d Uploaded
bgruening
parents:
diff changeset
11 -d "--input:input_${var_count},%(file_type)s" "${variant.input_variants}" "${variant.input_variants.ext}" "input_variants_${var_count}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
12 #end for
340633249b3d Uploaded
bgruening
parents:
diff changeset
13 -p '
340633249b3d Uploaded
bgruening
parents:
diff changeset
14 @JAR_PATH@
340633249b3d Uploaded
bgruening
parents:
diff changeset
15 -T "VariantRecalibrator"
340633249b3d Uploaded
bgruening
parents:
diff changeset
16 \$GATK2_SITE_OPTIONS
340633249b3d Uploaded
bgruening
parents:
diff changeset
17
340633249b3d Uploaded
bgruening
parents:
diff changeset
18 @THREADS@
340633249b3d Uploaded
bgruening
parents:
diff changeset
19
340633249b3d Uploaded
bgruening
parents:
diff changeset
20 #if $reference_source.reference_source_selector != "history":
340633249b3d Uploaded
bgruening
parents:
diff changeset
21 -R "${reference_source.ref_file.fields.path}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
22 #end if
340633249b3d Uploaded
bgruening
parents:
diff changeset
23 --recal_file "${output_recal}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
24 --tranches_file "${output_tranches}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
25 --rscript_file "${output_rscript}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
26 '
340633249b3d Uploaded
bgruening
parents:
diff changeset
27
340633249b3d Uploaded
bgruening
parents:
diff changeset
28 #set $rod_binding_names = dict()
340633249b3d Uploaded
bgruening
parents:
diff changeset
29 #for $rod_binding in $rod_bind:
340633249b3d Uploaded
bgruening
parents:
diff changeset
30 #if str( $rod_binding.rod_bind_type.rod_bind_type_selector ) == 'custom':
340633249b3d Uploaded
bgruening
parents:
diff changeset
31 #set $rod_bind_name = $rod_binding.rod_bind_type.custom_rod_name
340633249b3d Uploaded
bgruening
parents:
diff changeset
32 #elif str( $rod_binding.rod_bind_type.rod_bind_type_selector ) == 'comp':
340633249b3d Uploaded
bgruening
parents:
diff changeset
33 #set $rod_bind_name = "comp" + $rod_binding.rod_bind_type.custom_rod_name
340633249b3d Uploaded
bgruening
parents:
diff changeset
34 #else
340633249b3d Uploaded
bgruening
parents:
diff changeset
35 #set $rod_bind_name = $rod_binding.rod_bind_type.rod_bind_type_selector
340633249b3d Uploaded
bgruening
parents:
diff changeset
36 #end if
340633249b3d Uploaded
bgruening
parents:
diff changeset
37 #set $rod_binding_names[$rod_bind_name] = $rod_binding_names.get( $rod_bind_name, -1 ) + 1
340633249b3d Uploaded
bgruening
parents:
diff changeset
38 #if $rod_binding.rod_bind_type.rod_training_type.rod_training_type_selector == "not_training_truth_known":
340633249b3d Uploaded
bgruening
parents:
diff changeset
39 -d "--resource:${rod_bind_name},%(file_type)s" "${rod_binding.rod_bind_type.input_rod}" "${rod_binding.rod_bind_type.input_rod.ext}" "input_${rod_bind_name}_${rod_binding_names[$rod_bind_name]}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
40 #else:
340633249b3d Uploaded
bgruening
parents:
diff changeset
41 -d "--resource:${rod_bind_name},%(file_type)s,known=${rod_binding.rod_bind_type.rod_training_type.known},training=${rod_binding.rod_bind_type.rod_training_type.training},truth=${rod_binding.rod_bind_type.rod_training_type.truth},bad=${rod_binding.rod_bind_type.rod_training_type.bad},prior=${rod_binding.rod_bind_type.rod_training_type.prior}" "${rod_binding.rod_bind_type.input_rod}" "${rod_binding.rod_bind_type.input_rod.ext}" "input_${rod_bind_name}_${rod_binding_names[$rod_bind_name]}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
42 #end if
340633249b3d Uploaded
bgruening
parents:
diff changeset
43 #end for
340633249b3d Uploaded
bgruening
parents:
diff changeset
44
340633249b3d Uploaded
bgruening
parents:
diff changeset
45 #include source=$standard_gatk_options#
340633249b3d Uploaded
bgruening
parents:
diff changeset
46
340633249b3d Uploaded
bgruening
parents:
diff changeset
47 ##start analysis specific options
340633249b3d Uploaded
bgruening
parents:
diff changeset
48 -p '
340633249b3d Uploaded
bgruening
parents:
diff changeset
49 #if str( $annotations ) != "None":
340633249b3d Uploaded
bgruening
parents:
diff changeset
50 #for $annotation in str( $annotations.fields.gatk_value ).split( ',' ):
340633249b3d Uploaded
bgruening
parents:
diff changeset
51 --use_annotation "${annotation}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
52 #end for
340633249b3d Uploaded
bgruening
parents:
diff changeset
53 #end if
340633249b3d Uploaded
bgruening
parents:
diff changeset
54 #for $additional_annotation in $additional_annotations:
340633249b3d Uploaded
bgruening
parents:
diff changeset
55 --use_annotation "${additional_annotation.additional_annotation_name}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
56 #end for
340633249b3d Uploaded
bgruening
parents:
diff changeset
57 --mode "${mode}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
58 '
340633249b3d Uploaded
bgruening
parents:
diff changeset
59
340633249b3d Uploaded
bgruening
parents:
diff changeset
60 #if $analysis_param_type.analysis_param_type_selector == "advanced":
340633249b3d Uploaded
bgruening
parents:
diff changeset
61 -p '
340633249b3d Uploaded
bgruening
parents:
diff changeset
62 --maxGaussians "${analysis_param_type.max_gaussians}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
63 --maxIterations "${analysis_param_type.max_iterations}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
64 --numKMeans "${analysis_param_type.num_k_means}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
65 --stdThreshold "${analysis_param_type.std_threshold}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
66 --qualThreshold "${analysis_param_type.qual_threshold}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
67 --shrinkage "${analysis_param_type.shrinkage}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
68 --dirichlet "${analysis_param_type.dirichlet}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
69 --priorCounts "${analysis_param_type.prior_counts}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
70 #if str( $analysis_param_type.bad_variant_selector.bad_variant_selector_type ) == 'percent':
340633249b3d Uploaded
bgruening
parents:
diff changeset
71 --percentBadVariants "${analysis_param_type.bad_variant_selector.percent_bad_variants}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
72 #else:
340633249b3d Uploaded
bgruening
parents:
diff changeset
73 --minNumBadVariants "${analysis_param_type.bad_variant_selector.min_num_bad_variants}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
74 #end if
340633249b3d Uploaded
bgruening
parents:
diff changeset
75 --target_titv "${analysis_param_type.target_titv}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
76 #for $tranche in [ $tranche.strip() for $tranche in str( $analysis_param_type.ts_tranche ).split( ',' ) if $tranche.strip() ]
340633249b3d Uploaded
bgruening
parents:
diff changeset
77 --TStranche "${tranche}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
78 #end for
340633249b3d Uploaded
bgruening
parents:
diff changeset
79 #for $ignore_filter in $analysis_param_type.ignore_filters:
340633249b3d Uploaded
bgruening
parents:
diff changeset
80 #set $ignore_filter_name = str( $ignore_filter.ignore_filter_type.ignore_filter_type_selector )
340633249b3d Uploaded
bgruening
parents:
diff changeset
81 #if $ignore_filter_name == "custom":
340633249b3d Uploaded
bgruening
parents:
diff changeset
82 #set $ignore_filter_name = str( $ignore_filter.ignore_filter_type.filter_name )
340633249b3d Uploaded
bgruening
parents:
diff changeset
83 #end if
340633249b3d Uploaded
bgruening
parents:
diff changeset
84 --ignore_filter "${ignore_filter_name}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
85 #end for
340633249b3d Uploaded
bgruening
parents:
diff changeset
86 --ts_filter_level "${analysis_param_type.ts_filter_level}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
87 '
340633249b3d Uploaded
bgruening
parents:
diff changeset
88 #end if
340633249b3d Uploaded
bgruening
parents:
diff changeset
89
340633249b3d Uploaded
bgruening
parents:
diff changeset
90
340633249b3d Uploaded
bgruening
parents:
diff changeset
91 &amp;&amp;
340633249b3d Uploaded
bgruening
parents:
diff changeset
92 mv "${output_rscript}.pdf" "${output_tranches_pdf}"
340633249b3d Uploaded
bgruening
parents:
diff changeset
93
340633249b3d Uploaded
bgruening
parents:
diff changeset
94 </command>
340633249b3d Uploaded
bgruening
parents:
diff changeset
95 <inputs>
340633249b3d Uploaded
bgruening
parents:
diff changeset
96 <conditional name="reference_source">
340633249b3d Uploaded
bgruening
parents:
diff changeset
97 <expand macro="reference_source_selector_param" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
98 <when value="cached">
340633249b3d Uploaded
bgruening
parents:
diff changeset
99 <repeat name="variants" title="Variant" min="1" help="-input,--input &amp;lt;input&amp;gt;">
340633249b3d Uploaded
bgruening
parents:
diff changeset
100 <param name="input_variants" type="data" format="vcf" label="Variant file to recalibrate" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
101 </repeat>
340633249b3d Uploaded
bgruening
parents:
diff changeset
102 <param name="ref_file" type="select" label="Using reference genome" help="-R,--reference_sequence &amp;lt;reference_sequence&amp;gt;">
340633249b3d Uploaded
bgruening
parents:
diff changeset
103 <options from_data_table="gatk_picard_indexes">
340633249b3d Uploaded
bgruening
parents:
diff changeset
104 <!-- <filter type="data_meta" key="dbkey" ref="variants[0].input_variants" column="dbkey"/> -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
105 </options>
340633249b3d Uploaded
bgruening
parents:
diff changeset
106 <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
107 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
108 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
109 <when value="history"> <!-- FIX ME!!!! -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
110 <repeat name="variants" title="Variant" min="1" help="-input,--input &amp;lt;input&amp;gt;">
340633249b3d Uploaded
bgruening
parents:
diff changeset
111 <param name="input_variants" type="data" format="vcf" label="Variant file to recalibrate" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
112 </repeat>
340633249b3d Uploaded
bgruening
parents:
diff changeset
113 <param name="ref_file" type="data" format="fasta" label="Using reference file" help="-R,--reference_sequence &amp;lt;reference_sequence&amp;gt;" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
114 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
115 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
116
340633249b3d Uploaded
bgruening
parents:
diff changeset
117 <repeat name="rod_bind" title="Binding for reference-ordered data" help="-resource,--resource &amp;lt;resource&amp;gt;">
340633249b3d Uploaded
bgruening
parents:
diff changeset
118 <conditional name="rod_bind_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
119 <param name="rod_bind_type_selector" type="select" label="Binding Type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
120 <option value="dbsnp" selected="True">dbSNP</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
121 <option value="variant">Variants</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
122 <option value="snps">SNPs</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
123 <option value="indels">INDELs</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
124 <option value="hapmap">HapMap</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
125 <option value="omni">OMNI</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
126 <option value="mask">Mask</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
127 <option value="custom">Custom</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
128 <option value="comp">Comp</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
129 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
130 <when value="variant">
340633249b3d Uploaded
bgruening
parents:
diff changeset
131 <param name="input_rod" type="data" format="vcf" label="Variant ROD file" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
132 <conditional name="rod_training_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
133 <param name="rod_training_type_selector" type="select" label="Use as training/truth/known sites">
340633249b3d Uploaded
bgruening
parents:
diff changeset
134 <option value="is_training_truth_known">Set training/truth/known sites</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
135 <option value="not_training_truth_known" selected="True">Don't Set options</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
136 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
137 <when value="not_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
138 <!-- do nothing here -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
139 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
140 <when value="is_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
141 <param name="known" type="boolean" label="Is Known Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
142 <param name="training" type="boolean" label="Is Training Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
143 <param name="truth" type="boolean" label="Is Truth Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
144 <param name="bad" type="boolean" label="Is Bad Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
145 <param name="prior" type="float" label="prior probability of being true" value="12.0"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
146 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
147 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
148 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
149 <when value="comp">
340633249b3d Uploaded
bgruening
parents:
diff changeset
150 <param name="input_rod" type="data" format="vcf" label="ROD file" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
151 <param name="custom_rod_name" type="text" value="Unnamed" label="ROD Name"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
152 <conditional name="rod_training_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
153 <param name="rod_training_type_selector" type="select" label="Use as training/truth/known sites">
340633249b3d Uploaded
bgruening
parents:
diff changeset
154 <option value="is_training_truth_known">Set training/truth/known sites</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
155 <option value="not_training_truth_known" selected="True">Don't Set options</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
156 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
157 <when value="not_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
158 <!-- do nothing here -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
159 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
160 <when value="is_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
161 <param name="known" type="boolean" label="Is Known Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
162 <param name="training" type="boolean" label="Is Training Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
163 <param name="truth" type="boolean" label="Is Truth Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
164 <param name="bad" type="boolean" label="Is Bad Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
165 <param name="prior" type="float" label="prior probability of being true" value="12.0"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
166 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
167 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
168 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
169 <when value="mask">
340633249b3d Uploaded
bgruening
parents:
diff changeset
170 <param name="input_rod" type="data" format="vcf" label="ROD file" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
171 <conditional name="rod_training_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
172 <param name="rod_training_type_selector" type="select" label="Use as training/truth/known sites">
340633249b3d Uploaded
bgruening
parents:
diff changeset
173 <option value="is_training_truth_known">Set training/truth/known sites</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
174 <option value="not_training_truth_known" selected="True">Don't Set options</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
175 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
176 <when value="not_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
177 <!-- do nothing here -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
178 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
179 <when value="is_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
180 <param name="known" type="boolean" label="Is Known Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
181 <param name="training" type="boolean" label="Is Training Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
182 <param name="truth" type="boolean" label="Is Truth Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
183 <param name="bad" type="boolean" label="Is Bad Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
184 <param name="prior" type="float" label="prior probability of being true" value="12.0"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
185 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
186 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
187 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
188 <when value="dbsnp">
340633249b3d Uploaded
bgruening
parents:
diff changeset
189 <param name="input_rod" type="data" format="vcf" label="ROD file" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
190 <conditional name="rod_training_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
191 <param name="rod_training_type_selector" type="select" label="Use as training/truth/known sites">
340633249b3d Uploaded
bgruening
parents:
diff changeset
192 <option value="is_training_truth_known">Set training/truth/known sites</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
193 <option value="not_training_truth_known" selected="True">Don't Set options</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
194 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
195 <when value="not_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
196 <!-- do nothing here -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
197 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
198 <when value="is_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
199 <param name="known" type="boolean" label="Is Known Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
200 <param name="training" type="boolean" label="Is Training Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
201 <param name="truth" type="boolean" label="Is Truth Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
202 <param name="bad" type="boolean" label="Is Bad Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
203 <param name="prior" type="float" label="prior probability of being true" value="12.0"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
204 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
205 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
206 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
207 <when value="snps">
340633249b3d Uploaded
bgruening
parents:
diff changeset
208 <param name="input_rod" type="data" format="vcf" label="ROD file" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
209 <conditional name="rod_training_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
210 <param name="rod_training_type_selector" type="select" label="Use as training/truth/known sites">
340633249b3d Uploaded
bgruening
parents:
diff changeset
211 <option value="is_training_truth_known">Set training/truth/known sites</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
212 <option value="not_training_truth_known" selected="True">Don't Set options</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
213 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
214 <when value="not_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
215 <!-- do nothing here -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
216 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
217 <when value="is_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
218 <param name="known" type="boolean" label="Is Known Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
219 <param name="training" type="boolean" label="Is Training Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
220 <param name="truth" type="boolean" label="Is Truth Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
221 <param name="bad" type="boolean" label="Is Bad Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
222 <param name="prior" type="float" label="prior probability of being true" value="12.0"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
223 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
224 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
225 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
226 <when value="hapmap">
340633249b3d Uploaded
bgruening
parents:
diff changeset
227 <param name="input_rod" type="data" format="vcf" label="ROD file" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
228 <conditional name="rod_training_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
229 <param name="rod_training_type_selector" type="select" label="Use as training/truth/known sites">
340633249b3d Uploaded
bgruening
parents:
diff changeset
230 <option value="is_training_truth_known">Set training/truth/known sites</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
231 <option value="not_training_truth_known" selected="True">Don't Set options</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
232 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
233 <when value="not_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
234 <!-- do nothing here -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
235 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
236 <when value="is_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
237 <param name="known" type="boolean" label="Is Known Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
238 <param name="training" type="boolean" label="Is Training Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
239 <param name="truth" type="boolean" label="Is Truth Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
240 <param name="bad" type="boolean" label="Is Bad Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
241 <param name="prior" type="float" label="prior probability of being true" value="12.0"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
242 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
243 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
244 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
245 <when value="omni">
340633249b3d Uploaded
bgruening
parents:
diff changeset
246 <param name="input_rod" type="data" format="vcf" label="ROD file" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
247 <conditional name="rod_training_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
248 <param name="rod_training_type_selector" type="select" label="Use as training/truth/known sites">
340633249b3d Uploaded
bgruening
parents:
diff changeset
249 <option value="is_training_truth_known">Set training/truth/known sites</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
250 <option value="not_training_truth_known" selected="True">Don't Set options</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
251 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
252 <when value="not_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
253 <!-- do nothing here -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
254 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
255 <when value="is_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
256 <param name="known" type="boolean" label="Is Known Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
257 <param name="training" type="boolean" label="Is Training Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
258 <param name="truth" type="boolean" label="Is Truth Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
259 <param name="bad" type="boolean" label="Is Bad Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
260 <param name="prior" type="float" label="prior probability of being true" value="12.0"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
261 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
262 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
263 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
264 <when value="indels">
340633249b3d Uploaded
bgruening
parents:
diff changeset
265 <param name="input_rod" type="data" format="vcf" label="ROD file" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
266 <conditional name="rod_training_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
267 <param name="rod_training_type_selector" type="select" label="Use as training/truth/known sites">
340633249b3d Uploaded
bgruening
parents:
diff changeset
268 <option value="is_training_truth_known">Set training/truth/known sites</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
269 <option value="not_training_truth_known" selected="True">Don't Set options</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
270 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
271 <when value="not_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
272 <!-- do nothing here -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
273 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
274 <when value="is_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
275 <param name="known" type="boolean" label="Is Known Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
276 <param name="training" type="boolean" label="Is Training Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
277 <param name="truth" type="boolean" label="Is Truth Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
278 <param name="bad" type="boolean" label="Is Bad Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
279 <param name="prior" type="float" label="prior probability of being true" value="12.0"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
280 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
281 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
282 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
283 <when value="custom">
340633249b3d Uploaded
bgruening
parents:
diff changeset
284 <param name="custom_rod_name" type="text" value="Unknown" label="ROD Name"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
285 <param name="input_rod" type="data" format="vcf" label="ROD file" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
286 <conditional name="rod_training_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
287 <param name="rod_training_type_selector" type="select" label="Use as training/truth/known sites">
340633249b3d Uploaded
bgruening
parents:
diff changeset
288 <option value="is_training_truth_known">Set training/truth/known sites</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
289 <option value="not_training_truth_known" selected="True">Don't Set options</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
290 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
291 <when value="not_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
292 <!-- do nothing here -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
293 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
294 <when value="is_training_truth_known">
340633249b3d Uploaded
bgruening
parents:
diff changeset
295 <param name="known" type="boolean" label="Is Known Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
296 <param name="training" type="boolean" label="Is Training Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
297 <param name="truth" type="boolean" label="Is Truth Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
298 <param name="bad" type="boolean" label="Is Bad Site" truevalue="true" falsevalue="false"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
299 <param name="prior" type="float" label="prior probability of being true" value="12.0"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
300 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
301 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
302 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
303 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
304 </repeat>
340633249b3d Uploaded
bgruening
parents:
diff changeset
305
340633249b3d Uploaded
bgruening
parents:
diff changeset
306 <param name="annotations" type="select" multiple="True" display="checkboxes" label="annotations which should used for calculations" help="-an,--use_annotation &amp;lt;use_annotation&amp;gt;">
340633249b3d Uploaded
bgruening
parents:
diff changeset
307 <!-- load the available annotations from an external configuration file, since additional ones can be added to local installs -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
308 <options from_data_table="gatk2_annotations">
340633249b3d Uploaded
bgruening
parents:
diff changeset
309 <filter type="multiple_splitter" column="tools_valid_for" separator=","/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
310 <filter type="static_value" value="VariantRecalibrator" column="tools_valid_for"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
311 </options>
340633249b3d Uploaded
bgruening
parents:
diff changeset
312 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
313
340633249b3d Uploaded
bgruening
parents:
diff changeset
314 <repeat name="additional_annotations" title="Additional annotation" help="-an,--use_annotation &amp;lt;use_annotation&amp;gt;">
340633249b3d Uploaded
bgruening
parents:
diff changeset
315 <param name="additional_annotation_name" type="text" value="" label="Annotation name" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
316 </repeat>
340633249b3d Uploaded
bgruening
parents:
diff changeset
317
340633249b3d Uploaded
bgruening
parents:
diff changeset
318 <param name="mode" type="select" label="Recalibration mode" help="-mode,--mode &amp;lt;mode&amp;gt;">
340633249b3d Uploaded
bgruening
parents:
diff changeset
319 <option value="SNP" selected="True">SNP</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
320 <option value="INDEL">INDEL</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
321 <option value="BOTH">BOTH</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
322 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
323
340633249b3d Uploaded
bgruening
parents:
diff changeset
324 <expand macro="gatk_param_type_conditional" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
325
340633249b3d Uploaded
bgruening
parents:
diff changeset
326 <expand macro="analysis_type_conditional">
340633249b3d Uploaded
bgruening
parents:
diff changeset
327 <param name="max_gaussians" type="integer" label="maximum number of Gaussians to try during variational Bayes Algorithm" value="10" help="-mG,--maxGaussians &amp;lt;maxGaussians&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
328 <param name="max_iterations" type="integer" label="maximum number of maximum number of VBEM iterations to be performed in variational Bayes Algorithm" value="100" help="-mI,--maxIterations &amp;lt;maxIterations&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
329 <param name="num_k_means" type="integer" label="number of k-means iterations to perform in order to initialize the means of the Gaussians in the Gaussian mixture model" value="30" help="-nKM,--numKMeans &amp;lt;numKMeans&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
330 <param name="std_threshold" type="float" label="If a variant has annotations more than -std standard deviations away from mean then don't use it for building the Gaussian mixture model." value="8.0" help="-std,--stdThreshold &amp;lt;stdThreshold&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
331 <param name="qual_threshold" type="float" label="If a known variant has raw QUAL value less than -qual then don't use it for building the Gaussian mixture model." value="80.0" help="-qual,--qualThreshold &amp;lt;qualThreshold&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
332 <param name="shrinkage" type="float" label="shrinkage parameter in variational Bayes algorithm" value="1.0" help="-shrinkage,--shrinkage &amp;lt;shrinkage&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
333 <param name="dirichlet" type="float" label="dirichlet parameter in variational Bayes algorithm" value="0.001" help="-dirichlet,--dirichlet &amp;lt;dirichlet&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
334 <param name="prior_counts" type="float" label="number of prior counts to use in variational Bayes algorithm" value="20.0" help="-priorCounts,--priorCounts &amp;lt;priorCounts&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
335 <conditional name="bad_variant_selector">
340633249b3d Uploaded
bgruening
parents:
diff changeset
336 <param name="bad_variant_selector_type" type="select" label="How to specify bad variants">
340633249b3d Uploaded
bgruening
parents:
diff changeset
337 <option value="percent" selected="True">Percent</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
338 <option value="min_num">Number</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
339 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
340 <when value="percent">
340633249b3d Uploaded
bgruening
parents:
diff changeset
341 <param name="percent_bad_variants" type="float" label="percentage of the worst scoring variants to use when building the Gaussian mixture model of bad variants. 0.07 means bottom 7 percent." value="0.03" help="-percentBad,--percentBadVariants &amp;lt;percentBadVariants&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
342 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
343 <when value="min_num">
340633249b3d Uploaded
bgruening
parents:
diff changeset
344 <param name="min_num_bad_variants" type="integer" label="minimum amount of worst scoring variants to use when building the Gaussian mixture model of bad variants. Will override -percentBad arugment if necessary" value="2000" help="-minNumBad,--minNumBadVariants &amp;lt;minNumBadVariants&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
345 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
346 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
347 <param name="target_titv" type="float" label="expected novel Ti/Tv ratio to use when calculating FDR tranches and for display on optimization curve output figures. (approx 2.15 for whole genome experiments). ONLY USED FOR PLOTTING PURPOSES!" value="2.15" help="-titv,--target_titv &amp;lt;target_titv&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
348 <param name="ts_tranche" type="text" label="levels of novel false discovery rate (FDR, implied by ti/tv) at which to slice the data. (in percent, that is 1.0 for 1 percent)" value="100.0, 99.9, 99.0, 90.0" help="-tranche,--TStranche &amp;lt;TStranche&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
349 <repeat name="ignore_filters" title="Ignore Filter" help="-ignoreFilter,--ignore_filter &amp;lt;ignore_filter&amp;gt;">
340633249b3d Uploaded
bgruening
parents:
diff changeset
350 <conditional name="ignore_filter_type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
351 <param name="ignore_filter_type_selector" type="select" label="Filter Type">
340633249b3d Uploaded
bgruening
parents:
diff changeset
352 <option value="HARD_TO_VALIDATE">HARD_TO_VALIDATE</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
353 <option value="LowQual" >LowQual</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
354 <option value="custom" selected="True">Other</option>
340633249b3d Uploaded
bgruening
parents:
diff changeset
355 </param>
340633249b3d Uploaded
bgruening
parents:
diff changeset
356 <when value="custom">
340633249b3d Uploaded
bgruening
parents:
diff changeset
357 <param name="filter_name" type="text" value="" label="Filter name"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
358 </when>
340633249b3d Uploaded
bgruening
parents:
diff changeset
359 <when value="HARD_TO_VALIDATE" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
360 <when value="LowQual" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
361 </conditional>
340633249b3d Uploaded
bgruening
parents:
diff changeset
362 </repeat>
340633249b3d Uploaded
bgruening
parents:
diff changeset
363 <param name="ts_filter_level" type="float" label="truth sensitivity level at which to start filtering, used here to indicate filtered variants in plots" value="99.0" help="-ts_filter_level,--ts_filter_level &amp;lt;ts_filter_level&amp;gt;"/>
340633249b3d Uploaded
bgruening
parents:
diff changeset
364 </expand>
340633249b3d Uploaded
bgruening
parents:
diff changeset
365 </inputs>
340633249b3d Uploaded
bgruening
parents:
diff changeset
366 <outputs>
340633249b3d Uploaded
bgruening
parents:
diff changeset
367 <data format="gatk_recal" name="output_recal" label="${tool.name} on ${on_string} (Recalibration File)" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
368 <data format="gatk_tranche" name="output_tranches" label="${tool.name} on ${on_string} (Tranches File)" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
369 <data format="txt" name="output_rscript" label="${tool.name} on ${on_string} (RScript File)" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
370 <data format="pdf" name="output_tranches_pdf" label="${tool.name} on ${on_string} (PDF File)" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
371 <data format="txt" name="output_log" label="${tool.name} on ${on_string} (log)" />
340633249b3d Uploaded
bgruening
parents:
diff changeset
372 </outputs>
340633249b3d Uploaded
bgruening
parents:
diff changeset
373 <tests>
340633249b3d Uploaded
bgruening
parents:
diff changeset
374 <!-- ADD TESTS -->
340633249b3d Uploaded
bgruening
parents:
diff changeset
375 </tests>
340633249b3d Uploaded
bgruening
parents:
diff changeset
376 <help>
340633249b3d Uploaded
bgruening
parents:
diff changeset
377 **What it does**
340633249b3d Uploaded
bgruening
parents:
diff changeset
378
340633249b3d Uploaded
bgruening
parents:
diff changeset
379 Takes variant calls as .vcf files, learns a Gaussian mixture model over the variant annotations and evaluates the variant -- assigning an informative lod score
340633249b3d Uploaded
bgruening
parents:
diff changeset
380
340633249b3d Uploaded
bgruening
parents:
diff changeset
381 For more information on using the VariantRecalibrator module, see this `tool specific page &lt;http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_walkers_variantrecalibration_VariantRecalibrator.html&gt;`_.
340633249b3d Uploaded
bgruening
parents:
diff changeset
382
340633249b3d Uploaded
bgruening
parents:
diff changeset
383 To learn about best practices for variant detection using GATK, see this `overview &lt;http://www.broadinstitute.org/gatk/guide/topic?name=best-practices&gt;`_.
340633249b3d Uploaded
bgruening
parents:
diff changeset
384
340633249b3d Uploaded
bgruening
parents:
diff changeset
385 If you encounter errors, please view the `GATK FAQ &lt;http://www.broadinstitute.org/gatk/guide/topic?name=faqs&gt;`_.
340633249b3d Uploaded
bgruening
parents:
diff changeset
386
340633249b3d Uploaded
bgruening
parents:
diff changeset
387 ------
340633249b3d Uploaded
bgruening
parents:
diff changeset
388
340633249b3d Uploaded
bgruening
parents:
diff changeset
389 **Inputs**
340633249b3d Uploaded
bgruening
parents:
diff changeset
390
340633249b3d Uploaded
bgruening
parents:
diff changeset
391 GenomeAnalysisTK: VariantRecalibrator accepts a variant input file.
340633249b3d Uploaded
bgruening
parents:
diff changeset
392
340633249b3d Uploaded
bgruening
parents:
diff changeset
393
340633249b3d Uploaded
bgruening
parents:
diff changeset
394 **Outputs**
340633249b3d Uploaded
bgruening
parents:
diff changeset
395
340633249b3d Uploaded
bgruening
parents:
diff changeset
396 The output is in VCF format.
340633249b3d Uploaded
bgruening
parents:
diff changeset
397
340633249b3d Uploaded
bgruening
parents:
diff changeset
398
340633249b3d Uploaded
bgruening
parents:
diff changeset
399 Go `here &lt;http://www.broadinstitute.org/gatk/guide/topic?name=intro&gt;`_ for details on GATK file formats.
340633249b3d Uploaded
bgruening
parents:
diff changeset
400
340633249b3d Uploaded
bgruening
parents:
diff changeset
401 -------
340633249b3d Uploaded
bgruening
parents:
diff changeset
402
340633249b3d Uploaded
bgruening
parents:
diff changeset
403 **Settings**::
340633249b3d Uploaded
bgruening
parents:
diff changeset
404
340633249b3d Uploaded
bgruening
parents:
diff changeset
405
340633249b3d Uploaded
bgruening
parents:
diff changeset
406 tranches_file The output tranches file used by ApplyRecalibration
340633249b3d Uploaded
bgruening
parents:
diff changeset
407 use_annotation The names of the annotations which should used for calculations
340633249b3d Uploaded
bgruening
parents:
diff changeset
408 mode Recalibration mode to employ: 1.) SNP for recalibrating only snps (emitting indels untouched in the output VCF); 2.) INDEL for indels; and 3.) BOTH for recalibrating both snps and indels simultaneously. (SNP|INDEL|BOTH)
340633249b3d Uploaded
bgruening
parents:
diff changeset
409 maxGaussians The maximum number of Gaussians to try during variational Bayes algorithm
340633249b3d Uploaded
bgruening
parents:
diff changeset
410 maxIterations The maximum number of VBEM iterations to be performed in variational Bayes algorithm. Procedure will normally end when convergence is detected.
340633249b3d Uploaded
bgruening
parents:
diff changeset
411 numKMeans The number of k-means iterations to perform in order to initialize the means of the Gaussians in the Gaussian mixture model.
340633249b3d Uploaded
bgruening
parents:
diff changeset
412 stdThreshold If a variant has annotations more than -std standard deviations away from mean then don't use it for building the Gaussian mixture model.
340633249b3d Uploaded
bgruening
parents:
diff changeset
413 qualThreshold If a known variant has raw QUAL value less than -qual then don't use it for building the Gaussian mixture model.
340633249b3d Uploaded
bgruening
parents:
diff changeset
414 shrinkage The shrinkage parameter in variational Bayes algorithm.
340633249b3d Uploaded
bgruening
parents:
diff changeset
415 dirichlet The dirichlet parameter in variational Bayes algorithm.
340633249b3d Uploaded
bgruening
parents:
diff changeset
416 priorCounts The number of prior counts to use in variational Bayes algorithm.
340633249b3d Uploaded
bgruening
parents:
diff changeset
417 percentBadVariants What percentage of the worst scoring variants to use when building the Gaussian mixture model of bad variants. 0.07 means bottom 7 percent.
340633249b3d Uploaded
bgruening
parents:
diff changeset
418 minNumBadVariants The minimum amount of worst scoring variants to use when building the Gaussian mixture model of bad variants. Will override -percentBad arugment if necessary.
340633249b3d Uploaded
bgruening
parents:
diff changeset
419 recal_file The output recal file used by ApplyRecalibration
340633249b3d Uploaded
bgruening
parents:
diff changeset
420 target_titv The expected novel Ti/Tv ratio to use when calculating FDR tranches and for display on optimization curve output figures. (approx 2.15 for whole genome experiments). ONLY USED FOR PLOTTING PURPOSES!
340633249b3d Uploaded
bgruening
parents:
diff changeset
421 TStranche The levels of novel false discovery rate (FDR, implied by ti/tv) at which to slice the data. (in percent, that is 1.0 for 1 percent)
340633249b3d Uploaded
bgruening
parents:
diff changeset
422 ignore_filter If specified the optimizer will use variants even if the specified filter name is marked in the input VCF file
340633249b3d Uploaded
bgruening
parents:
diff changeset
423 path_to_Rscript The path to your implementation of Rscript. For Broad users this is maybe /broad/tools/apps/R-2.6.0/bin/Rscript
340633249b3d Uploaded
bgruening
parents:
diff changeset
424 rscript_file The output rscript file generated by the VQSR to aid in visualization of the input data and learned model
340633249b3d Uploaded
bgruening
parents:
diff changeset
425 path_to_resources Path to resources folder holding the Sting R scripts.
340633249b3d Uploaded
bgruening
parents:
diff changeset
426 ts_filter_level The truth sensitivity level at which to start filtering, used here to indicate filtered variants in plots
340633249b3d Uploaded
bgruening
parents:
diff changeset
427
340633249b3d Uploaded
bgruening
parents:
diff changeset
428 @CITATION_SECTION@
340633249b3d Uploaded
bgruening
parents:
diff changeset
429 </help>
340633249b3d Uploaded
bgruening
parents:
diff changeset
430 </tool>