Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

merge command collapsing non-identical variants to single record #2063

Open
Logan-Wallace opened this issue Dec 21, 2023 · 1 comment
Open

Comments

@Logan-Wallace
Copy link

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
@pd3
Copy link
Member

pd3 commented Dec 27, 2023

How do you distinguish the duplications, by their ID? Would --merge id work in such case?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

2 participants