Question: Human GRCh38 and dbSNP VCF and GATK
0
gravatar for win
2.9 years ago by
win810
India
win810 wrote:

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 • 3.2k views
ADD COMMENTlink modified 2.9 years ago by Carlos Borroto1.8k • written 2.9 years ago by win810
2
gravatar for Pierre Lindenbaum
2.9 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum119k wrote:

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 COMMENTlink written 2.9 years ago by Pierre Lindenbaum119k
1
gravatar for Carlos Borroto
2.9 years ago by
Carlos Borroto1.8k
Washington Metropolitan Area
Carlos Borroto1.8k wrote:

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 COMMENTlink written 2.9 years ago by Carlos Borroto1.8k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1361 users visited in the last hour