bcftools icon indicating copy to clipboard operation
bcftools copied to clipboard

merge command collapsing non-identical variants to single record

Open Logan-Wallace opened this issue 1 year ago • 1 comments

Hi there, When using the merge command to combine VCFs from 60 individuals I am getting an output which includes non-identical variants being output to the same record without a way for me to parse how the allelic ratios for each individual accord to the collapsed variant calls. For example, multiple duplication events output to the same record because they share the same start location but there isn't a way to tell to which variant the allelic ratio for a given individual accords. The argument --collapse none would seem to solve this problem but isn't available for the merge command. Is there a method to stop non-identical variants from outputting to the same record in the merged file or could the --collapse option be added to the merge command?

Command - bcftools merge –file-list vcf_filenames -Oz -o merged_files.vcf.gz

Screenshot 2023-12-21 at 11 35 59 AM

Logan-Wallace avatar Dec 21 '23 20:12 Logan-Wallace