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

Make bcf_remove_allele_set() set work on unpacked records #1358

Merged
merged 1 commit into from
Nov 10, 2021

Conversation

pd3
Copy link
Member

@pd3 pd3 commented Nov 10, 2021

Resolves #1353

vcfutils.c Outdated
@@ -256,10 +256,10 @@ int bcf_remove_allele_set(const bcf_hdr_t *header, bcf1_t *line, const struct kb
int *map = (int*) calloc(line->n_allele, sizeof(int));
uint8_t *dat = NULL;

bcf_unpack(line, BCF_UN_FMT);
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Should this be BCF_UN_ALL? Looking at bcf_unpack(), BCF_UN_FMT will only unpack the indiv data, but this needs to also access info and allele.

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Of course it should, duh

Copy link
Member Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

This is now replaced with the correct BCF_UN_ALL

@daviesrob daviesrob merged commit 0495748 into samtools:develop Nov 10, 2021
@pd3 pd3 deleted the fix-1353 branch February 28, 2023 13:39
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

vcfutils::bcf_remove_allele_set should unpack to BCF_UN_INFO
2 participants