VCFTOOLS 的一些常用参数

转换为PLINK格式:PLINK1.9

vcftools --gzvcf xxx.vcf --plink --out xxx

plink --vcf xxx.vcf --double-id A --out step1
plink --tfile step1 --recode --out step2
plink --file step2 --geno 0.1 --maf 0.05 --hwe 0.000001 --make-bed --out step3

VCFTOOLS的常用处理:VCF Manual

vcftools --gzvcf xxx.vcf.gz --keep id.txt --recode --stdout | gzip -c > yyy.vcf.gz
vcftools --gzvcf xxx.vcf.gz --positions pos.list --recode --out yyy
vcftools --gzvcf xxx.vcf.gz --chr 1 --from-bp a --to-bp b --recode --recode-INFO-all --out yyy
vcftools --gzvcf xxx.vcf.gz --freq --out xxx 

猜你喜欢

转载自blog.csdn.net/weixin_42284584/article/details/89154773