seq 比对

合并VCF Merge VCF 的方法:

2022-02-17  本文已影响0人  11的雾

Merge VCF 的方法:

1. gatk mergVcf

gatk MergeVcfs -I M11189W.chr21.common.g.vcf.gz -I N11189.chr21.common.g.vcf.gz -O merged.picard.vcf.gz

Error: has sample entries that don't match the other files.

参考GATK mergevcf模块,注意两个vcf要相互match

2. gatk GatherVcfs

gatk GatherVcfs I=35.vcf I=36.vcf O=all.vcf

3. bcftools

bcftools merge --merge all file1.vcf.gz file2.vcf.gz -O v -o merged.vcf.gz

4. python: bioinfokit package

from bioinfokit.analys import marker

# concatenate VCF files. You can provide multiple VCF files separated by comma.

marker.concatvcf("file_1.vcf,file_2.vcf,file_3.vcf,file_4.vcf")# merged VCF files will be stored in same directory (concat_vcf.vcf)

5. snpsift

java -jar SnpSift.jar split -j *.vcf > combined.vcf

上一篇下一篇

猜你喜欢

热点阅读