Repository 'bcftools_mpileup'
hg clone https://toolshed.g2.bx.psu.edu/repos/iuc/bcftools_mpileup

Changeset 2:cf06b44624c7 (2019-02-21)
Previous changeset 1:9c711df3258d (2017-12-12) Next changeset 3:4e1d23e5c691 (2019-05-27)
Commit message:
planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bcftools commit 2684e1443f03bfe2ae20c31d23817415ec8f7e69
modified:
bcftools_mpileup.xml
macros.xml
test-data/view.bcf
added:
test-data/summary.pdf
test-data/view.vcf_bgzip
removed:
test-data/view.bcf.csi
b
diff -r 9c711df3258d -r cf06b44624c7 bcftools_mpileup.xml
--- a/bcftools_mpileup.xml Tue Dec 12 14:05:31 2017 -0500
+++ b/bcftools_mpileup.xml Thu Feb 21 16:09:57 2019 -0500
[
@@ -1,5 +1,5 @@
 <?xml version='1.0' encoding='utf-8'?>
-<tool name="bcftools @EXECUTABLE@" id="bcftools_@EXECUTABLE@" version="@VERSION@.0">
+<tool name="bcftools @EXECUTABLE@" id="bcftools_@EXECUTABLE@" version="@TOOL_VERSION@">
     <description>Generate VCF or BCF containing genotype likelihoods for one or multiple alignment (BAM or CRAM) files</description>
     <macros>
         <token name="@EXECUTABLE@">mpileup</token>
@@ -16,7 +16,7 @@
             <option value="256">The alignment or this read is not primary</option>
             <option value="512">The read fails platform/vendor quality checks</option>
             <option value="1024">The read is a PCR or optical duplicate</option>
-         </xml>
+        </xml>
     </macros>
     <expand macro="requirements">
         <expand macro="samtools_requirement"/>
@@ -37,9 +37,9 @@
     #silent $bam_list.append($input_name)
     ln -s '${input.input_bam}' ${input_name} &&
     #if $input.input_bam.ext == 'bam':
-      ln -s '${input.input_bam.metadata.bam_index}' ${input_name}.${idx_ext} &&
+        ln -s '${input.input_bam.metadata.bam_index}' ${input_name}.${idx_ext} &&
     #else:
-      ln -s '${input.input_bam.metadata.cram_index}' ${input_name}.${idx_ext} &&
+        ln -s '${input.input_bam.metadata.cram_index}' ${input_name}.${idx_ext} &&
     #end if
 #else:
     #for $bam_count, $input_bam in enumerate( $input.input_bams ):
@@ -59,11 +59,11 @@
 
 #set $input_fa_ref = None
 #if $reference_source.reference_source_selector == "history":
-  #set $input_fa_ref = 'ref.fa'
-  ln -s '${reference_source.ref_file}' $input_fa_ref &&
-  samtools faidx $input_fa_ref &&
+    #set $input_fa_ref = 'ref.fa'
+    ln -s '${reference_source.ref_file}' $input_fa_ref &&
+    samtools faidx $input_fa_ref &&
 #elif $reference_source.reference_source_selector == "cached":
-  #set $input_fa_ref = $reference_source.ref_file.fields.path
+    #set $input_fa_ref = $reference_source.ref_file.fields.path
 #end if
 
 #set $section = $sec_restrict
@@ -81,53 +81,53 @@
 ## Indel Calling section
 #set $section = $sec_indel
 #if $section.perform_indel_calling.perform_indel_calling_selector == 'do_not_perform_indel_calling':
-  --skip-indels
+    --skip-indels
 #elif $section.perform_indel_calling.perform_indel_calling_selector == 'perform_indel_calling':
-  -o "${section.perform_indel_calling.gap_open_sequencing_error_probability}"
-  -e "${section.perform_indel_calling.gap_extension_sequencing_error_probability}"
-  -h "${section.perform_indel_calling.coefficient_for_modeling_homopolymer_errors}"
-  -L "${section.perform_indel_calling.skip_indel_calling_above_sample_depth}"
-  -m "${section.perform_indel_calling.minimum_gapped_reads_for_indel_candidates}"
-  --open-prob "${section.perform_indel_calling.open_seq_error_probability}"
-  -F "${section.perform_indel_calling.minimum_gapped_read_fraction}"
-  ${section.perform_indel_calling.gapped_read_per_sample}
-  #if len( $section.perform_indel_calling.platform_list_repeat ):
-    -P "${ ",".join( [ str( platform.platform_entry ) for platform in $section.perform_indel_calling.platform_list_repeat ] ) }"
-  #end if
+    -o "${section.perform_indel_calling.gap_open_sequencing_error_probability}"
+    -e "${section.perform_indel_calling.gap_extension_sequencing_error_probability}"
+    -h "${section.perform_indel_calling.coefficient_for_modeling_homopolymer_errors}"
+    -L "${section.perform_indel_calling.skip_indel_calling_above_sample_depth}"
+    -m "${section.perform_indel_calling.minimum_gapped_reads_for_indel_candidates}"
+    --open-prob "${section.perform_indel_calling.open_seq_error_probability}"
+    -F "${section.perform_indel_calling.minimum_gapped_read_fraction}"
+    ${section.perform_indel_calling.gapped_read_per_sample}
+    #if len( $section.perform_indel_calling.platform_list_repeat ):
+        -P "${ ",".join( [ str( platform.platform_entry ) for platform in $section.perform_indel_calling.platform_list_repeat ] ) }"
+    #end if
 #end if
 
 ## Filter section
 #set $section = $sec_filtering
 #if str( $section.filter_by_flags.filter_flags ) == "filter":
-  #if $section.filter_by_flags.require_flags:
-    --rf ${sum([int(flag) for flag in str($section.filter_by_flags.require_flags).split(',')])}
-  #end if
-  #if $section.filter_by_flags.exclude_flags:
-    --ff ${sum([int(flag) for flag in str($section.filter_by_flags.exclude_flags).split(',')])}
-  #end if
+    #if $section.filter_by_flags.require_flags:
+        --rf ${sum([int(flag) for flag in str($section.filter_by_flags.require_flags).split(',')])}
+    #end if
+    #if $section.filter_by_flags.exclude_flags:
+        --ff ${sum([int(flag) for flag in str($section.filter_by_flags.exclude_flags).split(',')])}
+    #end if
 #end if
 -d "${section.max_reads_per_bam}"
 ${section.skip_anomalous_read_pairs}
 #if str( $section.quality.quality_settings ) == "adjust":
-  $section.quality.baq
-  -q "${section.quality.minimum_mapping_quality}"
-  -Q "${section.quality.minimum_base_quality}"
-  -C "${section.quality.coefficient_for_downgrading}"
+    $section.quality.baq
+    -q "${section.quality.minimum_mapping_quality}"
+    -Q "${section.quality.minimum_base_quality}"
+    -C "${section.quality.coefficient_for_downgrading}"
 #end if
 #if str( $section.read_groups.read_groups_selector ) == "ignore_rg":
-  --ignore-RG
+    --ignore-RG
 #elif str( $section.read_groups.read_groups_selector ) == "paste":
-  -G "${section.read_groups.rg_action}${read_groups_file}"
+    -G "${section.read_groups.rg_action}${read_groups_file}"
 #elif str( $section.read_groups.read_groups_selector ) == "history"
-  -G "${section.read_groups.rg_action}${section.read_groups.read_groups}"
+    -G "${section.read_groups.rg_action}${section.read_groups.read_groups}"
 #end if
 
 #set $section = $sec_output_options
 #if $section.output_tags:
- --annotate "$section.output_tags"
+    --annotate "$section.output_tags"
 #end if
 #if $section.gvcf:
- --gvcf "$section.gvcf"
+    --gvcf "$section.gvcf"
 #end if
 
 ## Subset section
@@ -147,18 +147,17 @@
 #echo ' '.join($bam_list)#
 > '$output_file'
 #if str( $sec_filtering.read_groups.read_groups_selector ) == "paste":
-&& echo 'read-groups:'
-&& cat ${read_groups_file}
+    && echo 'read-groups:'
+    && cat ${read_groups_file}
 #end if
-]]>
-    </command>
+    ]]></command>
     <configfiles>
         <configfile name="read_groups_file">
 <![CDATA[#slurp
 #set pasted_data = ''
 #set $section = $sec_filtering
 #if str( $section.read_groups.read_groups_selector ) == "paste":
-  #set pasted_data = '\t'.join( str( $section.read_groups.group_paste).split() )
+    #set pasted_data = '\t'.join( str( $section.read_groups.group_paste).split() )
 #end if
 #slurp
 ${pasted_data}
@@ -444,7 +443,6 @@
 @BCFTOOLS_MANPAGE@#@EXECUTABLE@
 
 @BCFTOOLS_WIKI@
-]]>
-    </help>
+    ]]></help>
     <expand macro="citations" />
 </tool>
