Rare variants from annotated vcf file

[rohan_1925]

Hi,

I have annotated my vcf file of 20 samples from GATK's Unified genotyper using the following steps.

Unified genotyper->Variantrecalibration->Applyrecalibration->VariantAnnotator

My question is how should I proceed if I have to select rare variants (MAF<1%) for the candidate genes that I have,for each of these 20 samples?

[m_two]

You may want to download the dbSNP VCF file containing the 1000 Genome's project minor allele frequencies.
ftp://ftp.ncbi.nlm.nih.gov/snp/organ...uman_9606/VCF/

You can use joinx to process the data:

https://github.com/genome/joinx

##INFO=<ID=CAF,Number=.,Type=String,Description="An ordered, comma delimited list of allele frequencies based on 1000Genomes, starting with the reference allele followed by alternate alleles as ordered in the ALT column. Where a 1000Genomes alternate allele is not in the dbSNPs alternate allele set, the allele is added to the ALT column. The minor allele is the second largest value in the list, and was previuosly reported in VCF as the GMAF. This is the GMAF reported on the RefSNP and EntrezSNP pages and VariationReporter">

[m_two]

In addition, this may be useful: http://evs.gs.washington.edu/EVS/