retrieving from ExAC in VCF format
2
1
Entering edit mode
4.0 years ago

Hello! Is it possible to retrieve information about variants of a single gene from the ExAC Browser in VCF format? I can only see the option to "Export table to CSV".

ExAC VCF • 2.6k views
ADD COMMENT
0
Entering edit mode

thank you for your advice!

ADD REPLY
6
Entering edit mode
4.0 years ago
Vivek ★ 2.5k

You can use tabix and query directly from the FTP:

tabix -h ftp://ftp.broadinstitute.org/pub/ExAC_release/current/ExAC.r0.3.1.sites.vep.vcf.gz 2:39967768-39967768 > exported.vcf
ADD COMMENT
1
Entering edit mode

Thanks a lot! It works perfectly

ADD REPLY
1
Entering edit mode

Hello, is it possible to generate vcf files with tabix and exclude the 7601 samples deriving from the TCGA cohort in ExAC?

ADD REPLY
2
Entering edit mode

There are tools like vcftools, GATK or bcftools that can achieve this:

This thread might be helpful: Extract subset of samples from multigenome vcf file

ADD REPLY
0
Entering edit mode

I now have a VCF file and I would like to keep only the non-Finn European (NFE). I don't have a list of samples to keep, but I think this can be done by "INFO FIELD FILTERING" with VCFtools. In fact, I read lines like this: "##INFO=<id=an_nfe ..="" which="" are="" what="" i="" need.="" please="" let="" me="" know.="" thanks<="" p="">

ADD REPLY
1
Entering edit mode

If its in the info fields, you can use GATK's SelectVariants utility. Here are some examples on using JEXL expressions to filter your VCF file: https://gatkforums.broadinstitute.org/gatk/discussion/1255/using-jexl-to-apply-hard-filters-or-select-variants-based-on-annotation-values

ADD REPLY
0
Entering edit mode

Apologies, but maybe I have done some confusion above: I have looked at your link, but I'm afraid this is not very helpful for what I have to do. Basically I must keep only one subpopulation in my VCF file and in the info fields I have information about the populations. Now suppose I only want to keep one population in my vcf file (non Finn Europeans, NFE): is there an easy way to remove all the rest and keep only non finn europeans (maybe with vcftools)?

ADD REPLY
0
Entering edit mode
4.0 years ago
daniko07 ▴ 10

I use the following command to annotate vcf files: table_annovar.pl example/ex2.vcf humandb/ -buildver hg19 -out myanno -remove -protocol refGene,cytoBand,exac03,avsnp147,dbnsfp30a -operation g,r,f,f,f -nastring . -vcfinput

but unfortunately I see the COSMIC and CLINVAR annotations are missing. How can I include them? Please let me know. Thanks

ADD COMMENT

Login before adding your answer.

Traffic: 2590 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