vcf文件去除非变异的基因型(use GATK exclude nonvariant in vcf format,0|0,0/0)

对于某些特殊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

  参考链接:

https://www.biostars.org/p/203809/

https://software.broadinstitute.org/gatk/documentation/tooldocs/3.8-0/org_broadinstitute_gatk_tools_walkers_variantutils_SelectVariants.php

posted @ 2018-06-22 14:37  橙子牛奶糖  阅读(1390)  评论(0编辑  收藏  举报