annotate samtools_mpileup.xml @ 3:973fea5b4bdf draft

Uploaded tool and dependency definitions that specify samtools version 0.1.19.
author devteam
date Thu, 27 Mar 2014 15:27:36 -0400
parents dddde948c3e5
children c6fdfe3331d6
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
3
973fea5b4bdf Uploaded tool and dependency definitions that specify samtools version 0.1.19.
devteam
parents: 1
diff changeset
1 <tool id="samtools_mpileup" name="MPileup" version="0.0.3">
0
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
2 <description>SNP and indel caller</description>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
3 <requirements>
3
973fea5b4bdf Uploaded tool and dependency definitions that specify samtools version 0.1.19.
devteam
parents: 1
diff changeset
4 <requirement type="package" version="0.1.19">samtools</requirement>
0
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
5 </requirements>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
6 <command interpreter="python">samtools_wrapper.py
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
7 -p 'samtools mpileup'
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
8 --stdout "${output_log}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
9 #if $reference_source.reference_source_selector != "history":
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
10 -p '-f "${reference_source.ref_file.fields.path}"'
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
11 #else:
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
12 -d "-f" "${reference_source.ref_file}" "fa" "reference_input"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
13 #end if
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
14 #for $i, $input_bam in enumerate( $reference_source.input_bams ):
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
15 -d " " "${input_bam.input_bam}" "${input_bam.input_bam.ext}" "bam_input_${i}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
16 -d "" "${input_bam.input_bam.metadata.bam_index}" "bam_index" "bam_input_${i}" ##hardcode galaxy ext type as bam_index
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
17 #end for
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
18 -p '
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
19 #if str( $advanced_options.advanced_options_selector ) == "advanced":
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
20 ${advanced_options.skip_anomalous_read_pairs}
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
21 ${advanced_options.disable_probabilistic_realignment}
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
22 -C "${advanced_options.coefficient_for_downgrading}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
23 -d "${advanced_options.max_reads_per_bam}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
24 ${advanced_options.extended_BAQ_computation}
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
25 #if str( $advanced_options.position_list ) != 'None':
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
26 -l "${advanced_options.position_list}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
27 #end if
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
28 -q "${advanced_options.minimum_mapping_quality}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
29 -Q "${advanced_options.minimum_base_quality}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
30 #if str( $advanced_options.region_string ):
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
31 -r "${advanced_options.region_string}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
32 #end if
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
33 ${advanced_options.output_per_sample_read_depth}
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
34 ${advanced_options.output_per_sample_strand_bias_p_value}
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
35 #end if
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
36 #if str( $genotype_likelihood_computation_type.genotype_likelihood_computation_type_selector ) == 'perform_genotype_likelihood_computation':
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
37 ##-g or -u
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
38 -g
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
39 -e "${genotype_likelihood_computation_type.gap_extension_sequencing_error_probability}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
40 -h "${genotype_likelihood_computation_type.coefficient_for_modeling_homopolymer_errors}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
41 #if str( $genotype_likelihood_computation_type.perform_indel_calling.perform_indel_calling_selector ) == 'perform_indel_calling':
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
42 -L "${genotype_likelihood_computation_type.perform_indel_calling.skip_indel_calling_above_sample_depth}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
43 #else:
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
44 -I
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
45 #end if
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
46 -o "${genotype_likelihood_computation_type.gap_open_sequencing_error_probability}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
47 #if len( $genotype_likelihood_computation_type.platform_list_repeat ):
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
48 -P "${ ",".join( [ str( platform.platform_entry ) for platform in $genotype_likelihood_computation_type.platform_list_repeat ] ) }"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
49 #end if
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
50 #end if
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
51 &gt; "${output_mpileup}"
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
52 '
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
53 </command>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
54 <inputs>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
55 <conditional name="reference_source">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
56 <param name="reference_source_selector" type="select" label="Choose the source for the reference list">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
57 <option value="cached">Locally cached</option>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
58 <option value="history">History</option>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
59 </param>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
60 <when value="cached">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
61 <repeat name="input_bams" title="BAM file" min="1">
1
dddde948c3e5 Update samtools_mpileup to use the fasta_indexes data table.
Dave Bouvier <dave@bx.psu.edu>
parents: 0
diff changeset
62 <param name="input_bam" type="data" format="bam" label="BAM file">
dddde948c3e5 Update samtools_mpileup to use the fasta_indexes data table.
Dave Bouvier <dave@bx.psu.edu>
parents: 0
diff changeset
63 <validator type="unspecified_build" />
dddde948c3e5 Update samtools_mpileup to use the fasta_indexes data table.
Dave Bouvier <dave@bx.psu.edu>
parents: 0
diff changeset
64 <validator type="dataset_metadata_in_data_table" table_name="fasta_indexes" metadata_name="dbkey" metadata_column="1" message="Sequences are not currently available for the specified build." /> <!-- fixme!!! this needs to be a select -->
dddde948c3e5 Update samtools_mpileup to use the fasta_indexes data table.
Dave Bouvier <dave@bx.psu.edu>
parents: 0
diff changeset
65 </param>
0
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
66 </repeat>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
67 <param name="ref_file" type="select" label="Using reference genome">
1
dddde948c3e5 Update samtools_mpileup to use the fasta_indexes data table.
Dave Bouvier <dave@bx.psu.edu>
parents: 0
diff changeset
68 <options from_data_table="fasta_indexes">
dddde948c3e5 Update samtools_mpileup to use the fasta_indexes data table.
Dave Bouvier <dave@bx.psu.edu>
parents: 0
diff changeset
69 <!-- <filter type="data_meta" ref="input_bam" key="dbkey" column="1" /> does not yet work in a repeat...-->
0
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
70 </options>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
71 </param>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
72 </when>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
73 <when value="history"> <!-- FIX ME!!!! -->
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
74 <repeat name="input_bams" title="BAM file" min="1">
1
dddde948c3e5 Update samtools_mpileup to use the fasta_indexes data table.
Dave Bouvier <dave@bx.psu.edu>
parents: 0
diff changeset
75 <param name="input_bam" type="data" format="bam" label="BAM file">
dddde948c3e5 Update samtools_mpileup to use the fasta_indexes data table.
Dave Bouvier <dave@bx.psu.edu>
parents: 0
diff changeset
76 <validator type="metadata" check="bam_index" message="Metadata missing, click the pencil icon in the history item and use the auto-detect feature to correct this issue." />
dddde948c3e5 Update samtools_mpileup to use the fasta_indexes data table.
Dave Bouvier <dave@bx.psu.edu>
parents: 0
diff changeset
77 </param>
0
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
78 </repeat>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
79 <param name="ref_file" type="data" format="fasta" label="Using reference file" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
80 </when>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
81 </conditional>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
82
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
83
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
84 <conditional name="genotype_likelihood_computation_type">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
85 <param name="genotype_likelihood_computation_type_selector" type="select" label="Genotype Likelihood Computation">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
86 <option value="perform_genotype_likelihood_computation">Perform genotype likelihood computation</option>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
87 <option value="do_not_perform_genotype_likelihood_computation" selected="True">Do not perform genotype likelihood computation</option>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
88 </param>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
89 <when value="perform_genotype_likelihood_computation">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
90 <param name="gap_extension_sequencing_error_probability" type="integer" value="20" label="Phred-scaled gap extension sequencing error probability" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
91 <param name="coefficient_for_modeling_homopolymer_errors" type="integer" value="100" label="Coefficient for modeling homopolymer errors." />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
92 <conditional name="perform_indel_calling">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
93 <param name="perform_indel_calling_selector" type="select" label="Perform INDEL calling">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
94 <option value="perform_indel_calling" selected="True">Perform INDEL calling</option>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
95 <option value="do_not_perform_indel_calling">Do not perform INDEL calling</option>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
96 </param>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
97 <when value="perform_indel_calling">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
98 <param name="skip_indel_calling_above_sample_depth" type="integer" value="250" label="Skip INDEL calling if the average per-sample depth is above" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
99 </when>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
100 <when value="do_not_perform_indel_calling" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
101 </conditional>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
102 <param name="gap_open_sequencing_error_probability" type="integer" value="40" label="Phred-scaled gap open sequencing error probability" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
103 <repeat name="platform_list_repeat" title="Platform for INDEL candidates">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
104 <param name="platform_entry" type="text" value="" label="Platform to use for INDEL candidates" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
105 </repeat>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
106 </when>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
107 <when value="do_not_perform_genotype_likelihood_computation">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
108 <!-- Do nothing here -->
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
109 </when>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
110 </conditional>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
111 <conditional name="advanced_options">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
112 <param name="advanced_options_selector" type="select" label="Set advanced options">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
113 <option value="basic" selected="True">Basic</option>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
114 <option value="advanced">Advanced</option>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
115 </param>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
116 <when value="advanced">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
117 <param name="skip_anomalous_read_pairs" type="boolean" truevalue="-A" falsevalue="" checked="False" label="Do not skip anomalous read pairs in variant calling" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
118 <param name="disable_probabilistic_realignment" type="boolean" truevalue="-B" falsevalue="" checked="False" label=" Disable probabilistic realignment for the computation of base alignment quality (BAQ)" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
119 <param name="coefficient_for_downgrading" type="integer" value="0" label="Coefficient for downgrading mapping quality for reads containing excessive mismatches" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
120 <param name="max_reads_per_bam" type="integer" value="250" label="Max reads per BAM" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
121 <param name="extended_BAQ_computation" type="boolean" truevalue="-E" falsevalue="" checked="False" label="Extended BAQ computation" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
122 <param name="position_list" type="data" format="bed" label="List of regions or sites on which to operate" optional="True" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
123 <param name="minimum_mapping_quality" type="integer" value="0" label="Minimum mapping quality for an alignment to be used" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
124 <param name="minimum_base_quality" type="integer" value="13" label="Minimum base quality for a base to be considered" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
125 <param name="region_string" type="text" value="" label="Only generate pileup in region" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
126 <param name="output_per_sample_read_depth" type="boolean" truevalue="-D" falsevalue="" checked="False" label="Output per-sample read depth" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
127 <param name="output_per_sample_strand_bias_p_value" type="boolean" truevalue="-S" falsevalue="" checked="False" label="Output per-sample Phred-scaled strand bias P-value" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
128 </when>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
129 <when value="basic" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
130 </conditional>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
131 </inputs>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
132 <outputs>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
133 <data format="pileup" name="output_mpileup" label="${tool.name} on ${on_string}">
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
134 <change_format>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
135 <when input="genotype_likelihood_computation_type.genotype_likelihood_computation_type_selector" value="perform_genotype_likelihood_computation" format="bcf" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
136 </change_format>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
137 </data>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
138 <data format="txt" name="output_log" label="${tool.name} on ${on_string} (log)" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
139 </outputs>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
140 <tests>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
141 <test>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
142 <param name="reference_source_selector" value="history" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
143 <param name="ref_file" value="phiX.fasta" ftype="fasta" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
144 <param name="input_bam" value="gatk/gatk_table_recalibration/gatk_table_recalibration_out_1.bam" ftype="bam" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
145 <param name="genotype_likelihood_computation_type_selector" value="do_not_perform_genotype_likelihood_computation" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
146 <param name="advanced_options_selector" value="basic" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
147 <output name="output_mpileup" file="samtools/mpileup/samtools_mpileup_out_1.pileup" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
148 <output name="output_log" file="samtools/mpileup/samtools_mpileup_out_1.log" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
149 </test>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
150 <test>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
151 <param name="reference_source_selector" value="history" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
152 <param name="ref_file" value="phiX.fasta" ftype="fasta" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
153 <param name="input_bam" value="gatk/gatk_table_recalibration/gatk_table_recalibration_out_1.bam" ftype="bam" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
154 <param name="genotype_likelihood_computation_type_selector" value="perform_genotype_likelihood_computation" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
155 <param name="gap_extension_sequencing_error_probability" value="20" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
156 <param name="coefficient_for_modeling_homopolymer_errors" value="100" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
157 <param name="perform_indel_calling_selector" value="perform_indel_calling" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
158 <param name="skip_indel_calling_above_sample_depth" value="250" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
159 <param name="gap_open_sequencing_error_probability" value="40" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
160 <param name="platform_list_repeat" value="0" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
161 <param name="advanced_options_selector" value="basic" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
162 <output name="output_mpileup" file="samtools/mpileup/samtools_mpileup_out_2.bcf" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
163 <output name="output_log" file="samtools/mpileup/samtools_mpileup_out_1.log" />
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
164 </test>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
165 </tests>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
166 <help>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
167 **What it does**
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
168
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
169 Generate BCF or pileup for one or multiple BAM files. Alignment records are grouped by sample identifiers in @RG header lines. If sample identifiers are absent, each input file is regarded as one sample.
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
170
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
171 ------
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
172
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
173 **Settings**::
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
174
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
175 Input Options:
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
176 -6 Assume the quality is in the Illumina 1.3+ encoding.
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
177 -A Do not skip anomalous read pairs in variant calling.
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
178 -B Disable probabilistic realignment for the computation of base alignment quality (BAQ). BAQ is the Phred-scaled probability of a read base being misaligned. Applying this option greatly helps to reduce false SNPs caused by misalignments.
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
179 -b FILE List of input BAM files, one file per line [null]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
180 -C INT Coefficient for downgrading mapping quality for reads containing excessive mismatches. Given a read with a phred-scaled probability q of being generated from the mapped position, the new mapping quality is about sqrt((INT-q)/INT)*INT. A zero value disables this functionality; if enabled, the recommended value for BWA is 50. [0]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
181 -d INT At a position, read maximally INT reads per input BAM. [250]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
182 -E Extended BAQ computation. This option helps sensitivity especially for MNPs, but may hurt specificity a little bit.
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
183 -f FILE The faidx-indexed reference file in the FASTA format. The file can be optionally compressed by razip. [null]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
184 -l FILE BED or position list file containing a list of regions or sites where pileup or BCF should be generated [null]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
185 -q INT Minimum mapping quality for an alignment to be used [0]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
186 -Q INT Minimum base quality for a base to be considered [13]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
187 -r STR Only generate pileup in region STR [all sites]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
188 Output Options:
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
189
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
190 -D Output per-sample read depth
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
191 -g Compute genotype likelihoods and output them in the binary call format (BCF).
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
192 -S Output per-sample Phred-scaled strand bias P-value
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
193 -u Similar to -g except that the output is uncompressed BCF, which is preferred for piping.
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
194
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
195 Options for Genotype Likelihood Computation (for -g or -u):
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
196
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
197 -e INT Phred-scaled gap extension sequencing error probability. Reducing INT leads to longer indels. [20]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
198 -h INT Coefficient for modeling homopolymer errors. Given an l-long homopolymer run, the sequencing error of an indel of size s is modeled as INT*s/l. [100]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
199 -I Do not perform INDEL calling
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
200 -L INT Skip INDEL calling if the average per-sample depth is above INT. [250]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
201 -o INT Phred-scaled gap open sequencing error probability. Reducing INT leads to more indel calls. [40]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
202 -P STR Comma dilimited list of platforms (determined by @RG-PL) from which indel candidates are obtained. It is recommended to collect indel candidates from sequencing technologies that have low indel error rate such as ILLUMINA. [all]
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
203
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
204 ------
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
205
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
206 **Citation**
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
207
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
208 For the underlying tool, please cite `Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R; 1000 Genome Project Data Processing Subgroup. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009 Aug 15;25(16):2078-9. &lt;http://www.ncbi.nlm.nih.gov/pubmed/19505943&gt;`_
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
209
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
210 If you use this tool in Galaxy, please cite Blankenberg D, et al. *In preparation.*
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
211
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
212 </help>
44a18a94d7a9 Uploaded tool tarball.
devteam
parents:
diff changeset
213 </tool>