Hi all -
I have a VCF file containing genotypes for 12 samples at variant positions only. I'd like to convert those genotypes with GQ scores under a threshold to "no call" (i.e. ./.)
Does anyone know a utility to do this? My best thought at this point is to script it myself.
Thanks in advance!
I have a VCF file containing genotypes for 12 samples at variant positions only. I'd like to convert those genotypes with GQ scores under a threshold to "no call" (i.e. ./.)
Does anyone know a utility to do this? My best thought at this point is to script it myself.
Thanks in advance!