Entering edit mode
13.2 years ago
Curiosity
▴
130
Does any one know how to extract a SNP genotype from 1000 Genome data easily? VCF tools of GATK?
Does any one know how to extract a SNP genotype from 1000 Genome data easily? VCF tools of GATK?
I think you want:
vcftools --vcf vcf_file --snp rs_id --recode --keep-INFO-all --out output_prefix
awk '/^#/||$2=="rs123456"' in.vcf > out.vcf
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
what is gatk? Did you mean gawk?