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!