bactools合并plink格式不同染色体的数据

 

bactools合并plink格式不同染色体的数据

 

001、测试数据

[root@PC1 test]# ls
chr_10.vcf.gz      chr_2.vcf.gz      chr_4.vcf.gz      chr_6.vcf.gz      chr_8.vcf.gz
chr_10.vcf.gz.tbi  chr_2.vcf.gz.tbi  chr_4.vcf.gz.tbi  chr_6.vcf.gz.tbi  chr_8.vcf.gz.tbi
chr_1.vcf.gz       chr_3.vcf.gz      chr_5.vcf.gz      chr_7.vcf.gz      chr_9.vcf.gz
chr_1.vcf.gz.tbi   chr_3.vcf.gz.tbi  chr_5.vcf.gz.tbi  chr_7.vcf.gz.tbi  chr_9.vcf.gz.tbi
[root@PC1 test]# bcftools concat -O z -o merged.vcf.gz *.vcf.gz    ## 合并不同染色体的数据
Checking the headers and starting positions of 10 files
Concatenating chr_10.vcf.gz     0.009121 seconds
Concatenating chr_1.vcf.gz      0.027479 seconds
Concatenating chr_2.vcf.gz      0.028607 seconds
Concatenating chr_3.vcf.gz      0.025100 seconds
Concatenating chr_4.vcf.gz      0.014592 seconds
Concatenating chr_5.vcf.gz      0.011217 seconds
Concatenating chr_6.vcf.gz      0.010887 seconds
Concatenating chr_7.vcf.gz      0.011350 seconds
Concatenating chr_8.vcf.gz      0.011012 seconds
Concatenating chr_9.vcf.gz      0.009078 seconds
[root@PC1 test]# ls    ## 合并结果
chr_10.vcf.gz      chr_2.vcf.gz      chr_4.vcf.gz      chr_6.vcf.gz      chr_8.vcf.gz      merged.vcf.gz
chr_10.vcf.gz.tbi  chr_2.vcf.gz.tbi  chr_4.vcf.gz.tbi  chr_6.vcf.gz.tbi  chr_8.vcf.gz.tbi
chr_1.vcf.gz       chr_3.vcf.gz      chr_5.vcf.gz      chr_7.vcf.gz      chr_9.vcf.gz
chr_1.vcf.gz.tbi   chr_3.vcf.gz.tbi  chr_5.vcf.gz.tbi  chr_7.vcf.gz.tbi  chr_9.vcf.gz.tbi

 

 。