Extracting A Rs* Genotype Using Vcf Tools
Does any one know how to extract a SNP genotype from 1000 Genome data easily? VCF tools of GATK?
• 4,673 views
•
link
2 answers
I think you want:
vcftools --vcf vcf_file --snp rs_id --recode --keep-INFO-all --out output_prefix
• 2 views
•
link
awk '/^#/||$2=="rs123456"' in.vcf > out.vcf
• 0 views
•
link
Log in to answer this question.
what is gatk? Did you mean gawk?