Question: error using bqsr and incompatible contigs
0
gravatar for bioguy24
8 months ago by
bioguy24180
Chicago
bioguy24180 wrote:

Using the below command in GATK 4.0.6.0, I am getting:

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

Code:

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 • 421 views
ADD COMMENTlink modified 3 months ago by Biostar ♦♦ 20 • written 8 months ago by bioguy24180
1

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 8 months ago by RamRS21k
1

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 8 months ago by h.mon24k

Thank you both very much :).

ADD REPLYlink written 8 months ago by bioguy24180
0
gravatar for h.mon
8 months ago by
h.mon24k
Brazil
h.mon24k 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 8 months ago by h.mon24k
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: 1154 users visited in the last hour