annotate haplotype_caller.xml @ 5:84584664264c draft

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