Question: error using bqsr and incompatible contigs
gravatar for cmccabe
10 weeks ago by
cmccabe170 wrote:

Using the below command in GATK, I am getting:

A USER ERROR has occurred: Input files reference and features have incompatible contigs: No overlapping contigs found.


gatk BaseRecalibrator \ 
    -I NA12878_mark_duplicates.bam \ 
    -R /home/cmccabe/Desktop/fastq/resources/ucsc.hg19.fasta \
     --known-sites /home/cmccabe/Desktop/fastq/resources/dbsnp_138.b37.vcf \ 
     --known-sites /home/cmccabe/Desktop/fastq/resources/Mills_and_1000G_gold_standard.indels.hg19.vcf \
     -O /home/cmccabe/Desktop/fastq/NA12878_recal_data.table

I am using the resource bundle from the ftp as well as the fasta. Is there a way to ignore these errors or am I doing something wrong? I also aligned with the same fasta using bwa-mem. Thank you :).

bqsr gatk • 173 views
ADD COMMENTlink modified 10 weeks ago by h.mon20k • written 10 weeks ago by cmccabe170

Can you check the dbsnp file and ensure the contig names match those in the ucsc.hg19.fasta? I think the problem might lie there.

ADD REPLYlink written 10 weeks ago by RamRS18k

The warnings is hinting you have different chromosome names somewhere. Examine the fasta, bam and vcf and check if they all use chr1, chr2, chr3, and so on; or they all use 1, 2, 3, and so on.

ADD REPLYlink written 10 weeks ago by h.mon20k

Thank you both very much :).

ADD REPLYlink written 10 weeks ago by cmccabe170
gravatar for h.mon
10 weeks ago by
h.mon20k wrote:

Upon closer inspection, I see you are mixing hg19 and b37 references. Although they are largely the same, they are not identical, and there are more differences than just chromosome naming convention:

  • differences in mtDNA sequence
  • different REF bases at some SNPs

You should use the exactly same reference for all resources (reference genome for mapping and reference vcfs).

Further discussions:

What's the difference between b37 and hg19 resources?

problems in hg19 and b37 compatibility

ADD COMMENTlink written 10 weeks ago by h.mon20k
Please log in to add an answer.


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