If you happen to see a question you know the answer to, please do chime in and help your fellow community members. We appreciate your help!
Test-drive the GATK tools and Best Practices pipelines on Terra
Check out this blog post to learn how you can get started with GATK and try out the pipelines in preconfigured workspaces (with a user-friendly interface!) without having to install anything.
SelectVariants/getHomVarCount JEXL question
Had a JEXL syntax question that I was hoping to get some thoughts on. Let's say that I have a VCF file with 10 samples in it, S1 through S10. I would like to select the sites where among 5 specific samples, 3 are homozygous variant. So something like:
select ' vc.getHomVarCount(S1, S2, S3, S4, S5) >= 10 '
However this syntax doesn't work unfortunately. Is there a way to specify how to do something like this using vc.GetHomVarCount()?
Thanks as always!