Thanks a ton Nandini ... it works :)
How to extract SNPs from vcf file based on Population
Dear Friends,
My vcf file has SNPs available for different population(Africa, America, Europe,East Asia and South Asia ). I want to extract the data for Europe and East Asia together . Kindly let me know the possible ways.
Thanks in Advance
• 6,064 views
•
link
1 answer
You can do this easily using vcftools, GATK tools, plinkseq etc.
you first have to generate a text file with the list of samples that form the population of your choice, let's say "population_of_interest.txt" Then,
vcf-subset -e -c population_of_interest.txt input.vcf > output.vcf
or
vcftools --vcf input.vcf --keep population_of_interest.txt --recode > output.vcf
• 0 views
•
link
• 0 views
•
link
This code works fine when i run for one chromosome at a time. But, I want to extract SNPs for all chromosomes together ,please let me know if ithere is any other option ?
• 1 views
•
link
Log in to answer this question.