NCBI Human Variation Sets VCF validation

Question

[Jirapong]

The reference VCF file that available to download at http://www.ncbi.nlm.nih.gov/variation/docs/human_variation_vcf/ seems to have format issue.

This is happen when I try to run GATK BaseRecalibrator command. The error show that there is a duplicate allele added to VariantContext. I start to looks at the reference VCF file and found that there is duplicate Alternative Allele and also sometimes REF is showing in ALT column. I download vcftools to validate and get following waring

[tmp]$ vcf-validator homo_sapiens_GRCh37.vcf
1:2886090 .. REF allele listed in the ALT field??
1:4095845 .. Could not parse the allele(s) [AG], first base does not match the reference.
1:8121167 .. The alleles not unique: CAAT
1:8121167 .. The alleles not unique: CAAT
1:9127042 .. The alleles not unique: TAA
1:11408760 .. The alleles not unique: CTATGTATG
1:13177471 .. The alleles not unique: CTT
1:13894414 .. REF allele listed in the ALT field??
1:15015689 .. REF allele listed in the ALT field??

Is there other source where I can get a better Variation Sets file?

Answer 1

[iraun]

GATK provides variation sets files for human. I would recommend you to use them in order to prevent this kind of errors.

• Some interesting information about available VCF files and databases in GATK: http://gatkforums.broadinstitute.org/discussion/1213/whats-in-the-resource-bundle-and-how-can-i-get-it
• Link to access to ftp to download data: https://www.broadinstitute.org/gatk/guide/article?id=1215

Hope it helps.

[Jirapong]

@airan - Thank you very much. I able to get it from the second link ftp.broadinstitute.org/bundle/2.8/hg19/dbsnp_138.hg19.vcf.gz

[iraun]

Glad to help :)