b
diff -r 9c711df3258d -r cf06b44624c7 macros.xml
--- a/macros.xml Tue Dec 12 14:05:31 2017 -0500
+++ b/macros.xml Thu Feb 21 16:09:57 2019 -0500
[
b'@@ -1,5 +1,5 @@\n <macros>\n-  <token name="@VERSION@">1.4.0</token>\n+  <token name="@TOOL_VERSION@">1.9</token>\n   <xml name="stdio">\n     <stdio>\n       <exit_code range="1:" />\n@@ -10,13 +10,13 @@\n   </xml>\n   <xml name="requirements">\n     <requirements>\n-      <requirement type="package" version="1.4">bcftools</requirement>\n-      <requirement type="package" version="1.4">htslib</requirement>\n+      <requirement type="package" version="@TOOL_VERSION@">bcftools</requirement>\n+      <requirement type="package" version="1.9">htslib</requirement>\n       <yield />\n     </requirements>\n   </xml>\n   <xml name="samtools_requirement">\n-      <requirement type="package" version="1.3.1">samtools</requirement>\n+      <requirement type="package" version="1.9">samtools</requirement>\n   </xml>\n   <xml name="version_command">\n     <version_command>bcftools 2&gt;&amp;1 | grep \'Version:\'</version_command>\n@@ -39,7 +39,7 @@\n ]]>\n   </token>\n   <xml name="macro_input">\n-    <param name="input_file" type="data" format="vcf,vcf_bgzip,bcf,bcf_bgzip" label="VCF/BCF Data" />\n+    <param name="input_file" type="data" format="vcf,vcf_bgzip,bcf" label="VCF/BCF Data" />\n   </xml>\n   <token name="@PREPARE_INPUT_FILE@">\n <![CDATA[\n@@ -50,6 +50,11 @@\n   bcftools index $input_vcf &&\n #elif $input_file.is_of_type(\'vcf_bgzip\')\n   ln -s \'$input_file\' $input_vcf &&\n+  #if $input_file.metadata.tabix_index:\n+    ln -s \'${input_file.metadata.tabix_index}\' ${input_vcf}.tbi &&\n+  #else\n+    bcftools index $input_vcf &&\n+  #end if\n #elif $input_file.is_of_type(\'bcf\')\n   #set $input_vcf = \'input.bcf\'\n   ln -s \'$input_file\' $input_vcf &&\n@@ -58,8 +63,6 @@\n   #else\n     bcftools index $input_vcf &&\n   #end if\n-#elif $input_file.is_of_type(\'bcf_bgzip\')\n-  ln -s \'$input_file\' $input_vcf &&\n #end if\n ]]>\n   </token>\n@@ -68,7 +71,7 @@\n   </token>\n \n   <xml name="macro_inputs">\n-    <param name="input_files" type="data" format="vcf,bcf" label="Other VCF/BCF Datasets" multiple="True" />\n+    <param name="input_files" type="data" format="vcf,vcf_bgzip,bcf" label="Other VCF/BCF Datasets" multiple="True" />\n   </xml>\n   <token name="@PREPARE_INPUT_FILES@">\n <![CDATA[\n@@ -80,8 +83,13 @@\n   #if $input_file.is_of_type(\'vcf\')\n     bgzip -c \'$input_file\' > $input_vcf &&\n     bcftools index $input_vcf &&\n-  #elif $input_file.is_of_type(\'vcf_bgz\')\n-    ln -s \'$input_file\' $input_vcf\n+  #elif $input_file.is_of_type(\'vcf_bgzip\')\n+    ln -s \'$input_file\' $input_vcf &&\n+    #if $input_file.metadata.tabix_index:\n+      ln -s \'${input_file.metadata.tabix_index}\' ${input_vcf}.tbi &&\n+    #else\n+      bcftools index $input_vcf &&\n+    #end if\n   #elif $input_file.is_of_type(\'bcf\')\n     #set $input_vcf = \'input\' + str($i) + \'.bcf.gz\'\n     ln -s \'$input_file\' $input_vcf &&\n@@ -90,8 +98,6 @@\n     #else\n       bcftools index $input_vcf &&\n     #end if\n-  #elif $input_file.is_of_type(\'bcfvcf_bgz\')\n-    ln -s \'$input_file\' $input_vcf &&\n   #end if\n   echo \'$input_vcf\' >> $vcfs_list_file &&\n   $input_vcfs.append($input_vcf)\n@@ -106,7 +112,7 @@\n   </token>\n \n   <xml name="macro_fasta_ref">\n-    <param name="fasta_ref" argument="--fasta-ref" type="data" format="data" label="Reference sequence in FASTA format" optional="True" />\n+    <param name="fasta_ref" argument="--fasta-ref" type="data" format="data" optional="true" label="Reference sequence in FASTA format" />\n   </xml>\n   <token name="@PREPARE_FASTA_REF@">\n <![CDATA[\n@@ -148,7 +154,7 @@\n \n \n   <xml name="macro_AF_file">\n-    <param name="AF_file" argument="--AF-file" type="data" format="tabular" label="Allele frequencies file" optional="True" help="Tab-delimited file containing the columns CHR,POS,REF,ALT,AF" />\n+    <param name="AF_file" argument="--AF-file" type="data" format="tabular" optional="true" label="Allele frequencies file" help="Tab-delimited file containing the columns CHR,POS,REF,ALT,AF" />\n   </xml>\n   <!-- This may need to bgzip and tabix the file -->\n   <token name="@PREPARE_AF_FILE@">\n@@ -165,7 +171,7 @@\n   </token>\n \n   <xml '..b'ne__"/>\n         <when value="targets">\n-            <param name="targets" type="text" value="" label="Restrict to comma-separated list of targets" optional="true"\n+            <param name="targets" type="text" value="" optional="true" label="Restrict to comma-separated list of targets"\n                    help="Each target is specifed as: chr or chr:pos or chr:from-to">\n                  <validator type="regex" message="">^(\\w+(:\\d+(-\\d+)?)?(,\\w+(:\\d+(-\\d+)?)?)*)?$</validator>\n             </param>\n@@ -373,16 +379,16 @@\n   </token>\n \n   <xml name="macro_samples">\n-      <param name="samples" type="text" value="" label="Samples" optional="true"\n-             help="(-s) comma separated list of samples to annotate (or exclude)">\n+      <param argument="--samples" type="text" value="" optional="true" label="Samples"\n+             help="Comma separated list of samples to annotate (or exclude)">\n           <validator type="regex" message="">^(\\w+(,\\w+)*)?$</validator>\n       </param>\n       <param name="invert_samples" type="boolean" truevalue="^" falsevalue="" checked="false" label="Invert Samples"\n-             help="inverts the query/filtering applied by Samples (adds &quot;^&quot; prefix to exclude)" />\n-      <param name="samples_file" type="data" format="tabular" label="Samples File" optional="True"\n-             help="(-S) file of samples to include" />\n-      <param name="invert_samples_file" type="boolean" truevalue="^" falsevalue="" checked="false" label="Invert Samples File"\n-             help="inverts the query/filtering applied by Samples File" />\n+             help="Inverts the query/filtering applied by Samples (adds &quot;^&quot; prefix to exclude)" />\n+      <param argument="--samples_file" type="data" format="tabular" optional="true" label="Samples file"\n+             help="File of samples to include" />\n+      <param name="invert_samples_file" type="boolean" truevalue="^" falsevalue="" checked="false" label="Invert Samples file"\n+             help="inverts the query/filtering applied by Samples file" />\n   </xml>\n   <token name="@SAMPLES@">\n #set $samples_defined = False\n@@ -397,7 +403,7 @@\n   </token>\n \n   <xml name="macro_sample">\n-      <param name="sample" type="text" label="Sample" optional="True" help="apply variants of the given sample" />\n+      <param name="sample" type="text" optional="true" label="Sample" help="Apply variants of the given sample" />\n   </xml>\n   <token name="@SAMPLE@">\n #if $section.sample:\n@@ -407,7 +413,7 @@\n \n \n   <xml name="macro_include">\n-    <param name="include" type="text" label="Include" optional="True" help="(-i) select sites for which the expression is true">\n+    <param argument="--include" type="text" optional="true" label="Include" help="Select sites for which the expression is true">\n         <validator type="regex" message="Single quote not allowed">^[^\']*$</validator>\n         <sanitizer sanitize="False"/>\n     </param>\n@@ -419,7 +425,7 @@\n   </token>\n \n   <xml name="macro_exclude">\n-    <param name="exclude" type="text" label="Exclude" optional="True" help="(-e) exclude sites for which the expression is true">\n+    <param argument="--exclude" type="text" optional="true" label="Exclude" help="Exclude sites for which the expression is true">\n         <validator type="regex" message="Single quote not allowed">^[^\']*$</validator>\n         <sanitizer sanitize="False"/>\n     </param>\n@@ -431,8 +437,8 @@\n   </token>\n \n   <xml name="macro_columns">\n-    <param name="columns" type="text" value="" label="Columns" optional="true"\n-            help="list of columns in the annotation file, e.g. CHROM,POS,REF,ALT,-,INFO/TAG. See man page for details">\n+    <param name="columns" type="text" value="" optional="true" label="Columns"\n+            help="List of columns in the annotation file, e.g. CHROM,POS,REF,ALT,-,INFO/TAG. See man page for details">\n         <validator type="regex" message="COLUMN names  separated by commas">^([^,]+(,[^,]+)*)?$</validator>\n     </param>\n   </xml>\n'
b
diff -r 9c711df3258d -r cf06b44624c7 test-data/summary.pdf
b
Binary file test-data/summary.pdf has changed
b
diff -r 9c711df3258d -r cf06b44624c7 test-data/view.bcf
b
Binary file test-data/view.bcf has changed
b
diff -r 9c711df3258d -r cf06b44624c7 test-data/view.bcf.csi
b
Binary file test-data/view.bcf.csi has changed
b
diff -r 9c711df3258d -r cf06b44624c7 test-data/view.vcf_bgzip
b
Binary file test-data/view.vcf_bgzip has changed