annotate haplotype_caller.xml @ 9:fc0e2a67ace8 draft

Uploaded
author iuc
date Tue, 07 Apr 2015 22:46:05 -0400
parents b80301676614
children 0d369d08ad6e
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
6
b80301676614 Uploaded
iuc
parents: 5
diff changeset
1 <tool id="gatk2_haplotype_caller" name="Haplotype Caller" version="@VERSION@.1">
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
2 <description>Call SNPs and indels simultaneously via local de-novo assembly of haplotypes in an active region</description>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
3 <expand macro="requirements" />
6
b80301676614 Uploaded
iuc
parents: 5
diff changeset
4 <expand macro="version_command" />
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
5 <macros>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
6 <import>gatk2_macros.xml</import>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
7 </macros>
1
c8752e31f496 Uploaded
bgruening
parents: 0
diff changeset
8 <command interpreter="python">
c8752e31f496 Uploaded
bgruening
parents: 0
diff changeset
9 gatk2_wrapper.py
c8752e31f496 Uploaded
bgruening
parents: 0
diff changeset
10 --stdout "${output_log}"
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
11 #for $i, $input_bam in enumerate( $reference_source.input_bams ):
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
12 -d "-I" "${input_bam.input_bam}" "${input_bam.input_bam.ext}" "gatk_input_${i}"
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
13 #if str( $input_bam.input_bam.metadata.bam_index ) != "None":
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
14 -d "" "${input_bam.input_bam.metadata.bam_index}" "bam_index" "gatk_input_${i}" ##hardcode galaxy ext type as bam_index
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
15 #end if
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
16 #end for
1
c8752e31f496 Uploaded
bgruening
parents: 0
diff changeset
17 -p '
c8752e31f496 Uploaded
bgruening
parents: 0
diff changeset
18 @JAR_PATH@
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
19 -T "HaplotypeCaller"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
20 -o "${output_vcf}"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
21
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
22 \$GATK2_SITE_OPTIONS
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
23
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
24 --num_cpu_threads_per_data_thread \${GALAXY_SLOTS:-4}
1
c8752e31f496 Uploaded
bgruening
parents: 0
diff changeset
25
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
26 #if $reference_source.reference_source_selector != "history":
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
27 -R "${reference_source.ref_file.fields.path}"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
28 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
29 #if str($input_recal) != 'None':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
30 --BQSR "${input_recal}"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
31 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
32 '
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
33 @DBSNP_OPTIONS@
6
b80301676614 Uploaded
iuc
parents: 5
diff changeset
34 $allow_n_cigar_reads
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
35 #include source=$standard_gatk_options#
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
36
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
37 ##start analysis specific options
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
38 #if $analysis_param_type.analysis_param_type_selector == "advanced":
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
39 -p '
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
40 #if $analysis_param_type.heterozygosity.__str__.strip() != '':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
41 --heterozygosity $analysis_param_type.heterozygosity
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
42 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
43 --genotyping_mode "${analysis_param_type.genotyping_mode_type.genotyping_mode}"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
44 #if str( $analysis_param_type.genotyping_mode_type.genotyping_mode ) == 'GENOTYPE_GIVEN_ALLELES':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
45 --alleles "${analysis_param_type.genotyping_mode_type.input_alleles_rod}"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
46 #end if
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
47 #if not $analysis_param_type.emitRefConfidence is None:
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
48 --emitRefConfidence $analysis_param_type.emitRefConfidence
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
49 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
50
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
51 ## files
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
52 #if str($analysis_param_type.activeRegionIn) != 'None':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
53 --activeRegionIn "$analysis_param_type.activeRegionIn"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
54 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
55 #if str($analysis_param_type.comp) != 'None':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
56 --comp "$analysis_param_type.comp"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
57 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
58 ##
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
59 #if str( $analysis_param_type.annotation ) != "None":
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
60 #for $annotation in str( $analysis_param_type.annotation.fields.gatk_value ).split( ','):
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
61 --annotation "${annotation}"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
62 #end for
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
63 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
64 #for $additional_annotation in $analysis_param_type.additional_annotations:
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
65 --annotation "${additional_annotation.additional_annotation_name}"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
66 #end for
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
67 #if str( $analysis_param_type.group ) != "None":
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
68 #for $group in str( $analysis_param_type.group ).split( ','):
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
69 --group "${group}"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
70 #end for
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
71 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
72 #if str( $analysis_param_type.exclude_annotations ) != "None":
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
73 #for $annotation in str( $analysis_param_type.exclude_annotations.fields.gatk_value ).split( ','):
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
74 --excludeAnnotation "${annotation}"
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
75 #end for
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
76 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
77
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
78 ## value setings
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
79 #if $analysis_param_type.contamination_fraction_to_filter.__str__.strip() != '':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
80 --contamination_fraction_to_filter $analysis_param_type.contamination_fraction_to_filter
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
81 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
82 #if $analysis_param_type.minPruning.__str__.strip() != '':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
83 --minPruning $analysis_param_type.minPruning
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
84 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
85 #if $analysis_param_type.standard_min_confidence_threshold_for_calling.__str__.strip() != '':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
86 --standard_min_confidence_threshold_for_calling $analysis_param_type.standard_min_confidence_threshold_for_calling
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
87 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
88 #if $analysis_param_type.standard_min_confidence_threshold_for_emitting.__str__.strip() != '':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
89 --standard_min_confidence_threshold_for_emitting $analysis_param_type.standard_min_confidence_threshold_for_emitting
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
90 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
91 #if $analysis_param_type.gcpHMM.__str__.strip() != '':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
92 --gcpHMM $analysis_param_type.gcpHMM
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
93 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
94 #if $analysis_param_type.max_alternate_alleles.__str__.strip() != '':
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
95 --max_alternate_alleles $analysis_param_type.max_alternate_alleles
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
96 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
97 ## mode selections
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
98
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
99 #if $analysis_param_type.pair_hmm_implementation.__str__ != "None" and len($analysis_param_type.pair_hmm_implementation.__str__) > 0:
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
100 --pair_hmm_implementation $analysis_param_type.pair_hmm_implementation
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
101 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
102 ## optional outputs
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
103 #if $analysis_param_type.activeRegionOut:
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
104 --activeRegionOut $active_region_out
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
105 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
106 #if $analysis_param_type.graphOutput:
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
107 --graphOutput $graph_out
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
108 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
109 ## flags
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
110 $analysis_param_type.useAllelesTrigger
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
111 $analysis_param_type.fullHaplotype
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
112 $analysis_param_type.genotypeFullActiveRegion
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
113 $analysis_param_type.debug
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
114 '
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
115 #end if
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
116 </command>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
117 <inputs>
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
118 <param name="input_recal" type="data" format="gatk_report" optional="true" label="Covariates table recalibration file" help="The input covariates table file which enables on-the-fly base quality score recalibration. Enables on-the-fly recalibrate of base qualities. The covariates tables are produced by the BaseQualityScoreRecalibrator tool. Please be aware that one should only run recalibration with the covariates file created on the same input bam(s) (-BQSR,--BQSR &amp;lt;recal_file&amp;gt;)" />
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
119 <conditional name="reference_source">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
120 <expand macro="reference_source_selector_param" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
121 <when value="cached">
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
122 <repeat name="input_bams" title="BAM file" min="1" help="-I,--input_file &amp;lt;input_file&amp;gt;">
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
123 <param name="input_bam" type="data" format="bam" label="BAM file">
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
124 <validator type="unspecified_build" />
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
125 <validator type="dataset_metadata_in_data_table" table_name="gatk2_picard_indexes" metadata_name="dbkey" metadata_column="dbkey" message="Sequences are not currently available for the specified build." /> <!-- fixme!!! this needs to be a select -->
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
126 </param>
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
127 </repeat>
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
128 <param name="ref_file" type="select" label="Using reference genome" help="-R,--reference_sequence &amp;lt;reference_sequence&amp;gt;" >
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
129 <options from_data_table="gatk2_picard_indexes">
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
130 <!-- <filter type="data_meta" key="dbkey" ref="input_bam" column="dbkey"/> does not yet work in a repeat...-->
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
131 </options>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
132 <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
133 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
134 </when>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
135 <when value="history">
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
136 <repeat name="input_bams" title="BAM file" min="1" help="-I,--input_file &amp;lt;input_file&amp;gt;">
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
137 <param name="input_bam" type="data" format="bam" label="BAM file" />
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
138 </repeat>
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
139 <param name="ref_file" type="data" format="fasta" label="Using reference file" help="-R,--reference_sequence &amp;lt;reference_sequence&amp;gt;" />
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
140 </when>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
141 </conditional>
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
142 <expand macro="dbsnp_param" />
6
b80301676614 Uploaded
iuc
parents: 5
diff changeset
143
b80301676614 Uploaded
iuc
parents: 5
diff changeset
144 <expand macro="allow_n_cigar_reads" />
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
145 <expand macro="gatk_param_type_conditional" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
146
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
147 <conditional name="analysis_param_type">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
148 <param name="analysis_param_type_selector" type="select" label="Basic or Advanced Analysis options">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
149 <option value="basic" selected="True">Basic</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
150 <option value="advanced">Advanced</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
151 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
152 <when value="basic">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
153 <!-- Do nothing here -->
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
154 </when>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
155 <when value="advanced">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
156
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
157 <param name="activeRegionIn" type="data" format="bed,gatk_interval,picard_interval_list,vcf" optional="true" label="activeRegionIn" help="--activeRegionIn / -AR Use this interval list file as the active regions to process"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
158 <param name="activeRegionOut" type="boolean" checked="False" truevalue="" falsevalue="" label="activeRegionOut" help="--activeRegionOut / -ARO Output the active region to an interval list file"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
159
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
160 <param name="annotation" type="select" multiple="True" display="checkboxes" label="Annotation Types" help="-A,--annotation &amp;lt;annotation&amp;gt;">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
161 <!-- load the available annotations from an external configuration file, since additional ones can be added to local installs -->
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
162 <options from_data_table="gatk2_annotations">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
163 <filter type="multiple_splitter" column="tools_valid_for" separator=","/>
5
c2645201ffae Uploaded
iuc
parents: 4
diff changeset
164 <filter type="static_value" value="HaplotypeCaller" column="tools_valid_for"/>
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
165 </options>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
166 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
167 <repeat name="additional_annotations" title="Additional annotation" help="-A,--annotation &amp;lt;annotation&amp;gt;">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
168 <param name="additional_annotation_name" type="text" value="" label="Annotation name" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
169 </repeat>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
170 <!--
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
171 <conditional name="snpEff_rod_bind_type">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
172 <param name="snpEff_rod_bind_type_selector" type="select" label="Provide a snpEff reference-ordered data file">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
173 <option value="set_snpEff">Set snpEff</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
174 <option value="exclude_snpEff" selected="True">Don't set snpEff</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
175 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
176 <when value="exclude_snpEff">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
177 </when>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
178 <when value="set_snpEff">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
179 <param name="snpEff_input_rod" type="data" format="vcf" label="ROD file" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
180 <param name="snpEff_rod_name" type="hidden" value="snpEff" label="ROD Name"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
181 </when>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
182 </conditional>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
183 -->
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
184 <param name="group" type="select" multiple="True" display="checkboxes" label="Annotation Interfaces/Groups" help="-G,--group &amp;lt;group&amp;gt;">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
185 <option value="RodRequiringAnnotation">RodRequiringAnnotation</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
186 <option value="Standard">Standard</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
187 <option value="Experimental">Experimental</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
188 <option value="WorkInProgress">WorkInProgress</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
189 <option value="RankSumTest">RankSumTest</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
190 <!-- <option value="none">none</option> -->
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
191 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
192 <!-- <param name="family_string" type="text" value="" label="Family String"/> -->
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
193 <param name="exclude_annotations" type="select" multiple="True" display="checkboxes" label="Annotations to exclude" help="-XA,--excludeAnnotation &amp;lt;excludeAnnotation&amp;gt;" >
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
194 <!-- load the available annotations from an external configuration file, since additional ones can be added to local installs -->
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
195 <options from_data_table="gatk2_annotations">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
196 <filter type="multiple_splitter" column="tools_valid_for" separator=","/>
5
c2645201ffae Uploaded
iuc
parents: 4
diff changeset
197 <filter type="static_value" value="HaplotypeCaller" column="tools_valid_for"/>
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
198 </options>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
199 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
200
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
201 <param name="comp" type="data" format="vcf" optional="true" label="comp" help="--comp / -comp comparison VCF file"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
202 <param name="contamination_fraction_to_filter" type="float" value="0.05" optional="true" label="contamination_fraction_to_filter" help="--contamination_fraction_to_filter / -contamination Fraction of contamination in sequencing data (for all samples) to aggressively remove">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
203 <validator type="in_range" message="value between 0.00 and 1.00" min="0" max="1"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
204 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
205 <param name="debug" type="boolean" checked="False" truevalue="-debug" falsevalue="" label="debug" help="--debug / -debug If specified, print out very verbose debug information about each triggering active region"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
206
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
207 <conditional name="genotyping_mode_type">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
208 <param name="genotyping_mode" type="select" label="How to determine the alternate allele to use for genotyping" help="-gt_mode,--genotyping_mode &amp;lt;genotyping_mode&amp;gt;">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
209 <option value="DISCOVERY" selected="True">DISCOVERY</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
210 <option value="GENOTYPE_GIVEN_ALLELES">GENOTYPE_GIVEN_ALLELES</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
211 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
212 <when value="DISCOVERY">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
213 <!-- Do nothing here -->
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
214 </when>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
215 <when value="GENOTYPE_GIVEN_ALLELES">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
216 <param name="input_alleles_rod" type="data" format="vcf" label="Alleles ROD file" help="-alleles,--alleles &amp;lt;alleles&amp;gt;" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
217 </when>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
218 </conditional>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
219
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
220
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
221 <param name="graphOutput" type="boolean" checked="False" truevalue="" falsevalue="" label="graphOutput" help="--graphOutput / -graph File to which debug assembly graph information should be written"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
222 <param name="heterozygosity" type="float" value="0.0010" optional="true" label="heterozygosity" help="--heterozygosity / -hets Heterozygosity value used to compute prior likelihoods for any locus"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
223 <param name="minPruning" type="integer" value="1" optional="true" label="minPruning" help="--minPruning / -minPruning The minimum allowed pruning factor in assembly graph. Paths with &gt;= X supporting kmers are pruned from the graph">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
224 <validator type="in_range" message="value between 0 and 127" min="0" max="127"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
225 </param>
4
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
226 <!-- http://www.broadinstitute.org/gatk/guide/article?id=2940 -->
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
227 <param name="emitRefConfidence" type="select" optional="true" label="Output confidence estimates" help="Emitting a per-bp or summarized confidence estimate for a site being strictly homozygous-reference (--emitRefConfidence)">
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
228 <option value="NONE" selected="True">don't emit anything</option>
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
229 <option value="BP_RESOLUTION">BP_RESOLUTION (emit detailed information for each BP)</option>
e67da4f2c9bf Uploaded
bgruening
parents: 1
diff changeset
230 <option value="GVCF">GVCF (emit a block summarized version of the BP_RESOLUTION data)</option>
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
231 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
232 <param name="pair_hmm_implementation" type="select" optional="true" label="pair_hmm_implementation" help="--pair_hmm_implementation / -pairHMM The PairHMM implementation to use for genotype likelihood calculations">
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
233 <option value="EXACT">EXACT</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
234 <option value="ORIGINAL">ORIGINAL</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
235 <option value="CACHING">CACHING</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
236 <option value="LOGLESS_CACHING" selected="True">LOGLESS_CACHING</option>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
237 </param>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
238 <param name="standard_min_confidence_threshold_for_calling" type="float" value="30.0" optional="true" label="standard_min_confidence_threshold_for_calling" help="--standard_min_confidence_threshold_for_calling / -stand_call_conf The minimum phred-scaled confidence threshold at which variants should be called"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
239 <param name="standard_min_confidence_threshold_for_emitting" type="float" value="30.0" optional="true" label="standard_min_confidence_threshold_for_emitting" help="--standard_min_confidence_threshold_for_emitting / -stand_emit_conf The minimum phred-scaled confidence threshold at which variants should be emitted (and filtered with LowQual if less than the calling threshold)"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
240 <param name="useAllelesTrigger" type="boolean" checked="False" truevalue="-allelesTrigger" falsevalue="" label="useAllelesTrigger" help="--useAllelesTrigger / -allelesTrigger If specified, use additional trigger on variants found in an external alleles file"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
241 <param name="fullHaplotype" type="boolean" checked="False" truevalue="-fullHaplotype" falsevalue="" label="fullHaplotype" help="--fullHaplotype / -fullHaplotype If specified, output the full haplotype sequence instead of converting to individual variants w.r.t. the reference"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
242 <param name="gcpHMM" type="integer" value="10" optional="true" label="gcpHMM" help="--gcpHMM / -gcpHMM Flat gap continuation penalty for use in the Pair HMM"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
243 <param name="genotypeFullActiveRegion" type="boolean" checked="False" truevalue="-genotypeFullActiveRegion" falsevalue="" label="genotypeFullActiveRegion" help="--genotypeFullActiveRegion / -genotypeFullActiveRegion If specified, alternate alleles are considered to be the full active region for the purposes of genotyping"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
244 <param name="max_alternate_alleles" type="integer" value="6" optional="true" label="max_alternate_alleles" help="--max_alternate_alleles / -maxAltAlleles Maximum number of alternate alleles to genotype"/>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
245 </when>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
246 </conditional>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
247 </inputs>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
248 <outputs>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
249 <data format="vcf" name="output_vcf" label="${tool.name} on ${on_string} (VCF)" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
250 <data format="vcf" name="graph_out" label="${tool.name} on ${on_string} graph" >
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
251 <filter>analysis_param_type['analysis_param_type_selector'] == "advanced" and analysis_param_type['graphOutput'] == True</filter>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
252 </data>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
253 <data format="vcf" name="active_region_out" label="${tool.name} on ${on_string} activeRegion" >
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
254 <filter>analysis_param_type['analysis_param_type_selector'] == "advanced" and analysis_param_type['activeRegionOut'] == True</filter>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
255 </data>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
256 <data format="txt" name="output_log" label="${tool.name} on ${on_string} (log)" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
257 </outputs>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
258 <tests>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
259 <test>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
260 <param name="input_recal" value="gatk/gatk_count_covariates/gatk_count_covariates_out_1.csv" ftype="csv" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
261 <param name="reference_source_selector" value="history" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
262 <param name="ref_file" value="phiX.fasta" ftype="fasta" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
263 <param name="input_bam" value="gatk/gatk_indel_realigner/gatk_indel_realigner_out_1.bam" ftype="bam" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
264 <param name="gatk_param_type_selector" value="basic" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
265 <param name="analysis_param_type_selector" value="basic" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
266 <output name="output_bam" file="gatk/gatk_table_recalibration/gatk_table_recalibration_out_1.bam" ftype="bam" lines_diff="4" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
267 <output name="output_log" file="gatk/gatk_table_recalibration/gatk_table_recalibration_out_1.log.contains" compare="contains" />
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
268 </test>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
269 </tests>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
270 <help>
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
271 **What it does**
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
272
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
273 **HaplotypeCaller**
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
274 calls SNPs and indels simultaneously via local de-novo assembly of haplotypes in an active region.
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
275 Haplotypes are evaluated using an affine gap penalty Pair HMM.
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
276
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
277 For more information on using read based compression in the GATK, see this `tool specific page &lt;http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_walkers_haplotypecaller_HaplotypeCaller.html&gt;`_.
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
278
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
279 To learn about best practices for variant detection using GATK, see this `overview &lt;http://www.broadinstitute.org/gatk/guide/topic?name=best-practices&gt;`_.
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
280
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
281 If you encounter errors, please view the `GATK FAQ &lt;http://www.broadinstitute.org/gatk/guide/topic?name=faqs&gt;`_.
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
282
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
283 ------
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
284
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
285 **Inputs**
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
286
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
287 GenomeAnalysisTK: PrintReads accepts aligned BAM files.
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
288
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
289
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
290 **Outputs**
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
291
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
292 The output is a VCF file with raw, unrecalibrated SNP and indel calls.
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
293
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
294
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
295 Go `here &lt;http://www.broadinstitute.org/gatk/guide/topic?name=intro&gt;`_ for details on GATK file formats.
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
296
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
297 -------
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
298
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
299 **Settings**::
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
300
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
301 activeRegionIn Use this interval list file as the active regions to process
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
302 activeRegionOut Output the active region to this interval list file
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
303 alleles The set of alleles at which to genotype when --genotyping_mode is GENOTYPE_GIVEN_ALLELES
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
304 annotation One or more specific annotations to apply to variant calls
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
305 comp comparison VCF file
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
306 contamination Fraction of contamination in sequencing data (for all samples) to aggressively remove
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
307 dbsnp dbSNP file
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
308 debug If specified, print out very verbose debug information about each triggering active region
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
309 excludeAnnotation One or more specific annotations to exclude
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
310 genotyping_mode Specifies how to determine the alternate alleles to use for genotyping
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
311 graphOutput File to which debug assembly graph information should be written
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
312 group One or more classes/groups of annotations to apply to variant calls
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
313 heterozygosity Heterozygosity value used to compute prior likelihoods for any locus
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
314 minPruning The minimum allowed pruning factor in assembly graph. Paths with less than or equal supporting kmers are pruned from the graph
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
315 pair_hmm_implementation The PairHMM implementation to use for genotype likelihood calculations
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
316 stand_call_conf The minimum phred-scaled confidence threshold at which variants should be called
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
317 stand_emit_conf The minimum phred-scaled confidence threshold at which variants should be emitted (and filtered with LowQual if less than the calling threshold)
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
318 useAllelesTrigger If specified, use additional trigger on variants found in an external alleles file
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
319 fullHaplotype If specified, output the full haplotype sequence instead of converting to individual variants w.r.t. the reference
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
320 gcpHMM Flat gap continuation penalty for use in the Pair HMM
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
321 genotypeFullActiveRegion If specified, alternate alleles are considered to be the full active region for the purposes of genotyping
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
322 max_alternate_alleles Maximum number of alternate alleles to genotype
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
323
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
324 @CITATION_SECTION@
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
325 </help>
6
b80301676614 Uploaded
iuc
parents: 5
diff changeset
326 <expand macro="citations" />
0
9ad15075b9fa Uploaded
iuc
parents:
diff changeset
327 </tool>