Adam Auton
-
2013-09-19
- assigned_to: Petr Danecek
vcf-merge fails when trying to merge blank vcf files. Trying to merge VCF fils that contain only a header, vcf-merge throws the following error:
close failed: test.vcf.gz
main::error('close failed: test.vcf.gz\x{a}') called at my path/vcf-merge line 543
Should it not just create a file with only a header containing an entry for all merged VCF samples (which is the behaviour exhibited by 'htscmd vcfmerge')?