Human GRCh38 and dbSNP VCF and GATK
2
0
Entering edit mode
8.0 years ago
win ▴ 980

Hi all, I used the GRCH38 reference file and aligned my FASTQ files to this, there were no issues here.

Then i used GATKs BaseRecalibration tool and used dbSNP144, at this step I get an error that states VCF and reference have incompatible contigs.

Could someone point where the issue might be or if there a GATK version of GRCh38 and dbSNP VCF.

Thanks

NGS • 5.6k views
ADD COMMENT
2
Entering edit mode
8.0 years ago

see http://gatkforums.broadinstitute.org/gatk/discussion/3843/plans-to-update-the-gatk-bundle

ftp://gsapubftp-anonymous@ftp.broadinstitute.org/bundle/hg38/hg38bundle/

(march 10) As we've mentioned previously, we are providing a beta version of the hg38 bundle but it is not yet officially supported, it is not guaranteed to be complete, and we are not taking requests for additional files. When the full official bundle is ready, we will make an announcement on the blog and its contents will be documented.

ADD COMMENT
1
Entering edit mode
8.0 years ago
Carlos Borroto ★ 2.1k

You problem seems to be you are using a dbSNP VCF with GRCh37 coordinates. The FTP location Pierre links in his answer has a GRCh38 version that should solve your problem.

However, you should know the current GATK version is not "alt" aware yet. While mappers like BWA are and can correctly map reads against a reference with alt contigs, current GATK will most likely fail to call variants in regions with alt contigs. If you must have GRCh38 and GATK at this time, you should strip the reference of any alt contigs before mapping. I believe this invalidates the main advantage of using GRCh38 and is the reason I'm holding the move to GRCh38.

ADD COMMENT

Login before adding your answer.

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