# HG changeset patch # User iuc # Date 1550782336 18000 # Node ID 4422d017fe79fda6a67ed86c5833096d895f3153 # Parent 410cdd7cf2bc5609f904ee27ee0dab8243de2b21 planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/bcftools commit 2684e1443f03bfe2ae20c31d23817415ec8f7e69 diff -r 410cdd7cf2bc -r 4422d017fe79 bcftools_mpileup.xml --- a/bcftools_mpileup.xml Tue Dec 12 13:53:56 2017 -0500 +++ b/bcftools_mpileup.xml Thu Feb 21 15:52:16 2019 -0500 @@ -1,5 +1,5 @@ - + Generate VCF or BCF containing genotype likelihoods for one or multiple alignment (BAM or CRAM) files mpileup @@ -16,7 +16,7 @@ - + @@ -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 -]]> - + ]]> - + ]]> diff -r 410cdd7cf2bc -r 4422d017fe79 macros.xml --- a/macros.xml Tue Dec 12 13:53:56 2017 -0500 +++ b/macros.xml Thu Feb 21 15:52:16 2019 -0500 @@ -1,5 +1,5 @@ - 1.4.0 + 1.9 @@ -10,13 +10,13 @@ - bcftools - htslib + bcftools + htslib - samtools + samtools bcftools 2>&1 | grep 'Version:' @@ -39,7 +39,7 @@ ]]> - + @@ -68,7 +71,7 @@ - + $input_vcf && bcftools index $input_vcf && - #elif $input_file.is_of_type('vcf_bgz') - ln -s '$input_file' $input_vcf + #elif $input_file.is_of_type('vcf_bgzip') + ln -s '$input_file' $input_vcf && + #if $input_file.metadata.tabix_index: + ln -s '${input_file.metadata.tabix_index}' ${input_vcf}.tbi && + #else + bcftools index $input_vcf && + #end if #elif $input_file.is_of_type('bcf') #set $input_vcf = 'input' + str($i) + '.bcf.gz' ln -s '$input_file' $input_vcf && @@ -90,8 +98,6 @@ #else bcftools index $input_vcf && #end if - #elif $input_file.is_of_type('bcfvcf_bgz') - ln -s '$input_file' $input_vcf && #end if echo '$input_vcf' >> $vcfs_list_file && $input_vcfs.append($input_vcf) @@ -106,7 +112,7 @@ - + - + @@ -165,7 +171,7 @@ - + #if 'estimate_AF' in $section and $section.estimate_AF: @@ -174,7 +180,7 @@ - + - + #if 'ploidy_file' in $section and $section.ploidy_file: @@ -208,7 +214,7 @@ - + @@ -224,8 +230,8 @@ - + ^([^ \t\n\r\f\v,]+(,[^ \t\n\r\f\v,]+)*)?$ @@ -271,13 +277,13 @@ - ^(\w+(:\d+(-\d+)?)?(,\w+(:\d+(-\d+)?)?)*)?$ - + @@ -346,7 +352,7 @@ - ^(\w+(:\d+(-\d+)?)?(,\w+(:\d+(-\d+)?)?)*)?$ @@ -373,16 +379,16 @@ - + ^(\w+(,\w+)*)?$ - - + help="Inverts the query/filtering applied by Samples (adds "^" prefix to exclude)" /> + + #set $samples_defined = False @@ -397,7 +403,7 @@ - + #if $section.sample: @@ -407,7 +413,7 @@ - + ^[^']*$ @@ -419,7 +425,7 @@ - + ^[^']*$ @@ -431,8 +437,8 @@ - + ^([^,]+(,[^,]+)*)?$ diff -r 410cdd7cf2bc -r 4422d017fe79 test-data/summary.pdf Binary file test-data/summary.pdf has changed diff -r 410cdd7cf2bc -r 4422d017fe79 test-data/view.bcf Binary file test-data/view.bcf has changed diff -r 410cdd7cf2bc -r 4422d017fe79 test-data/view.bcf.csi Binary file test-data/view.bcf.csi has changed diff -r 410cdd7cf2bc -r 4422d017fe79 test-data/view.vcf_bgzip Binary file test-data/view.vcf_bgzip has changed