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

[RFE] haplotypecaller: support for custom GVCF GQ slices #52

Open
ntm opened this issue Jun 30, 2021 · 1 comment
Open

[RFE] haplotypecaller: support for custom GVCF GQ slices #52

ntm opened this issue Jun 30, 2021 · 1 comment

Comments

@ntm
Copy link

ntm commented Jun 30, 2021

Hi,

GATK4 allows to specify custom boundaries for the GQ slices used to group consecutive non-variant positions into a non-variant block, with --gvcf-gq-bands . This feature is quite useful IMO, as the resulting GVCFs can be much smaller (and much faster to work with downstream). The only "loss" is fine-grained resolution of GQ slices for non-variant blocks, but only in the GQ ranges that the user decided he doesn't care about.
It would be great if a similar feature could be implemented in elPrep!

Regards,
Nicolas

@matthdsm
Copy link
Contributor

This seems like something that can be done using GATK ReblockGVCF, which will be part of the updated Broad WGS guidelines

M

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

No branches or pull requests

2 participants