annotate freebayes.xml @ 20:8fb829f953d6 draft

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