This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Filtering down 1000 Genomes data just for a set of SNPs

I have some SNP data, and I want to download 1000 Genomes VCF files, so that I can isolate out ancestral alleles for my sites. I know the data slicer tool can help me filter out populations and filter out individuals (to avoid having to download everything). I was wondering if anyone knows of a way to do this except for certain sites in the genome?

(My alternative would be get all 22 VCFs onto HPC and use a script to filter them out... and then use another script to isolate the ancestral allele)

vcf snp
• 2,724 views • link
updated by Ram Ram • written by devenvyas devenvyas

1 answer

if you only need the ancestral alleles, you can get them from the wgs sites file. a fairly simple way of getting them would be to use bcftools query:

bcftools query -f '%CHROM\t%POS\t%ID\t%AA\n' ftp://ftp.1000genomes.ebi.ac.uk/vol1/ftp/release/20130502/ALL.wgs.phase3_shapeit2_mvncall_integrated_v5b.20130502.sites.vcf.gz > 1000g.ancestral_alleles.txt

note that the ancestral allele on that last column comes in this format:

##INFO=<ID=AA,Number=1,Type=String,Description="Ancestral Allele. Format: AA|REF|ALT|IndelType. AA: Ancestral allele, REF:Reference Allele, ALT:Alternate Allele, IndelType:Type of Indel (REF, ALT and IndelType are only defined for indels)">

you can add any filtering option (by regions, by rs ids,...) to make this command as complex as needed. you'll be querying the file remotely and retrieving only what you need.

• 1 views • link
written by Jorge Amigo Jorge Amigo

Thanks! That worked!

• 0 views • link
written by devenvyas devenvyas

Log in to answer this question.

More posts like this