Question: error using bqsr and incompatible contigs
gravatar for cmccabe
5 months 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 • 296 views
ADD COMMENTlink modified 3 days ago by Biostar ♦♦ 20 • written 5 months 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 5 months ago by RamRS20k

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 5 months ago by h.mon22k

Thank you both very much :).

ADD REPLYlink written 5 months ago by cmccabe170
gravatar for h.mon
5 months ago by
h.mon22k 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 5 months ago by h.mon22k
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: 1341 users visited in the last hour