GWAS前面的部分,打包在这里
链接:https://pan.baidu.com/s/1cp4DY2Q0rYuNPT7skp1Dww
提取码:zwbi
然后找到一个基因,PMEL
plink –file BB –chr 34 –from-bp 1339431 –to-bp 1343416 –make-bed –out output_BB –chr-set 39
拆分
bcftools view -s IDGGB_RJF01,IDGGB_RJF02,IDGGB_RJF03 GRCg6a.928.SNP_2.vcf -Ov -o species2.vcf