将 Control 和 Case vcf 文件转换为 plink 1.9.0 可以处理的格式

将 Control 和 Case vcf 文件转换为 plink 1.9.0 可以处理的格式

我有两个 vcf 文件,一个包含有 TP53 突变的样本,另一个包含没有 TP53 突变的样本。没有突变的vcf是对照,有突变的vcf是情况。如何将这两个文件组合成可以由 plink 1.9.0 中的 --assoc 函数处理的文件?

答案1

我不确定我是否明白你想做什么,但我认为你需要的是合并两个文件。您需要使用 tabix 对它们进行索引,然后尝试:

bcftools merge file1.vcf.gz file2.vcf.gz -o merged_files.vcf.gz

然后我会尝试:

plink --vcf merged_files.vcf.gz --assoc [params]

如果这不能回答您的问题,请尝试概述您的输入是什么以及您的输出 VCF 应该是什么样子。

相关内容