# HG changeset patch # User wrtz # Date 1452811843 18000 # Node ID 3ebaa7bd2773278a87534cd218eaea95b014b14d # Parent 4c0acb64c5fc56ec04cc6fc9f7accb3c6fbbaf8b Uploaded diff -r 4c0acb64c5fc -r 3ebaa7bd2773 bcftools_merge.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/bcftools_merge.xml Thu Jan 14 17:50:43 2016 -0500 @@ -0,0 +1,75 @@ + + Merge multiple VCF/BCF files + + #for $i, $inp in enumerate( $inputs ): + #if str( $inp.input_file_index ) != "None": + ln -s -f "${inp.input}" "input_file_${i}.${inp.input.ext}" && + ln -s -f "${inp.input_file_index}" "input_file_${i}.${inp.input.ext}.tbi" && + #end if + #end for + bcftools merge + --output "${output_file}" + --output-type "${output_type}" + #if str( $regions_file ) != "None": + --regions-file "${regions_file}" + #end if + #if str( $regions ) != "": + --regions "${regions}" + #end if + ${force_samples} + #if str( $info_rules ) != "": + --info-rules "${info_rules}" + #end if + #if str( $merge ) != "None": + --merge "${merge}" + #end if + + #for $i, $inp in enumerate( $inputs ): + #if str( $inp.input_file_index ) != "None": + "input_file_${i}.${inp.input.ext}" + #else + "${inp.input}" + #end if + #end for + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +**About this tool** + +**BCFtools merge**: Merge multiple VCF/BCF files from non-overlapping sample sets to create one multi-sample file. For example, when merging file A.vcf.gz containing samples S1, S2 and S3 and file B.vcf.gz containing samples S3 and S4, the output file will contain four samples named S1, S2, S3, 2:S3 and S4. + +Please see https://samtools.github.io/bcftools/bcftools.html for more info on options. + +