对于某些特殊vcf,想去除没有变异的基因型(主要形式为0|0或者0/0),则需要用到GATK的--excludeNonVariants参数,命令行如下:
java -Xmx8g -jar GenomeAnalysisTK.jar -T SelectVariants -R GRCh37.fa -V 1000Genomes.vcf --excludeNonVariants -o 1000Genomes_onlyvariants.vcf
参考链接:
/p/203809/
/gatk/documentation/tooldocs/3.8-0/org_broadinstitute_gatk_tools_walkers_variantutils_SelectVariants.php
如果觉得《vcf文件去除非变异的基因型(use GATK exclude nonvariant in vcf format 0|0 0/0)》对你有帮助,请点赞、收藏,并留下你的观点哦!