基因型过滤标准:保留双等位基因(biallelic sites)以及次等位基因频率大于0.05的位点:

bcftools view --types snps -m 2 -M 2 -q 0.05:minor genotypes.chr22.vcf | bgzip -c > genotypes.chr22.vcf.gz

  • --types snps -m 2 -M 2 指的是保留双等位基因(biallelic sites);
  • 0.05:minor 指的是次等位基因频率大于0.05;
  • genotypes.chr22.vcf.gz 指的是过滤后的输出文件;

相关文章:

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