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

相关文章:

  • 2021-05-19
  • 2022-12-23
  • 2022-12-23
  • 2022-02-06
  • 2021-12-25
  • 2022-12-23
  • 2021-04-29
猜你喜欢
  • 2021-11-09
  • 2021-11-09
  • 2022-12-23
  • 2022-01-19
  • 2021-04-07
  • 2022-01-08
  • 2022-12-23
相关资源
相似解决方案