annotate freebayes.xml @ 13:2dd40b601766 draft

Uploaded
author devteam
date Fri, 20 Jun 2014 17:18:49 -0400
parents
children 59731e950e23
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
13
2dd40b601766 Uploaded
devteam
parents:
diff changeset
1 <?xml version="1.0"?>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
2 <tool id="freebayes" name="FreeBayes" version="freebayes-0.9.14">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
3 <requirements>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
4 <requirement type="package" version="freebayes-0.9.14_8a407cf5f4">freebayes</requirement>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
5 <requirement type="package" version="0.1.18">samtools</requirement>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
6 </requirements>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
7 <description> - Bayesian genetic variant detector</description>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
8 <command>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
9 ##set up input files
2dd40b601766 Uploaded
devteam
parents:
diff changeset
10
2dd40b601766 Uploaded
devteam
parents:
diff changeset
11 #set $reference_fasta_filename = "localref.fa"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
12
2dd40b601766 Uploaded
devteam
parents:
diff changeset
13 #if str( $reference_source.reference_source_selector ) == "history":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
14 ln -s "${reference_source.ref_file}" "${reference_fasta_filename}" &amp;&amp;
2dd40b601766 Uploaded
devteam
parents:
diff changeset
15 samtools faidx "${reference_fasta_filename}" 2&gt;&amp;1 || echo "Error running samtools faidx for FreeBayes" &gt;&amp;2 &amp;&amp;
2dd40b601766 Uploaded
devteam
parents:
diff changeset
16 #else:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
17 #set $reference_fasta_filename = str( $reference_source.ref_file.fields.path )
2dd40b601766 Uploaded
devteam
parents:
diff changeset
18 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
19
2dd40b601766 Uploaded
devteam
parents:
diff changeset
20 #for $bam_count, $input_bam in enumerate( $reference_source.input_bams ):
2dd40b601766 Uploaded
devteam
parents:
diff changeset
21 ln -s "${input_bam.input_bam}" "localbam_${bam_count}.bam" &amp;&amp;
2dd40b601766 Uploaded
devteam
parents:
diff changeset
22 ln -s "${input_bam.input_bam.metadata.bam_index}" "localbam_${bam_count}.bam.bai" &amp;&amp;
2dd40b601766 Uploaded
devteam
parents:
diff changeset
23 #end for
2dd40b601766 Uploaded
devteam
parents:
diff changeset
24
2dd40b601766 Uploaded
devteam
parents:
diff changeset
25 ## Tabixize optional input_varinat_vcf file (for --variant-input option)
2dd40b601766 Uploaded
devteam
parents:
diff changeset
26
2dd40b601766 Uploaded
devteam
parents:
diff changeset
27 #if ( str( $options_type.options_type_selector ) == 'cline' or str( $options_type.options_type_selector ) == 'full' ) and str( $options_type.optional_inputs.optional_inputs_selector ) == 'set' and str( $options_type.optional_inputs.input_variant_type.input_variant_type_selector ) == "provide_vcf":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
28 ln -s "${options_type.optional_inputs.input_variant_type.input_variant_vcf}" input_variant_vcf.vcf.gz &amp;&amp;
2dd40b601766 Uploaded
devteam
parents:
diff changeset
29 ln -s "${Tabixized_input}" input_variant_vcf.vcf.gz.tbi &amp;&amp;
2dd40b601766 Uploaded
devteam
parents:
diff changeset
30 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
31
2dd40b601766 Uploaded
devteam
parents:
diff changeset
32 ##finished setting up inputs
2dd40b601766 Uploaded
devteam
parents:
diff changeset
33
2dd40b601766 Uploaded
devteam
parents:
diff changeset
34 ##COMMAND LINE STARTS HERE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
35
2dd40b601766 Uploaded
devteam
parents:
diff changeset
36 freebayes
2dd40b601766 Uploaded
devteam
parents:
diff changeset
37 #for $bam_count, $input_bam in enumerate( $reference_source.input_bams ):
2dd40b601766 Uploaded
devteam
parents:
diff changeset
38 --bam "localbam_${bam_count}.bam"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
39 #end for
2dd40b601766 Uploaded
devteam
parents:
diff changeset
40 --fasta-reference "${reference_fasta_filename}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
41
2dd40b601766 Uploaded
devteam
parents:
diff changeset
42 ##outputs
2dd40b601766 Uploaded
devteam
parents:
diff changeset
43 --vcf "${output_vcf}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
44
2dd40b601766 Uploaded
devteam
parents:
diff changeset
45 #if str( $target_limit_type.target_limit_type_selector ) == "limit_by_target_file":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
46 --targets "${target_limit_type.input_target_bed}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
47 #elif str( $target_limit_type.target_limit_type_selector ) == "limit_by_region":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
48 --region "${target_limit_type.region_chromosome}:${target_limit_type.region_start}..${target_limit_type.region_end}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
49 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
50
2dd40b601766 Uploaded
devteam
parents:
diff changeset
51 ##advanced options
2dd40b601766 Uploaded
devteam
parents:
diff changeset
52 #if str( $options_type.options_type_selector ) == "simple":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
53 ##do nothing as command like build up to this point is sufficinet for simple diploid calling
2dd40b601766 Uploaded
devteam
parents:
diff changeset
54
2dd40b601766 Uploaded
devteam
parents:
diff changeset
55 #elif str( $options_type.options_type_selector ) == "simple_w_filters":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
56
2dd40b601766 Uploaded
devteam
parents:
diff changeset
57 --standard-filters
2dd40b601766 Uploaded
devteam
parents:
diff changeset
58 --min-coverage "${options_type.min_coverage}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
59
2dd40b601766 Uploaded
devteam
parents:
diff changeset
60 #elif str( $options_type.options_type_selector ) == "naive":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
61
2dd40b601766 Uploaded
devteam
parents:
diff changeset
62 --haplotype-length 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
63 --min-alternate-count 1
2dd40b601766 Uploaded
devteam
parents:
diff changeset
64 --min-alternate-fraction 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
65 --pooled-continuous
2dd40b601766 Uploaded
devteam
parents:
diff changeset
66 --report-monomorphic
2dd40b601766 Uploaded
devteam
parents:
diff changeset
67
2dd40b601766 Uploaded
devteam
parents:
diff changeset
68 #elif str( $options_type.options_type_selector ) == "naive_w_filters":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
69
2dd40b601766 Uploaded
devteam
parents:
diff changeset
70 --haplotype-length 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
71 --min-alternate-count 1
2dd40b601766 Uploaded
devteam
parents:
diff changeset
72 --min-alternate-fraction 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
73 --pooled-continuous
2dd40b601766 Uploaded
devteam
parents:
diff changeset
74 --report-monomorphic
2dd40b601766 Uploaded
devteam
parents:
diff changeset
75 --standard-filters
2dd40b601766 Uploaded
devteam
parents:
diff changeset
76 --min-coverage "${options_type.min_coverage}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
77
2dd40b601766 Uploaded
devteam
parents:
diff changeset
78 #elif str( $options_type.options_type_selector ) == "cline":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
79
2dd40b601766 Uploaded
devteam
parents:
diff changeset
80 ${options_type.cline}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
81
2dd40b601766 Uploaded
devteam
parents:
diff changeset
82 @optional_inputs_outputs@
2dd40b601766 Uploaded
devteam
parents:
diff changeset
83
2dd40b601766 Uploaded
devteam
parents:
diff changeset
84 #elif str( $options_type.options_type_selector ) == "full":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
85
2dd40b601766 Uploaded
devteam
parents:
diff changeset
86 ##optional inputs and outputs
2dd40b601766 Uploaded
devteam
parents:
diff changeset
87
2dd40b601766 Uploaded
devteam
parents:
diff changeset
88 @optional_inputs_outputs@
2dd40b601766 Uploaded
devteam
parents:
diff changeset
89
2dd40b601766 Uploaded
devteam
parents:
diff changeset
90 ## REPORTING
2dd40b601766 Uploaded
devteam
parents:
diff changeset
91
2dd40b601766 Uploaded
devteam
parents:
diff changeset
92 #if str( $options_type.reporting.reporting_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
93 --pvar ${options_type.reporting.pvar}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
94 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
95
2dd40b601766 Uploaded
devteam
parents:
diff changeset
96 ## POPULATION MODEL
2dd40b601766 Uploaded
devteam
parents:
diff changeset
97
2dd40b601766 Uploaded
devteam
parents:
diff changeset
98 #if str( $options_type.population_model.population_model_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
99 --theta "${options_type.population_model.T}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
100 --ploidy "${options_type.population_model.P}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
101 ${options_type.population_model.J}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
102 ${options_type.population_model.K}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
103
2dd40b601766 Uploaded
devteam
parents:
diff changeset
104 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
105
2dd40b601766 Uploaded
devteam
parents:
diff changeset
106 ## REFERENCE ALLELE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
107
2dd40b601766 Uploaded
devteam
parents:
diff changeset
108 #if str( $options_type.reference_allele.reference_allele_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
109 ${options_type.reference_allele.Z}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
110 --reference-quality "${options_type.reference_allele.reference_quality}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
111 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
112
2dd40b601766 Uploaded
devteam
parents:
diff changeset
113 ## ALLELE SCOPE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
114
2dd40b601766 Uploaded
devteam
parents:
diff changeset
115 #if str( $options_type.allele_scope.allele_scope_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
116 ${options_type.allele_scope.I}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
117 ${options_type.allele_scope.i}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
118 ${options_type.allele_scope.X}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
119 ${options_type.allele_scope.u}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
120 -n "${options_type.allele_scope.n}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
121 --haplotype-length "${options_type.allele_scope.haplotype_length}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
122 --min-repeat-length "${options_type.allele_scope.min_repeat_length}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
123 --min-repeat-entropy "${options_type.allele_scope.min_repeat_entropy}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
124 ${options_type.allele_scope.no_partial_observations}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
125 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
126
2dd40b601766 Uploaded
devteam
parents:
diff changeset
127 ## REALIGNMENT
2dd40b601766 Uploaded
devteam
parents:
diff changeset
128
2dd40b601766 Uploaded
devteam
parents:
diff changeset
129 ${options_type.O}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
130
2dd40b601766 Uploaded
devteam
parents:
diff changeset
131 ##INPUT FILTERS
2dd40b601766 Uploaded
devteam
parents:
diff changeset
132
2dd40b601766 Uploaded
devteam
parents:
diff changeset
133 #if str( $options_type.input_filters.input_filters_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
134 ${options_type.input_filters.use_duplicate_reads}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
135 -m "${options_type.input_filters.m}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
136 -q "${options_type.input_filters.q}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
137 -R "${options_type.input_filters.R}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
138 -Y "${options_type.input_filters.Y}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
139
2dd40b601766 Uploaded
devteam
parents:
diff changeset
140 #if str( $options_type.input_filters.mismatch_filters.mismatch_filters_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
141 -Q "${options_type.input_filters.mismatch_filters.Q}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
142 -U "${options_type.input_filters.mismatch_filters.U}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
143 -z "${options_type.input_filters.mismatch_filters.z}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
144 --read-snp-limit "${options_type.input_filters.mismatch_filters.read_snp_limit}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
145 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
146
2dd40b601766 Uploaded
devteam
parents:
diff changeset
147 -e "${options_type.input_filters.e}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
148 -F "${options_type.input_filters.F}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
149 -C "${options_type.input_filters.C}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
150 --min-alternate-qsum "${options_type.input_filters.min_alternate_qsum}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
151 -G "${options_type.input_filters.G}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
152 --min-coverage "${options_type.input_filters.min_coverage}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
153 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
154
2dd40b601766 Uploaded
devteam
parents:
diff changeset
155 ## POPULATION AND MAPPABILITY PRIORS
2dd40b601766 Uploaded
devteam
parents:
diff changeset
156
2dd40b601766 Uploaded
devteam
parents:
diff changeset
157 #if str( $options_type.population_mappability_priors.population_mappability_priors_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
158 ${options_type.population_mappability_priors.k}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
159 ${options_type.population_mappability_priors.w}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
160 ${options_type.population_mappability_priors.V}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
161 ${options_type.population_mappability_priors.a}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
162 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
163
2dd40b601766 Uploaded
devteam
parents:
diff changeset
164 ## GENOTYPE LIKELIHOODS
2dd40b601766 Uploaded
devteam
parents:
diff changeset
165
2dd40b601766 Uploaded
devteam
parents:
diff changeset
166 #if str( $options_type.genotype_likelihoods.genotype_likelihoods_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
167 --base-quality-cap "${$options_type.genotype_likelihoods.base_quality_cap}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
168 ${$options_type.genotype_likelihoods.experimental_gls}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
169 --prob_contamination "${$options_type.genotype_likelihoods.prob_contamination}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
170 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
171
2dd40b601766 Uploaded
devteam
parents:
diff changeset
172 ## ALGORITHMIC FEATURES
2dd40b601766 Uploaded
devteam
parents:
diff changeset
173
2dd40b601766 Uploaded
devteam
parents:
diff changeset
174 #if str( $options_type.algorithmic_features.algorithmic_features_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
175 v
2dd40b601766 Uploaded
devteam
parents:
diff changeset
176 -B "${options_type.algorithmic_features.B}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
177 --genotyping-max-banddepth "${options_type.algorithmic_features.genotyping_max_banddepth}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
178 -W "${options_type.algorithmic_features.W}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
179 ${options_type.algorithmic_features.N}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
180
2dd40b601766 Uploaded
devteam
parents:
diff changeset
181 #if str( $options_type.algorithmic_features.genotype_variant_threshold.genotype_variant_threshold_selector ) == "True":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
182 -S "${options_type.algorithmic_features.genotype_variant_threshold.S}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
183 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
184
2dd40b601766 Uploaded
devteam
parents:
diff changeset
185 ${options_type.algorithmic_features.j}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
186 ${options_type.algorithmic_features.H}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
187 -D "${options_type.algorithmic_features.D}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
188 ${options_type.algorithmic_features.genotype_qualities}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
189 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
190 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
191
2dd40b601766 Uploaded
devteam
parents:
diff changeset
192 </command>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
193
2dd40b601766 Uploaded
devteam
parents:
diff changeset
194 <macros>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
195 <token name="@optional_inputs_outputs@">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
196 ## This token gets injected in commane in two instances: when options_type.options_type_selector == "full" and "cline"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
197
2dd40b601766 Uploaded
devteam
parents:
diff changeset
198 #if $options_type.optional_inputs.optional_inputs_selector:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
199
2dd40b601766 Uploaded
devteam
parents:
diff changeset
200 #if $options_type.optional_inputs.output_trace_option:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
201 --trace "${output_trace}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
202 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
203
2dd40b601766 Uploaded
devteam
parents:
diff changeset
204 #if $options_type.optional_inputs.output_failed_alleles_option:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
205 --failed-alleles "${output_failed_alleles_bed}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
206 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
207
2dd40b601766 Uploaded
devteam
parents:
diff changeset
208 #if $options_type.optional_inputs.samples:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
209 --samples "${options_type.optional_inputs.samples}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
210 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
211
2dd40b601766 Uploaded
devteam
parents:
diff changeset
212 #if $options_type.optional_inputs.populations:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
213 --populations "${options_type.optional_inputs.populations}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
214 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
215
2dd40b601766 Uploaded
devteam
parents:
diff changeset
216 #if $options_type.optional_inputs.A:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
217 --cnv-map "${options_type.optional_inputs.A}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
218 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
219
2dd40b601766 Uploaded
devteam
parents:
diff changeset
220 #if str( $options_type.optional_inputs.input_variant_type.input_variant_type_selector ) == "provide_vcf":
2dd40b601766 Uploaded
devteam
parents:
diff changeset
221 --variant-input input_variant_vcf.vcf.gz ## input_variant_vcf.vcf.gz is symlinked to a galaxy-generated dataset in "Tabixize optional input_varinat_vcf file" section of the command line above
2dd40b601766 Uploaded
devteam
parents:
diff changeset
222 ${options_type.optional_inputs.input_variant_type.only_use_input_alleles}
2dd40b601766 Uploaded
devteam
parents:
diff changeset
223 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
224
2dd40b601766 Uploaded
devteam
parents:
diff changeset
225 #if $options_type.optional_inputs.haplotype_basis_alleles:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
226 --haplotype-basis-alleles "${options_type.optional_inputs.haplotype_basis_alleles}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
227 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
228
2dd40b601766 Uploaded
devteam
parents:
diff changeset
229 #if $options_type.optional_inputs.observation_bias:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
230 --observation-bias "${options_type.optional_inputs.observation_bias}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
231 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
232
2dd40b601766 Uploaded
devteam
parents:
diff changeset
233 #if $options_type.optional_inputs.contamination_estimates:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
234 --contamination-estimates "${options_type.optional_inputs.contamination_estimates}"
2dd40b601766 Uploaded
devteam
parents:
diff changeset
235 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
236
2dd40b601766 Uploaded
devteam
parents:
diff changeset
237 #end if
2dd40b601766 Uploaded
devteam
parents:
diff changeset
238 </token>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
239 <xml name="optional_file_inputs">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
240 <conditional name="optional_inputs">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
241 <param name="optional_inputs_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Do you want to provide additional inputs?" help="Sets --samples, --populations, --cnv-map, --trace, --failed-alleles, --varinat-input, --only-use-input-alleles, --haplotype-basis-alleles, --report-all-haplotype-alleles, --report-monomorphic options, --observation-bias, and --contamination-estimates" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
242 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
243 <param name="output_failed_alleles_option" type="boolean" truevalue="--failed-alleles" falsevalue="" checked="False" label="Write out failed alleles file" help="--failed-alleles" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
244 <param name="output_trace_option" type="boolean" truevalue="--trace" falsevalue="" checked="False" label="Write out algorithm trace file" help="--trace"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
245 <param name="samples" type="data" format="txt" label="Limit analysis to samples listed (one per line) in the FILE" optional="True" help="-s --samples; default=By default FreeBayes will analyze all samples in its input BAM files"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
246 <param name="populations" type="data" format="txt" label="Populations File" optional="True" help="--populations; default=False. Each line of FILE should list a sample and a population which it is part of. The population-based bayesian inference model will then be partitioned on the basis of the populations" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
247 <param name="A" type="data" format="bed" label="Read a copy number map from the BED file FILE" optional="True" help="-A --cnv-map; default=copy number is set to as specified by --ploidy. Read a copy number map from the BED file FILE, which has the format: reference sequence, start, end, sample name, copy number ... for each region in each sample which does not have the default copy number as set by --ploidy."/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
248 <conditional name="input_variant_type">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
249 <param name="input_variant_type_selector" type="select" label="Provide variants file">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
250 <option value="do_not_provide" selected="True">Do not provide</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
251 <option value="provide_vcf">Provide VCF file</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
252 </param>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
253 <when value="do_not_provide">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
254 <!-- Do nothing here -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
255 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
256 <when value="provide_vcf">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
257 <param name="input_variant_vcf" type="data" format="vcf_bgzip" label="Use variants reported in VCF file as input to the algorithm">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
258 <conversion name="Tabixized_input" type="tabix" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
259 </param>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
260 <param name="only_use_input_alleles" type="boolean" truevalue="--only-use-input-alleles" falsevalue="" checked="False" label="Only provide variant calls and genotype likelihoods for sites in VCF" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
261 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
262 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
263 <param name="haplotype_basis_alleles" type="data" format="vcf" label="Only use variant alleles provided in this input VCF for the construction of complex or haplotype alleles" optional="True" help="--haplotype-basis-alleles" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
264 <param name="report_monomorphic" type="boolean" truevalue="--report-monomorphic" falsevalue="" checked="False" label="Report even loci which appear to be monomorphic, and report all considered alleles, even those which are not in called genotypes." help="--report-monomorphic " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
265 <param name="observation_bias" optional="True" type="data" format="tabular" label="Load read length-dependent allele observation biases from" help="--observation-bias; The format is [length] [alignment efficiency relative to reference] where the efficiency is 1 if there is no relative observation bias" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
266 <param name="contamination_estimates" optional="True" type="data" format="tabular" label="Upload per-sample estimates of contamination from" help="--contamination-estimates; The format should be: sample p(read=R|genotype=AR) p(read=A|genotype=AA) Sample '*' can be used to set default contamination estimates." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
267 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
268 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
269 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
270 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
271 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
272 </xml>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
273 </macros>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
274
2dd40b601766 Uploaded
devteam
parents:
diff changeset
275 <inputs>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
276 <conditional name="reference_source">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
277 <param name="reference_source_selector" type="select" label="Load reference genome from">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
278 <option value="cached">Local cache</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
279 <option value="history">History</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
280 </param>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
281 <when value="cached">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
282 <repeat name="input_bams" title="Sample BAM file" min="1">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
283 <param name="input_bam" type="data" format="bam" label="BAM file">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
284 <validator type="unspecified_build" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
285 <validator type="dataset_metadata_in_data_table" table_name="sam_fa_indexes" metadata_name="dbkey" metadata_column="1" message="Sequences are not currently available for the specified build." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
286 </param>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
287 </repeat>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
288 <param name="ref_file" type="select" label="Using reference genome">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
289 <options from_data_table="sam_fa_indexes">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
290 <!-- <filter type="sam_fa_indexes" key="dbkey" ref="input_bam" column="value"/> does not yet work in a repeat...-->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
291 </options>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
292 <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
293 </param>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
294 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
295 <when value="history"> <!-- FIX ME!!!! -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
296 <repeat name="input_bams" title="Sample BAM file" min="1">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
297 <param name="input_bam" type="data" format="bam" label="BAM file" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
298 </repeat>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
299 <param name="ref_file" type="data" format="fasta" label="Use the folloing dataset as the reference sequence" help="You can upload a FASTA sequence to the history and use it as reference" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
300 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
301 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
302
2dd40b601766 Uploaded
devteam
parents:
diff changeset
303 <conditional name="target_limit_type">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
304 <param name="target_limit_type_selector" type="select" label="Limit variant calling to a set of regions?" help="Sets --targets or --region options">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
305 <option value="do_not_limit" selected="True">Do not limit</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
306 <option value="limit_by_target_file">Limit by target file</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
307 <option value="limit_by_region">Limit to region</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
308 </param>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
309 <when value="do_not_limit">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
310 <!-- Do nothing here -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
311 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
312 <when value="limit_by_target_file">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
313 <param name="input_target_bed" type="data" format="bed" label="Limit analysis to targets listed in the BED-format FILE." help="-t --targets"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
314 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
315 <when value="limit_by_region">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
316 <param name="region_chromosome" type="text" label="Region Chromosome" value="" help="-r --region"/> <!--only once? -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
317 <param name="region_start" type="integer" label="Region Start" value="" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
318 <param name="region_end" type="integer" label="Region End" value="" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
319 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
320 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
321
2dd40b601766 Uploaded
devteam
parents:
diff changeset
322 <conditional name="options_type">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
323 <param name="options_type_selector" type="select" label="Choose parameter selection level" help="Select how much control over the freebayes run you need" >
2dd40b601766 Uploaded
devteam
parents:
diff changeset
324 <option value="simple" selected="True">1:Simple diploid calling</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
325 <option value="simple_w_filters">2:Simple diploid calling with filtering and coverage</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
326 <option value="naive">3:Frequency-based pooled calling</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
327 <option value="naive_w_filters">4:Frequency-based pooled calling with filtering and coverage</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
328 <option value="full">5:Complete list of all options</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
329 <option value="cline">6:Input parameters on the command line</option>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
330 </param>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
331 <when value="full">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
332
2dd40b601766 Uploaded
devteam
parents:
diff changeset
333 <expand macro="optional_file_inputs" /> <!-- see macros section -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
334
2dd40b601766 Uploaded
devteam
parents:
diff changeset
335 <!-- reporting -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
336
2dd40b601766 Uploaded
devteam
parents:
diff changeset
337 <conditional name="reporting">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
338 <param name="reporting_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Set reporting option?" help="Sets -P --pvar option" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
339 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
340 <param name="pvar" type="float" value="0.0" label="Report sites if the probability that there is a polymorphism at the site is greater than" help="-P --pvar; default=0.0. Note that post-filtering is generally recommended over the use of this parameter. " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
341 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
342 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
343 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
344 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
345 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
346
2dd40b601766 Uploaded
devteam
parents:
diff changeset
347 <!-- population model -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
348
2dd40b601766 Uploaded
devteam
parents:
diff changeset
349 <conditional name="population_model">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
350 <param name="population_model_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Set population model?" help="Sets --theta, --ploidy, --pooled-discrete, and --pooled-continuous options " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
351 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
352 <param name="T" type="float" value="0.001" label="The expected mutation rate or pairwise nucleotide diversity among the population under analysis" help="-T --theta; default = 0.001. This serves as the single parameter to the Ewens Sampling Formula prior model." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
353 <param name="P" type="integer" value="2" label="Set ploidy for the analysis" help="-p --ploidy; default=2" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
354 <param name="J" type="boolean" truevalue="-J" falsevalue="" checked="False" label="Assume that samples result from pooled sequencing" help="-J --pooled-discrete; default=False. Model pooled samples using discrete genotypes across pools. When using this flag, set --ploidy to the number of alleles in each sample or use the --cnv-map to define per-sample ploidy." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
355 <param name="K" type="boolean" truevalue="-K" falsevalue="" checked="False" label="Output all alleles which pass input filters, regardles of genotyping outcome or model" help="-K, --poled-continuous; default=False. " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
356 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
357 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
358 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
359 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
360 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
361
2dd40b601766 Uploaded
devteam
parents:
diff changeset
362 <!-- reference allele -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
363
2dd40b601766 Uploaded
devteam
parents:
diff changeset
364 <conditional name="reference_allele">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
365 <param name="reference_allele_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Use reference allele?" help="Sets --use-reference-allele and --reference-quality options " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
366 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
367 <param name="Z" type="boolean" truevalue="-Z" falsevalue="" checked="False" label="Include the reference allele in the analysis as if it is another sample from the same population" help="-Z --use-reference-allele; default=False" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
368 <param name="reference_quality" type="text" size="8" value="100,60" label="Assign mapping quality of MQ (100) to the reference allele at each site and base quality of BQ (60)" help="--reference-quality; default=100,60 " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
369 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
370 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
371 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
372 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
373 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
374
2dd40b601766 Uploaded
devteam
parents:
diff changeset
375 <!-- allelic scope -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
376
2dd40b601766 Uploaded
devteam
parents:
diff changeset
377 <conditional name="allele_scope">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
378 <param name="allele_scope_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Set allelic scope?" help="Sets -I, i, -X, -u, -n, -E, --haplotype-length, --min-repeat-length, --min-repeat-entropy, and --no-partial-observations options " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
379 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
380 <param name="I" type="boolean" truevalue="-I" falsevalue="" checked="False" label="Ignore SNP alleles" help="-I --no-snps; default=False" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
381 <param name="i" type="boolean" truevalue="-i" falsevalue="" checked="False" label="Ignore indels alleles" help="-i --no-indels; default=False" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
382 <param name="X" type="boolean" truevalue="-X" falsevalue="" checked="False" label="Ignore multi-nucleotide polymorphisms, MNPs" help="-X --no-mnps; default=False" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
383 <param name="u" type="boolean" truevalue="-u" falsevalue="" checked="False" label="Ignore complex events (composites of other classes)." help="-u --no-complex; default=False" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
384 <param name="n" type="integer" value="0" label="How many best SNP alleles to evaluate" help="-n --use-best-n-alleles; default=0 (all). Alleles are ranked by the sum of supporting quality scores. Set to 0 to evaluate all" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
385 <param name="haplotype_length" type="integer" value="3" label="Allow haplotype calls with contiguous embedded matches of up to (nucleotides)" help="--haplotype-length; default=3." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
386 <param name="min_repeat_length" type="integer" value="5" label="When assembling observations across repeats, require the total repeat length at least this many bp" help="--min-repeat-length; default=5." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
387 <param name="min_repeat_entropy" type="integer" value="0" label="To detect interrupted repeats, build across sequence until it has entropy > (bits per bp)" help="--min-repeat-entrpy; default=0 (off)." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
388 <param name="no_partial_observations" type="boolean" truevalue="--no-partial-observations" falsevalue="" checked="False" label="Exclude observations which do not fully span the dynamically-determined detection window" help="--no-partial-observations; default=use all observations, dividing partial support across matching haplotypes when generating haplotypes. " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
389 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
390 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
391 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
392 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
393 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
394
2dd40b601766 Uploaded
devteam
parents:
diff changeset
395 <!-- indel realignment -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
396
2dd40b601766 Uploaded
devteam
parents:
diff changeset
397 <param name="O" type="boolean" truevalue="-O" falsevalue="" checked="False" label="Turn off left-alignment of indels?" help="-O --dont-left-align-indels; default=False (do left align). " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
398
2dd40b601766 Uploaded
devteam
parents:
diff changeset
399 <!-- input filters -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
400
2dd40b601766 Uploaded
devteam
parents:
diff changeset
401 <conditional name="input_filters">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
402 <param name="input_filters_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Set input filters?" help="Sets -4, -m, -q, -R, -Y, -Q, -U, -z, -&#36;, -e, -0, -F, -C, -3, -G, and -&#33; options " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
403 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
404 <param name="use_duplicate_reads" type="boolean" truevalue="--use-duplicate-reads" falsevalue="" checked="False" label="Include duplicate-marked alignments in the analysis." help="-4 --use-duplicate-reads; default=False (exclude duplicates marked as such in alignments)." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
405 <param name="m" type="integer" value="1" label="Exclude alignments from analysis if they have a mapping quality less than" help="-m --min-mapping-quality; default=1" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
406 <param name="q" type="integer" value="0" label="Exclude alleles from analysis if their supporting base quality less than" help="-q --min-base-quality; default=0" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
407 <param name="R" type="integer" value="0" label="Consider any allele in which the sum of qualities of supporting observations is at least" help="-R --min-supporting-allele-qsum; default=0" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
408 <param name="Y" type="integer" value="0" label="Consider any allele in which and the sum of mapping qualities of supporting reads is at least" help="-Y --min-supporting-mapping-qsum; default=0" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
409 <conditional name="mismatch_filters">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
410 <param name="mismatch_filters_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Perform mismatch filtering?" help="Sets -Q, -U, -z, and &#36; options" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
411 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
412 <param name="Q" type="integer" value="10" label="Count mismatches toward -U (option below) if the base quality of the mismatch is >=" help="-Q --mismatch-base-quality-threshold; default=10" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
413 <param name="U" type="integer" value="1000" optional="True" label="Exclude reads with more than N mismatches where each mismatch has base quality >= Q (option above)" help="-U --read-mismatch-limit; default=~unbound" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
414 <param name="z" type="float" value="1.0" min="0.0" max="1.0" label="Exclude reads with more than N [0,1] fraction of mismatches where each mismatch has base quality >= Q (second option above)" help="-z --read-max-mismatch-fraction; default=1.0" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
415 <param name="read_snp_limit" type="integer" value="1000" label="Exclude reads with more than N base mismatches, ignoring gaps with quality >= Q (third option abobe)" help="-$amp; --read-snp-limit N " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
416 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
417 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
418 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
419 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
420 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
421 <param name="e" type="integer" value="1000" label="Exclude reads with more than this number of separate gaps" help="-e --read-snp-limit; default=~unbounded" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
422 <param name="standard_filters" type="boolean" truevalue="-0" falsevalue="" checked="False" label="Use stringent input base and mapping quality filters" help="-0 --standard-filters; default=False. Equivalent to -m 30 -q 20 -R 0 -S 0" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
423 <param name="F" type="float" value="0.2" label="Require at least this fraction of observations supporting an alternate allele within a single individual in the in order to evaluate the position" help="-F --min-alternate-fraction; default=0.2" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
424 <param name="C" type="integer" value="2" label="Require at least this count of observations supporting an alternate allele within a single individual in order to evaluate the position" help="-C --min-alternate-count; default=2" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
425 <param name="min_alternate_qsum" type="integer" value="0" label="Require at least this sum of quality of observations supporting an alternate allele within a single individual in order to evaluate the position" help="-3 --min-alternate-qsum; default=0" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
426 <param name="G" type="integer" value="1" label="Require at least this count of observations supporting an alternate allele within the total population in order to use the allele in analysis" help="-G --min-alternate-total N; default=1" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
427 <param name="min_coverage" type="integer" value="0" label="Require at least this coverage to process a site" help="-! --min-coverage; default=0 " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
428 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
429 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
430 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
431 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
432 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
433
2dd40b601766 Uploaded
devteam
parents:
diff changeset
434 <!-- population and mappability priors -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
435
2dd40b601766 Uploaded
devteam
parents:
diff changeset
436 <conditional name="population_mappability_priors">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
437 <param name="population_mappability_priors_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Set population and mappability priors?" help="Sets -k, -w, -V, and -a options " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
438 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
439 <param name="k" type="boolean" truevalue="-k" falsevalue="" checked="False" label="No population priors" help="-k --no-population-priors; default=False. Equivalent to --pooled-discrete --hwe-priors-off and removal of Ewens Sampling Formula component of priors." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
440 <param name="w" type="boolean" truevalue="-w" falsevalue="" checked="False" label="Disable estimation of the probability of the combination arising under HWE given the allele frequency as estimated by observation frequency" help="-w --hwe-priors-off; default=False" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
441 <param name="V" type="boolean" truevalue="-V" falsevalue="" checked="False" label="Disable incorporation of prior expectations about observations" help="-V --binomial-obs-priors-off; default=False. Uses read placement probability, strand balance probability, and read position (5&#39;'-3&#39;') probability." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
442 <param name="a" type="boolean" truevalue="-a" falsevalue="" checked="False" label="isable use of aggregate probability of observation balance between alleles as a component of the priors" help="-a --allele-balance-priors-off; default=False " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
443 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
444 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
445 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
446 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
447 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
448
2dd40b601766 Uploaded
devteam
parents:
diff changeset
449 <!-- genotype likelihoods -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
450
2dd40b601766 Uploaded
devteam
parents:
diff changeset
451 <conditional name="genotype_likelihoods">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
452 <param name="genotype_likelihoods_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Tweak genotype likelihoods?" help="Sets --base-quality-cap, --experimental-gls, and --prob-contamination options. " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
453 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
454 <param name="base_quality_cap" type="integer" value="0" label="Limit estimated observation quality by capping base quality at" help="--base-quality-cap" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
455 <param name="experimental_gls" type="boolean" truevalue="--experimental-gls" falsevalue="" checked="False" label="Generate genotype likelihoods using 'effective base depth' metric qual = 1-BaseQual * 1-MapQual" help="--experimental-gls; Incorporate partial observations. This is the default when contamination estimates are provided. Optimized for diploid samples." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
456 <param name="prob_contamination" type="float" value="10e-9" label="An estimate of contamination to use for all samples. " help="--prob-contamination; default=10e-9." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
457 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
458 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
459 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
460 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
461 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
462
2dd40b601766 Uploaded
devteam
parents:
diff changeset
463 <!-- algorithmic features -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
464
2dd40b601766 Uploaded
devteam
parents:
diff changeset
465 <conditional name="algorithmic_features">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
466 <param name="algorithmic_features_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Tweak agrithmic features?" help="Sets --report-genotypes-likelihood-max, -B, --genotyping-max-banddepth, -W, -N, S, -j, -H, -D, -= options " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
467 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
468 <param name="report_genotype_likelihood_max" type="boolean" truevalue="--report-genotype-likelihood-max" falsevalue="" checked="False" label="Report genotypes using the maximum-likelihood estimate provided from genotype likelihoods." help="--report-genotype-likelihood-max; default=False" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
469 <param name="B" type="integer" value="1000" label="Iterate no more than N times during genotyping step" help="-B --genotyping-max-iterations; default=1000." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
470 <param name="genotyping_max_banddepth" type="integer" value="6" label="Integrate no deeper than the Nth best genotype by likelihood when genotyping" help="--genotyping-max-banddepth; default=6" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
471 <param name="W" type="text" size="8" value="1,3" label="Integrate all genotype combinations in our posterior space which include no more than N (1) samples with their Mth (3) best data likelihood" help="-W --posterior-integration-limits; default=1,3" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
472 <param name="N" type="boolean" truevalue="--exclude-unobserved-genotypes" falsevalue="" checked="False" label="Skip sample genotypings for which the sample has no supporting reads" help="-N --exclude-unobserved-genotypes; default=False" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
473 <conditional name="genotype_variant_threshold">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
474 <param name="genotype_variant_threshold_selector" type="boolean" truevalue="set" falsevalue="do_not_set" label="Do you want to to limit posterior integration" help="-S --genotype-variant-threshold" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
475 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
476 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
477 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
478 <when value="set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
479 <param name="S" value="" type="integer" label="Limit posterior integration to samples where the second-best genotype likelihood is no more than log(N) from the highest genotype likelihood for the sample." help="-S --genotype-variant-threshold; default=~unbounded" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
480 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
481 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
482 <param name="j" type="boolean" truevalue="-j" falsevalue="" checked="False" label="Use mapping quality of alleles when calculating data likelihoods" help="-j --use-mapping-quality; default=False" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
483 <param name="H" type="boolean" truevalue="-H" falsevalue="" checked="False" label="Use a weighted sum of base qualities around an indel, scaled by the distance from the indel" help="-H --harmonic-indel-quality; default=use a minimum Base Quality in flanking sequence." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
484 <param name="D" type="float" value="0.9" label="Incorporate non-independence of reads by scaling successive observations by this factor during data likelihood calculations" help="-D --read-dependence-factor; default=0.9." />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
485 <param name="genotype_qualities" type="boolean" truevalue="--genotype-qualities" falsevalue="" checked="False" label="Calculate the marginal probability of genotypes and report as GQ in each sample field in the VCF output" help="-= --genotype-qualities; default=False " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
486 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
487 <when value="do_not_set">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
488 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
489 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
490 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
491 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
492 <when value="simple">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
493 <!-- do nothing -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
494 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
495 <when value="simple_w_filters">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
496 <!-- add standard-filters to command line -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
497 <param name="min_coverage" type="integer" value="0" label="Require at least this coverage to process a site" help="-! --min-coverage; default=0 " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
498 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
499 <when value="naive">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
500 <!-- do nothing build command line using haplotype-length 0 min-alternate-count 1 min-alternate-fraction 0 pooled-continuous report-monomorphic -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
501 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
502 <when value="naive_w_filters">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
503 <!-- do nothing build command line using haplotype-length 0 min-alternate-count 1 min-alternate-fraction 0 pooled-continuous report-monomorphic standard-filters-->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
504 <param name="min_coverage" type="integer" value="0" label="Require at least this coverage to process a site" help="-! --min-coverage; default=0 " />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
505 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
506 <when value="cline">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
507
2dd40b601766 Uploaded
devteam
parents:
diff changeset
508 <expand macro="optional_file_inputs" /> <!-- see macros section -->
2dd40b601766 Uploaded
devteam
parents:
diff changeset
509
2dd40b601766 Uploaded
devteam
parents:
diff changeset
510 <param name="cline" size="60" type="text" value="-m 20 -q 30" label="Type command line tags here" help="All paremeters that DO NOT involve filenames can be typed here. Use &quot;Do you want to provide additional inputs?&quot; section above to control input and output files. For full syntax check help section below">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
511 <sanitizer>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
512 <valid initial="string.printable">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
513 <remove value="&apos;"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
514 </valid>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
515 <mapping initial="none">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
516 <add source="&apos;" target="__sq__"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
517 </mapping>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
518 </sanitizer>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
519 </param>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
520 </when>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
521
2dd40b601766 Uploaded
devteam
parents:
diff changeset
522 </conditional>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
523
2dd40b601766 Uploaded
devteam
parents:
diff changeset
524 </inputs>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
525 <outputs>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
526 <data format="vcf" name="output_vcf" label="${tool.name} on ${on_string} (variants)" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
527 <data format="bed" name="output_failed_alleles_bed" label="${tool.name} on ${on_string} (failed alleles)">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
528 <filter>( options_type['options_type_selector'] == 'cline' or options_type['options_type_selector'] == 'full' ) and options_type['optional_inputs']['optional_inputs_selector'] is True and options_type['optional_inputs']['output_failed_alleles_option'] is True</filter>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
529 </data>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
530 <data format="txt" name="output_trace" label="${tool.name} on ${on_string} (trace)">
2dd40b601766 Uploaded
devteam
parents:
diff changeset
531 <filter>( options_type['options_type_selector'] == 'cline' or options_type['options_type_selector'] == 'full' ) and options_type['optional_inputs']['optional_inputs_selector'] is True and options_type['optional_inputs']['output_trace_option'] is True</filter>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
532 </data>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
533 </outputs>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
534 <tests>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
535 <test>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
536 <param name="reference_source_selector" value="history" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
537 <param name="ref_file" ftype="fasta" value="freebayes-phix174.fasta"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
538 <param name="input_bam" ftype="bam" value="freebayes-phix174.bam"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
539 <param name="options_type_selector" value="simple"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
540 <output name="output_vcf" file="freebayes-phix174-test1.vcf" compare="contains"/>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
541 </test>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
542 </tests>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
543 <stdio>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
544 <exit_code range="1:" />
2dd40b601766 Uploaded
devteam
parents:
diff changeset
545 </stdio>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
546 <help>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
547 **What it does**
2dd40b601766 Uploaded
devteam
parents:
diff changeset
548
2dd40b601766 Uploaded
devteam
parents:
diff changeset
549 FreeBayes is a Bayesian genetic variant detector designed to find small polymorphisms, specifically SNPs (single-nucleotide polymorphisms), indels (insertions and deletions), MNPs (multi-nucleotide polymorphisms), and complex events (composite insertion and substitution events) smaller than the length of a short-read sequencing alignment.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
550
2dd40b601766 Uploaded
devteam
parents:
diff changeset
551 See https://github.com/ekg/freebayes for details on FreeBayes.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
552
2dd40b601766 Uploaded
devteam
parents:
diff changeset
553 This Galaxy instance of FreeBayes corresponds to release 8a407cf5f4416b5eba5bf27ca80144cd5e75bb80
2dd40b601766 Uploaded
devteam
parents:
diff changeset
554
2dd40b601766 Uploaded
devteam
parents:
diff changeset
555 ------
2dd40b601766 Uploaded
devteam
parents:
diff changeset
556
2dd40b601766 Uploaded
devteam
parents:
diff changeset
557 **Description**
2dd40b601766 Uploaded
devteam
parents:
diff changeset
558
2dd40b601766 Uploaded
devteam
parents:
diff changeset
559 Privided BAM file(s) and a reference. FreeBayes will provide VCF output on standard out describing SNPs, indels, and complex variants in samples in the input alignments.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
560
2dd40b601766 Uploaded
devteam
parents:
diff changeset
561 By default, FreeBayes will consider variants supported by at least 2 observations in a single sample (-C) and also by at least 20% of the reads from a single sample (-F). These settings are suitable to low to high depth sequencing in haploid and diploid samples, but users working with polyploid or pooled samples may wish to adjust them depending on the characteristics of their sequencing data.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
562
2dd40b601766 Uploaded
devteam
parents:
diff changeset
563 FreeBayes is capable of calling variant haplotypes shorter than a read length where multiple polymorphisms segregate on the same read. The maximum distance between polymorphisms phased in this way is determined by the --max-complex-gap, which defaults to 3bp. In practice, this can comfortably be set to half the read length.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
564
2dd40b601766 Uploaded
devteam
parents:
diff changeset
565 Ploidy may be set to any level (-p), but by default all samples are assumed to be diploid. FreeBayes can model per-sample and per-region variation in copy-number (-A) using a copy-number variation map.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
566
2dd40b601766 Uploaded
devteam
parents:
diff changeset
567 FreeBayes can act as a frequency-based pooled caller and describe variants and haplotypes in terms of observation frequency rather than called genotypes. To do so, use --pooled-continuous and set input filters to a suitable level. Allele observation counts will be described by AO and RO fields in the VCF output.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
568
2dd40b601766 Uploaded
devteam
parents:
diff changeset
569 -------
2dd40b601766 Uploaded
devteam
parents:
diff changeset
570
2dd40b601766 Uploaded
devteam
parents:
diff changeset
571 **Galaxy-specific options**
2dd40b601766 Uploaded
devteam
parents:
diff changeset
572
2dd40b601766 Uploaded
devteam
parents:
diff changeset
573 Galaxy allows six levels of control over FreeBayes options provided by **Choose parameter selection level** menu option. These are:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
574
2dd40b601766 Uploaded
devteam
parents:
diff changeset
575 1. *Simple diploid calling*: The simples possible FreeBayes application. Equvalent of using FreeBayes with only a BAM input and no other parameter options.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
576 2. *Simple diploid calling with filtering and coverage*: Same as #1 plus two additional options: -0 (standard filters: --min-mapping-quality 30 --min-base-quality 20 --min-supporting-allele-qsum 0 --genotype-varinat-threshold 0) and --min-coverage.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
577 3. *Frequency-based pooled calling*: This is equivalent to using FreeBayes with the following options: --haplotype-length 0 --min-alternate-count 1 --min-alternate-fraction 0 --pooled-continuous --report-monomorphic. This is the best choice for calling varinats in mixtures such as viral, bacterial, or organellar genomes.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
578 4. *Frequency-based pooled calling with filtering and coverage*: Same as #3 but adds -0 and --min-coverage like in #2.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
579 5. *Complete list of all options*: Gives you full control by exposing all FreeBayes options as Galaxy widgets.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
580 6. *Input parameters on the command line*: Similar to the choice above but for those who does not like clicking. Here options can be directly typed into a text box.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
581
2dd40b601766 Uploaded
devteam
parents:
diff changeset
582 -----
2dd40b601766 Uploaded
devteam
parents:
diff changeset
583
2dd40b601766 Uploaded
devteam
parents:
diff changeset
584 **FreeBayes options**
2dd40b601766 Uploaded
devteam
parents:
diff changeset
585
2dd40b601766 Uploaded
devteam
parents:
diff changeset
586 .. class:: infomark
2dd40b601766 Uploaded
devteam
parents:
diff changeset
587
2dd40b601766 Uploaded
devteam
parents:
diff changeset
588 Note that each Galaxy parameter widget corresponding to command line flags listed below:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
589
2dd40b601766 Uploaded
devteam
parents:
diff changeset
590 Input and output::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
591
2dd40b601766 Uploaded
devteam
parents:
diff changeset
592 -t --targets FILE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
593 Limit analysis to targets listed in the BED-format FILE.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
594 -r --region chrom:start_position-end_position
2dd40b601766 Uploaded
devteam
parents:
diff changeset
595 Limit analysis to the specified region, 0-base coordinates,
2dd40b601766 Uploaded
devteam
parents:
diff changeset
596 end_position included. Either '-' or '..' maybe used as a separator.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
597 -s --samples FILE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
598 Limit analysis to samples listed (one per line) in the FILE.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
599 By default FreeBayes will analyze all samples in its input
2dd40b601766 Uploaded
devteam
parents:
diff changeset
600 BAM files.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
601 --populations FILE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
602 Each line of FILE should list a sample and a population which
2dd40b601766 Uploaded
devteam
parents:
diff changeset
603 it is part of. The population-based bayesian inference model
2dd40b601766 Uploaded
devteam
parents:
diff changeset
604 will then be partitioned on the basis of the populations.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
605 -A --cnv-map FILE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
606 Read a copy number map from the BED file FILE, which has
2dd40b601766 Uploaded
devteam
parents:
diff changeset
607 the format:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
608 reference sequence, start, end, sample name, copy number
2dd40b601766 Uploaded
devteam
parents:
diff changeset
609 ... for each region in each sample which does not have the
2dd40b601766 Uploaded
devteam
parents:
diff changeset
610 default copy number as set by --ploidy.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
611 --trace FILE Output an algorithmic trace to FILE.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
612 --failed-alleles FILE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
613 Write a BED file of the analyzed positions which do not
2dd40b601766 Uploaded
devteam
parents:
diff changeset
614 pass --pvar to FILE.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
615 -@ --variant-input VCF
2dd40b601766 Uploaded
devteam
parents:
diff changeset
616 Use variants reported in VCF file as input to the algorithm.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
617 Variants in this file will be treated as putative variants
2dd40b601766 Uploaded
devteam
parents:
diff changeset
618 even if there is not enough support in the data to pass
2dd40b601766 Uploaded
devteam
parents:
diff changeset
619 input filters.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
620 -l --only-use-input-alleles
2dd40b601766 Uploaded
devteam
parents:
diff changeset
621 Only provide variant calls and genotype likelihoods for sites
2dd40b601766 Uploaded
devteam
parents:
diff changeset
622 and alleles which are provided in the VCF input, and provide
2dd40b601766 Uploaded
devteam
parents:
diff changeset
623 output in the VCF for all input alleles, not just those which
2dd40b601766 Uploaded
devteam
parents:
diff changeset
624 have support in the data.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
625 --haplotype-basis-alleles VCF
2dd40b601766 Uploaded
devteam
parents:
diff changeset
626 When specified, only variant alleles provided in this input
2dd40b601766 Uploaded
devteam
parents:
diff changeset
627 VCF will be used for the construction of complex or haplotype
2dd40b601766 Uploaded
devteam
parents:
diff changeset
628 alleles.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
629 --report-all-haplotype-alleles
2dd40b601766 Uploaded
devteam
parents:
diff changeset
630 At sites where genotypes are made over haplotype alleles,
2dd40b601766 Uploaded
devteam
parents:
diff changeset
631 provide information about all alleles in output, not only
2dd40b601766 Uploaded
devteam
parents:
diff changeset
632 those which are called.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
633 --report-monomorphic
2dd40b601766 Uploaded
devteam
parents:
diff changeset
634 Report even loci which appear to be monomorphic, and report all
2dd40b601766 Uploaded
devteam
parents:
diff changeset
635 considered alleles, even those which are not in called genotypes.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
636 Loci which do not have any potential alternates have '.' for ALT.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
637
2dd40b601766 Uploaded
devteam
parents:
diff changeset
638 Reporting::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
639
2dd40b601766 Uploaded
devteam
parents:
diff changeset
640 -P --pvar N Report sites if the probability that there is a polymorphism
2dd40b601766 Uploaded
devteam
parents:
diff changeset
641 at the site is greater than N. default: 0.0. Note that post-
2dd40b601766 Uploaded
devteam
parents:
diff changeset
642 filtering is generally recommended over the use of this parameter.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
643
2dd40b601766 Uploaded
devteam
parents:
diff changeset
644 Population model::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
645
2dd40b601766 Uploaded
devteam
parents:
diff changeset
646 -T --theta N The expected mutation rate or pairwise nucleotide diversity
2dd40b601766 Uploaded
devteam
parents:
diff changeset
647 among the population under analysis. This serves as the
2dd40b601766 Uploaded
devteam
parents:
diff changeset
648 single parameter to the Ewens Sampling Formula prior model
2dd40b601766 Uploaded
devteam
parents:
diff changeset
649 default: 0.001
2dd40b601766 Uploaded
devteam
parents:
diff changeset
650 -p --ploidy N Sets the default ploidy for the analysis to N. default: 2
2dd40b601766 Uploaded
devteam
parents:
diff changeset
651 -J --pooled-discrete
2dd40b601766 Uploaded
devteam
parents:
diff changeset
652 Assume that samples result from pooled sequencing.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
653 Model pooled samples using discrete genotypes across pools.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
654 When using this flag, set --ploidy to the number of
2dd40b601766 Uploaded
devteam
parents:
diff changeset
655 alleles in each sample or use the --cnv-map to define
2dd40b601766 Uploaded
devteam
parents:
diff changeset
656 per-sample ploidy.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
657 -K --pooled-continuous
2dd40b601766 Uploaded
devteam
parents:
diff changeset
658 Output all alleles which pass input filters, regardles of
2dd40b601766 Uploaded
devteam
parents:
diff changeset
659 genotyping outcome or model.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
660
2dd40b601766 Uploaded
devteam
parents:
diff changeset
661 Reference allele::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
662
2dd40b601766 Uploaded
devteam
parents:
diff changeset
663 -Z --use-reference-allele
2dd40b601766 Uploaded
devteam
parents:
diff changeset
664 This flag includes the reference allele in the analysis as
2dd40b601766 Uploaded
devteam
parents:
diff changeset
665 if it is another sample from the same population.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
666 --reference-quality MQ,BQ
2dd40b601766 Uploaded
devteam
parents:
diff changeset
667 Assign mapping quality of MQ to the reference allele at each
2dd40b601766 Uploaded
devteam
parents:
diff changeset
668 site and base quality of BQ. default: 100,60
2dd40b601766 Uploaded
devteam
parents:
diff changeset
669
2dd40b601766 Uploaded
devteam
parents:
diff changeset
670 Allele scope::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
671
2dd40b601766 Uploaded
devteam
parents:
diff changeset
672 -I --no-snps Ignore SNP alleles.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
673 -i --no-indels Ignore insertion and deletion alleles.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
674 -X --no-mnps Ignore multi-nuceotide polymorphisms, MNPs.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
675 -u --no-complex Ignore complex events (composites of other classes).
2dd40b601766 Uploaded
devteam
parents:
diff changeset
676 -n --use-best-n-alleles N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
677 Evaluate only the best N SNP alleles, ranked by sum of
2dd40b601766 Uploaded
devteam
parents:
diff changeset
678 supporting quality scores. (Set to 0 to use all; default: all)
2dd40b601766 Uploaded
devteam
parents:
diff changeset
679 -E --max-complex-gap N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
680 --haplotype-length N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
681 Allow haplotype calls with contiguous embedded matches of up
2dd40b601766 Uploaded
devteam
parents:
diff changeset
682 to this length. (default: 3)
2dd40b601766 Uploaded
devteam
parents:
diff changeset
683 --min-repeat-length N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
684 When assembling observations across repeats, require the total repeat
2dd40b601766 Uploaded
devteam
parents:
diff changeset
685 length at least this many bp. (default: 5)
2dd40b601766 Uploaded
devteam
parents:
diff changeset
686 --min-repeat-entropy N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
687 To detect interrupted repeats, build across sequence until it has
2dd40b601766 Uploaded
devteam
parents:
diff changeset
688 entropy > N bits per bp. (default: 0, off)
2dd40b601766 Uploaded
devteam
parents:
diff changeset
689 --no-partial-observations
2dd40b601766 Uploaded
devteam
parents:
diff changeset
690 Exclude observations which do not fully span the dynamically-determined
2dd40b601766 Uploaded
devteam
parents:
diff changeset
691 detection window. (default, use all observations, dividing partial
2dd40b601766 Uploaded
devteam
parents:
diff changeset
692 support across matching haplotypes when generating haplotypes.)
2dd40b601766 Uploaded
devteam
parents:
diff changeset
693
2dd40b601766 Uploaded
devteam
parents:
diff changeset
694 Indel realignment::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
695
2dd40b601766 Uploaded
devteam
parents:
diff changeset
696 -O --dont-left-align-indels
2dd40b601766 Uploaded
devteam
parents:
diff changeset
697 Turn off left-alignment of indels, which is enabled by default.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
698
2dd40b601766 Uploaded
devteam
parents:
diff changeset
699 Input filters::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
700
2dd40b601766 Uploaded
devteam
parents:
diff changeset
701 -4 --use-duplicate-reads
2dd40b601766 Uploaded
devteam
parents:
diff changeset
702 Include duplicate-marked alignments in the analysis.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
703 default: exclude duplicates marked as such in alignments
2dd40b601766 Uploaded
devteam
parents:
diff changeset
704 -m --min-mapping-quality Q
2dd40b601766 Uploaded
devteam
parents:
diff changeset
705 Exclude alignments from analysis if they have a mapping
2dd40b601766 Uploaded
devteam
parents:
diff changeset
706 quality less than Q. default: 1
2dd40b601766 Uploaded
devteam
parents:
diff changeset
707 -q --min-base-quality Q
2dd40b601766 Uploaded
devteam
parents:
diff changeset
708 Exclude alleles from analysis if their supporting base
2dd40b601766 Uploaded
devteam
parents:
diff changeset
709 quality is less than Q. default: 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
710 -R --min-supporting-allele-qsum Q
2dd40b601766 Uploaded
devteam
parents:
diff changeset
711 Consider any allele in which the sum of qualities of supporting
2dd40b601766 Uploaded
devteam
parents:
diff changeset
712 observations is at least Q. default: 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
713 -Y --min-supporting-mapping-qsum Q
2dd40b601766 Uploaded
devteam
parents:
diff changeset
714 Consider any allele in which and the sum of mapping qualities of
2dd40b601766 Uploaded
devteam
parents:
diff changeset
715 supporting reads is at least Q. default: 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
716 -Q --mismatch-base-quality-threshold Q
2dd40b601766 Uploaded
devteam
parents:
diff changeset
717 Count mismatches toward --read-mismatch-limit if the base
2dd40b601766 Uploaded
devteam
parents:
diff changeset
718 quality of the mismatch is >= Q. default: 10
2dd40b601766 Uploaded
devteam
parents:
diff changeset
719 -U --read-mismatch-limit N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
720 Exclude reads with more than N mismatches where each mismatch
2dd40b601766 Uploaded
devteam
parents:
diff changeset
721 has base quality >= mismatch-base-quality-threshold.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
722 default: ~unbounded
2dd40b601766 Uploaded
devteam
parents:
diff changeset
723 -z --read-max-mismatch-fraction N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
724 Exclude reads with more than N [0,1] fraction of mismatches where
2dd40b601766 Uploaded
devteam
parents:
diff changeset
725 each mismatch has base quality >= mismatch-base-quality-threshold
2dd40b601766 Uploaded
devteam
parents:
diff changeset
726 default: 1.0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
727 -$ --read-snp-limit N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
728 Exclude reads with more than N base mismatches, ignoring gaps
2dd40b601766 Uploaded
devteam
parents:
diff changeset
729 with quality >= mismatch-base-quality-threshold.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
730 default: ~unbounded
2dd40b601766 Uploaded
devteam
parents:
diff changeset
731 -e --read-indel-limit N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
732 Exclude reads with more than N separate gaps.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
733 default: ~unbounded
2dd40b601766 Uploaded
devteam
parents:
diff changeset
734 -0 --standard-filters Use stringent input base and mapping quality filters
2dd40b601766 Uploaded
devteam
parents:
diff changeset
735 Equivalent to -m 30 -q 20 -R 0 -S 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
736 -F --min-alternate-fraction N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
737 Require at least this fraction of observations supporting
2dd40b601766 Uploaded
devteam
parents:
diff changeset
738 an alternate allele within a single individual in the
2dd40b601766 Uploaded
devteam
parents:
diff changeset
739 in order to evaluate the position. default: 0.2
2dd40b601766 Uploaded
devteam
parents:
diff changeset
740 -C --min-alternate-count N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
741 Require at least this count of observations supporting
2dd40b601766 Uploaded
devteam
parents:
diff changeset
742 an alternate allele within a single individual in order
2dd40b601766 Uploaded
devteam
parents:
diff changeset
743 to evaluate the position. default: 2
2dd40b601766 Uploaded
devteam
parents:
diff changeset
744 -3 --min-alternate-qsum N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
745 Require at least this sum of quality of observations supporting
2dd40b601766 Uploaded
devteam
parents:
diff changeset
746 an alternate allele within a single individual in order
2dd40b601766 Uploaded
devteam
parents:
diff changeset
747 to evaluate the position. default: 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
748 -G --min-alternate-total N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
749 Require at least this count of observations supporting
2dd40b601766 Uploaded
devteam
parents:
diff changeset
750 an alternate allele within the total population in order
2dd40b601766 Uploaded
devteam
parents:
diff changeset
751 to use the allele in analysis. default: 1
2dd40b601766 Uploaded
devteam
parents:
diff changeset
752 -! --min-coverage N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
753 Require at least this coverage to process a site. default: 0
2dd40b601766 Uploaded
devteam
parents:
diff changeset
754
2dd40b601766 Uploaded
devteam
parents:
diff changeset
755 Population priors::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
756
2dd40b601766 Uploaded
devteam
parents:
diff changeset
757 -k --no-population-priors
2dd40b601766 Uploaded
devteam
parents:
diff changeset
758 Equivalent to --pooled-discrete --hwe-priors-off and removal of
2dd40b601766 Uploaded
devteam
parents:
diff changeset
759 Ewens Sampling Formula component of priors.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
760
2dd40b601766 Uploaded
devteam
parents:
diff changeset
761 Mappability priors::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
762
2dd40b601766 Uploaded
devteam
parents:
diff changeset
763 -w --hwe-priors-off
2dd40b601766 Uploaded
devteam
parents:
diff changeset
764 Disable estimation of the probability of the combination
2dd40b601766 Uploaded
devteam
parents:
diff changeset
765 arising under HWE given the allele frequency as estimated
2dd40b601766 Uploaded
devteam
parents:
diff changeset
766 by observation frequency.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
767 -V --binomial-obs-priors-off
2dd40b601766 Uploaded
devteam
parents:
diff changeset
768 Disable incorporation of prior expectations about observations.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
769 Uses read placement probability, strand balance probability,
2dd40b601766 Uploaded
devteam
parents:
diff changeset
770 and read position (5'-3') probability.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
771 -a --allele-balance-priors-off
2dd40b601766 Uploaded
devteam
parents:
diff changeset
772 Disable use of aggregate probability of observation balance between alleles
2dd40b601766 Uploaded
devteam
parents:
diff changeset
773 as a component of the priors.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
774
2dd40b601766 Uploaded
devteam
parents:
diff changeset
775 Genotype likelihoods::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
776
2dd40b601766 Uploaded
devteam
parents:
diff changeset
777 --observation-bias FILE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
778 Read length-dependent allele observation biases from FILE.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
779 The format is [length] [alignment efficiency relative to reference]
2dd40b601766 Uploaded
devteam
parents:
diff changeset
780 where the efficiency is 1 if there is no relative observation bias.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
781 --base-quality-cap Q
2dd40b601766 Uploaded
devteam
parents:
diff changeset
782 Limit estimated observation quality by capping base quality at Q.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
783 --experimental-gls
2dd40b601766 Uploaded
devteam
parents:
diff changeset
784 Generate genotype likelihoods using 'effective base depth' metric
2dd40b601766 Uploaded
devteam
parents:
diff changeset
785 qual = 1-BaseQual * 1-MapQual. Incorporate partial observations.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
786 This is the default when contamination estimates are provided.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
787 Optimized for diploid samples.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
788 --prob-contamination F
2dd40b601766 Uploaded
devteam
parents:
diff changeset
789 An estimate of contamination to use for all samples. default: 10e-9
2dd40b601766 Uploaded
devteam
parents:
diff changeset
790 --contamination-estimates FILE
2dd40b601766 Uploaded
devteam
parents:
diff changeset
791 A file containing per-sample estimates of contamination, such as
2dd40b601766 Uploaded
devteam
parents:
diff changeset
792 those generated by VerifyBamID. The format should be:
2dd40b601766 Uploaded
devteam
parents:
diff changeset
793 sample p(read=R|genotype=AR) p(read=A|genotype=AA)
2dd40b601766 Uploaded
devteam
parents:
diff changeset
794 Sample '*' can be used to set default contamination estimates.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
795
2dd40b601766 Uploaded
devteam
parents:
diff changeset
796 Algorithmic features::
2dd40b601766 Uploaded
devteam
parents:
diff changeset
797
2dd40b601766 Uploaded
devteam
parents:
diff changeset
798 --report-genotype-likelihood-max
2dd40b601766 Uploaded
devteam
parents:
diff changeset
799 Report genotypes using the maximum-likelihood estimate provided
2dd40b601766 Uploaded
devteam
parents:
diff changeset
800 from genotype likelihoods.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
801 -B --genotyping-max-iterations N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
802 Iterate no more than N times during genotyping step. default: 1000.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
803 --genotyping-max-banddepth N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
804 Integrate no deeper than the Nth best genotype by likelihood when
2dd40b601766 Uploaded
devteam
parents:
diff changeset
805 genotyping. default: 6.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
806 -W --posterior-integration-limits N,M
2dd40b601766 Uploaded
devteam
parents:
diff changeset
807 Integrate all genotype combinations in our posterior space
2dd40b601766 Uploaded
devteam
parents:
diff changeset
808 which include no more than N samples with their Mth best
2dd40b601766 Uploaded
devteam
parents:
diff changeset
809 data likelihood. default: 1,3.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
810 -N --exclude-unobserved-genotypes
2dd40b601766 Uploaded
devteam
parents:
diff changeset
811 Skip sample genotypings for which the sample has no supporting reads.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
812 -S --genotype-variant-threshold N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
813 Limit posterior integration to samples where the second-best
2dd40b601766 Uploaded
devteam
parents:
diff changeset
814 genotype likelihood is no more than log(N) from the highest
2dd40b601766 Uploaded
devteam
parents:
diff changeset
815 genotype likelihood for the sample. default: ~unbounded
2dd40b601766 Uploaded
devteam
parents:
diff changeset
816 -j --use-mapping-quality
2dd40b601766 Uploaded
devteam
parents:
diff changeset
817 Use mapping quality of alleles when calculating data likelihoods.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
818 -H --harmonic-indel-quality
2dd40b601766 Uploaded
devteam
parents:
diff changeset
819 Use a weighted sum of base qualities around an indel, scaled by the
2dd40b601766 Uploaded
devteam
parents:
diff changeset
820 distance from the indel. By default use a minimum BQ in flanking sequence.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
821 -D --read-dependence-factor N
2dd40b601766 Uploaded
devteam
parents:
diff changeset
822 Incorporate non-independence of reads by scaling successive
2dd40b601766 Uploaded
devteam
parents:
diff changeset
823 observations by this factor during data likelihood
2dd40b601766 Uploaded
devteam
parents:
diff changeset
824 calculations. default: 0.9
2dd40b601766 Uploaded
devteam
parents:
diff changeset
825 -= --genotype-qualities
2dd40b601766 Uploaded
devteam
parents:
diff changeset
826 Calculate the marginal probability of genotypes and report as GQ in
2dd40b601766 Uploaded
devteam
parents:
diff changeset
827 each sample field in the VCF output.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
828
2dd40b601766 Uploaded
devteam
parents:
diff changeset
829
2dd40b601766 Uploaded
devteam
parents:
diff changeset
830 ------
2dd40b601766 Uploaded
devteam
parents:
diff changeset
831
2dd40b601766 Uploaded
devteam
parents:
diff changeset
832 **Citation**
2dd40b601766 Uploaded
devteam
parents:
diff changeset
833
2dd40b601766 Uploaded
devteam
parents:
diff changeset
834 For the underlying tool, please cite `Erik Garrison and Gabor Marth. Haplotype-based variant detection from short-read sequencing &lt;http://arxiv.org/abs/1207.3907&gt;`_.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
835
2dd40b601766 Uploaded
devteam
parents:
diff changeset
836 The initial version of the wrapper was produced by Dan Blankenberg and upgraded by Anton Nekrutenko.
2dd40b601766 Uploaded
devteam
parents:
diff changeset
837
2dd40b601766 Uploaded
devteam
parents:
diff changeset
838 </help>
2dd40b601766 Uploaded
devteam
parents:
diff changeset
839 </tool>