I'm trying to get a vcf with the invariant sites i.e. select out ONLY the homozygous references sites in a multi-sample vcf (samples are either all 0/0 or a mix of 0/0 and ./.). I do not want the sites where the reference is N or all the sites are ./. (before filtering)
Is the only way to do this:
VariantFiltration with a filter QUAL > -1 (to remove the reference N sites)
Then
selectVariants with --selectTypeToInclude NO_VARIATION