Extracting A Rs* Genotype Using Vcf Tools
2
1
Entering edit mode
12.6 years ago
Curiosity ▴ 130

Does any one know how to extract a SNP genotype from 1000 Genome data easily? VCF tools of GATK?

gatk • 4.0k views
ADD COMMENT
0
Entering edit mode

what is gatk? Did you mean gawk?

ADD REPLY
2
Entering edit mode
12.6 years ago
Adam ★ 1.0k

I think you want:

vcftools --vcf vcf_file --snp rs_id --recode --keep-INFO-all --out output_prefix

ADD COMMENT
1
Entering edit mode
12.6 years ago
lh3 33k

awk '/^#/||$2=="rs123456"' in.vcf > out.vcf

ADD COMMENT

Login before adding your answer.

Traffic: 2159 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6