annotate print_reads.xml @ 15:01ff8dd37d4d draft default tip

Uploaded
author lz_hust
date Sat, 01 Jun 2019 07:20:41 -0400
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
15
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
1 <tool id="gatk2_print_reads" name="Print Reads" version="@VERSION@.0">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
2 <description>on BAM files</description>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
3 <macros>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
4 <import>gatk2_macros.xml</import>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
5 </macros>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
6 <expand macro="requirements" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
7 <expand macro="version_command" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
8 <command interpreter="python">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
9 gatk2_wrapper.py
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
10 --stdout "${output_log}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
11 -d "-I" "${reference_source.input_bam}" "${reference_source.input_bam.ext}" "gatk_input"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
12 #if str( $reference_source.input_bam.metadata.bam_index ) != "None":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
13 -d "" "${reference_source.input_bam.metadata.bam_index}" "bam_index" "gatk_input" ##hardcode galaxy ext type as bam_index
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
14 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
15 -p '
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
16 @JAR_PATH@
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
17 -T "PrintReads"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
18 -o "${output_bam}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
19 \$GATK2_SITE_OPTIONS
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
20
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
21 ## according to http://www.broadinstitute.org/gatk/guide/article?id=1975
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
22 --num_cpu_threads_per_data_thread \${GALAXY_SLOTS:-6}
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
23
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
24 #if $reference_source.reference_source_selector != "history":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
25 -R "${reference_source.ref_file.fields.path}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
26 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
27 #if str($input_recal) != 'None':
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
28 --BQSR "${input_recal}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
29 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
30 --disable_bam_indexing
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
31 '
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
32
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
33 #include source=$standard_gatk_options#
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
34
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
35 ##start analysis specific options
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
36 #if $analysis_param_type.analysis_param_type_selector == "advanced":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
37 -p '
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
38 #if $analysis_param_type.default_read_group_type.default_read_group_type_selector == "set":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
39 --default_read_group "${analysis_param_type.default_read_group_type.default_read_group}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
40 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
41 #if str( $analysis_param_type.default_platform ) != "default":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
42 --default_platform "${analysis_param_type.default_platform}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
43 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
44 #if str( $analysis_param_type.force_read_group_type.force_read_group_type_selector ) == "set":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
45 --force_read_group "${analysis_param_type.force_read_group_type.force_read_group}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
46 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
47 #if str( $analysis_param_type.force_platform ) != "default":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
48 --force_platform "${analysis_param_type.force_platform}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
49 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
50 ${analysis_param_type.exception_if_no_tile}
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
51 #if str( $analysis_param_type.solid_options_type.solid_options_type_selector ) == "set":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
52 #if str( $analysis_param_type.solid_options_type.solid_recal_mode ) != "default":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
53 --solid_recal_mode "${analysis_param_type.solid_options_type.solid_recal_mode}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
54 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
55 #if str( $analysis_param_type.solid_options_type.solid_nocall_strategy ) != "default":
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
56 --solid_nocall_strategy "${analysis_param_type.solid_options_type.solid_nocall_strategy}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
57 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
58 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
59 ${analysis_param_type.simplify_bam}
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
60 --preserve_qscores_less_than "${analysis_param_type.preserve_qscores_less_than}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
61 --smoothing "${analysis_param_type.smoothing}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
62 --max_quality_score "${analysis_param_type.max_quality_score}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
63 --window_size_nqs "${analysis_param_type.window_size_nqs}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
64 --homopolymer_nback "${analysis_param_type.homopolymer_nback}"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
65 ${analysis_param_type.do_not_write_original_quals}
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
66 '
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
67 #end if
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
68 </command>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
69 <inputs>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
70 <param name="input_recal" type="data" format="gatk_report" optional="true" label="Covariates table recalibration file"
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
71 help="The input covariates table file which enables on-the-fly base quality score recalibration (intended for use with BaseRecalibrator files) (-BQSR,--BQSR)" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
72 <conditional name="reference_source">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
73 <expand macro="reference_source_selector_param" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
74 <when value="cached">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
75 <param name="input_bam" type="data" format="bam" label="BAM file" help="-I,--input_file &amp;lt;input_file&amp;gt;">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
76 <validator type="unspecified_build" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
77 <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 -->
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
78 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
79 <param name="ref_file" type="select" label="Using reference genome" help="-R,--reference_sequence &amp;lt;reference_sequence&amp;gt;" >
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
80 <options from_data_table="gatk2_picard_indexes">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
81 <filter type="data_meta" key="dbkey" ref="input_bam" column="dbkey"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
82 </options>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
83 <validator type="no_options" message="A built-in reference genome is not available for the build associated with the selected input file"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
84 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
85 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
86 <when value="history">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
87 <param name="input_bam" type="data" format="bam" label="BAM file" help="-I,--input_file &amp;lt;input_file&amp;gt;" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
88 <param name="ref_file" type="data" format="fasta" label="Using reference file" help="-R,--reference_sequence &amp;lt;reference_sequence&amp;gt;">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
89 <options>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
90 <filter type="data_meta" key="dbkey" ref="input_bam" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
91 </options>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
92 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
93 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
94 </conditional>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
95
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
96 <expand macro="gatk_param_type_conditional" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
97
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
98 <conditional name="analysis_param_type">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
99 <param name="analysis_param_type_selector" type="select" label="Basic or Advanced Analysis options">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
100 <option value="basic" selected="True">Basic</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
101 <option value="advanced">Advanced</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
102 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
103 <when value="basic">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
104 <!-- Do nothing here -->
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
105 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
106 <when value="advanced">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
107 <conditional name="default_read_group_type">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
108 <param name="default_read_group_type_selector" type="select" label="Set default Read Group" help="--default_read_group">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
109 <option value="default" selected="True">Don't Set</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
110 <option value="set">Set</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
111 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
112 <when value="default">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
113 <!-- do nothing here -->
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
114 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
115 <when value="set">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
116 <param name="default_read_group" type="text" value="Unknown" label="If a read has no read group then default to the provided String"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
117 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
118 </conditional>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
119 <param name="default_platform" type="select" label="Set default Platform" help="--default_platform">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
120 <option value="default" selected="True">Don't Set</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
121 <option value="illumina">illumina</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
122 <option value="454">454</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
123 <option value="solid">solid</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
124 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
125 <conditional name="force_read_group_type">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
126 <param name="force_read_group_type_selector" type="select" label="Force Read Group" help="--force_read_group">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
127 <option value="default" selected="True">Don't Force</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
128 <option value="set">Force</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
129 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
130 <when value="default">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
131 <!-- do nothing here -->
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
132 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
133 <when value="set">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
134 <param name="force_read_group" type="text" value="Unknown" label="If provided, the read group ID of EVERY read will be forced to be the provided String."/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
135 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
136 </conditional>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
137 <param name="force_platform" type="select" label="Force Platform" help="--force_platform">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
138 <option value="default" selected="True">Don't Force</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
139 <option value="illumina">illumina</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
140 <option value="454">454</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
141 <option value="solid">solid</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
142 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
143 <param name="exception_if_no_tile" type="boolean" checked="False" truevalue="--exception_if_no_tile" falsevalue="" label="Throw an exception when no tile can be found" help="--exception_if_no_tile"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
144 <conditional name="solid_options_type">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
145 <param name="solid_options_type_selector" type="select" label="Set SOLiD specific options">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
146 <option value="default" selected="True">Don't Set</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
147 <option value="set">Set</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
148 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
149 <when value="default">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
150 <!-- do nothing here -->
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
151 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
152 <when value="set">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
153 <param name="solid_recal_mode" type="select" label="How should we recalibrate solid bases in which the reference was inserted" help="-sMode,--solid_recal_mode &amp;lt;solid_recal_mode&amp;gt;">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
154 <option value="default" selected="True">Don't set</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
155 <option value="DO_NOTHING">DO_NOTHING</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
156 <option value="SET_Q_ZERO">SET_Q_ZERO</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
157 <option value="SET_Q_ZERO_BASE_N">SET_Q_ZERO_BASE_N</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
158 <option value="REMOVE_REF_BIAS">REMOVE_REF_BIAS</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
159 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
160 <param name="solid_nocall_strategy" type="select" label="Behavior of the recalibrator when it encounters no calls" help="-solid_nocall_strategy,--solid_nocall_strategy &amp;lt;solid_nocall_strategy&amp;gt;">
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
161 <option value="default" selected="True">Don't set</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
162 <option value="THROW_EXCEPTION">THROW_EXCEPTION</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
163 <option value="LEAVE_READ_UNRECALIBRATED">LEAVE_READ_UNRECALIBRATED</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
164 <option value="PURGE_READ">PURGE_READ</option>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
165 </param>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
166 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
167 </conditional>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
168 <param name="simplify_bam" type="boolean" checked="False" truevalue="-simplifyBAM" falsevalue="" label="Simplify BAM" help="-simplifyBAM,--simplifyBAM"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
169 <param name="window_size_nqs" type="integer" value="5" label="Window size used by MinimumNQSCovariate" help="--window_size_nqs"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
170 <param name="homopolymer_nback" type="integer" value="7" label="Number of previous bases to look at in HomopolymerCovariate" help="-nback,--homopolymer_nback &amp;lt;homopolymer_nback&amp;gt;" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
171 <param name="preserve_qscores_less_than" type="integer" value="5" label="Bases with quality scores less than this threshold won't be recalibrated" help="-pQ,--preserve_qscores_less_than &amp;lt;preserve_qscores_less_than&amp;gt;"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
172 <param name="smoothing" type="integer" value="1" label="smoothing" help="-sm,--smoothing &amp;lt;smoothing&amp;gt;"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
173 <param name="max_quality_score" type="integer" value="50" label="Max quality score" help="-maxQ,--max_quality_score &amp;lt;max_quality_score&amp;gt;"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
174 <param name="do_not_write_original_quals" type="boolean" checked="False" truevalue="--doNotWriteOriginalQuals" falsevalue="" label="Do Not Write Original Quality tag" help="-noOQs,--doNotWriteOriginalQuals"/>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
175 </when>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
176 </conditional>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
177 </inputs>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
178 <outputs>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
179 <data format="bam" name="output_bam" label="${tool.name} on ${on_string} (BAM)" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
180 <data format="txt" name="output_log" label="${tool.name} on ${on_string} (log)" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
181 </outputs>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
182 <tests>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
183 <test>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
184 <param name="input_recal" value="gatk/gatk_count_covariates/gatk_count_covariates_out_1.csv" ftype="csv" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
185 <param name="reference_source_selector" value="history" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
186 <param name="ref_file" value="phiX.fasta" ftype="fasta" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
187 <param name="input_bam" value="gatk/gatk_indel_realigner/gatk_indel_realigner_out_1.bam" ftype="bam" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
188 <param name="gatk_param_type_selector" value="basic" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
189 <param name="analysis_param_type_selector" value="basic" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
190 <output name="output_bam" file="gatk/gatk_table_recalibration/gatk_table_recalibration_out_1.bam" ftype="bam" lines_diff="4" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
191 <output name="output_log" file="gatk/gatk_table_recalibration/gatk_table_recalibration_out_1.log.contains" compare="contains" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
192 </test>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
193 </tests>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
194 <help>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
195 **What it does**
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
196
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
197 This walker is designed to work as the second pass in a two-pass processing step, doing a by-read traversal. For each base in each read this walker calculates various user-specified covariates (such as read group, reported quality score, cycle, and dinuc) Using these values as a key in a large hashmap the walker calculates an empirical base quality score and overwrites the quality score currently in the read. This walker then outputs a new bam file with these updated (recalibrated) reads. Note: This walker expects as input the recalibration table file generated previously by CovariateCounterWalker. Note: This walker is designed to be used in conjunction with CovariateCounterWalker.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
198
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
199 For more information on base quality score recalibration using the GATK, see this `tool specific page &lt;http://www.broadinstitute.org/gatk/gatkdocs/org_broadinstitute_sting_gatk_walkers_readutils_PrintReads.html&gt;`_.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
200
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
201 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;`_.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
202
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
203 If you encounter errors, please view the `GATK FAQ &lt;http://www.broadinstitute.org/gatk/guide/topic?name=faqs&gt;`_.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
204
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
205 ------
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
206
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
207 **Inputs**
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
208
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
209 GenomeAnalysisTK: PrintReads accepts an aligned BAM and a recalibration (gatk_report) input files.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
210
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
211
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
212 **Outputs**
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
213
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
214 The output is in BAM format.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
215
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
216
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
217 Go `here &lt;http://www.broadinstitute.org/gatk/guide/topic?name=intro&gt;`_ for details on GATK file formats.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
218
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
219 -------
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
220
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
221 **Settings**::
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
222
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
223 default_read_group If a read has no read group then default to the provided String.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
224 default_platform If a read has no platform then default to the provided String. Valid options are illumina, 454, and solid.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
225 force_read_group If provided, the read group ID of EVERY read will be forced to be the provided String. This is useful to collapse all data into a single read group.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
226 force_platform If provided, the platform of EVERY read will be forced to be the provided String. Valid options are illumina, 454, and solid.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
227 window_size_nqs The window size used by MinimumNQSCovariate for its calculation
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
228 homopolymer_nback The number of previous bases to look at in HomopolymerCovariate
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
229 exception_if_no_tile If provided, TileCovariate will throw an exception when no tile can be found. The default behavior is to use tile = -1
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
230 solid_recal_mode How should we recalibrate solid bases in whichthe reference was inserted? Options = DO_NOTHING, SET_Q_ZERO, SET_Q_ZERO_BASE_N, or REMOVE_REF_BIAS (DO_NOTHING|SET_Q_ZERO|SET_Q_ZERO_BASE_N|REMOVE_REF_BIAS)
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
231 solid_nocall_strategy Defines the behavior of the recalibrator when it encounters no calls in the color space. Options = THROW_EXCEPTION, LEAVE_READ_UNRECALIBRATED, or PURGE_READ (THROW_EXCEPTION|LEAVE_READ_UNRECALIBRATED|PURGE_READ)
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
232 recal_file Filename for the input covariates table recalibration .gatk_report file
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
233 out The output BAM file
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
234 bam_compression Compression level to use for writing BAM files
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
235 disable_bam_indexing Turn off on-the-fly creation of indices for output BAM files.
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
236 simplifyBAM If provided, output BAM files will be simplified to include just key reads for downstream variation discovery analyses (removing duplicates, PF-, non-primary reads), as well stripping all extended tags from the kept reads except the read group identifier
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
237 preserve_qscores_less_than Bases with quality scores less than this threshold won't be recalibrated, default=5. In general it's unsafe to change qualities scores below &lt; 5, since base callers use these values to indicate random or bad bases
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
238 smoothing Number of imaginary counts to add to each bin bin order to smooth out bins with few data points, default=1
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
239 max_quality_score The integer value at which to cap the quality scores, default=50
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
240 doNotWriteOriginalQuals If true, we will not write the original quality (OQ) tag for each read
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
241
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
242 @CITATION_SECTION@
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
243 </help>
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
244 <expand macro="citations" />
01ff8dd37d4d Uploaded
lz_hust
parents:
diff changeset
245 </tool>