Entering edit mode
                    5.4 years ago
        MAPK
        
    
        ★
    
    2.1k
    I am tryiing to run GenomicsDBImport for exome region of chromosome 2 for whole genome gvcf (multiple samples provided as mylist2.list)
my bed file is in this format:
chr2.bed
chr2    41359   41765
chr2    45307   45927
chr2    45729   46454
mylist2.list
1_S    /544562991/1_S.chr2.g.vcf.gz
2_S    /544562991/2_S.chr2.g.vcf.gz
After running this command:
gatk --java-options "-Xms4G -Xmx32G -DGATK_STACKTRACE_ON_USER_EXCEPTION=true" \
 GenomicsDBImport \
--genomicsdb-workspace-path /home/exomes/SNPs/chr2 \
--batch-size 50 \
--intervals /08-JointCalling/SNPs/chr2.bed \
--sample-name-map /home/mylist2.list \
--tmp-dir=/data/tmp \
--reader-threads 5
I am getting this error below. Can someone please help me understand/resolve this error?
org.broadinstitute.hellbender.exceptions.GATKException: Cannot call query with different interval, expected:chr2:41360-41765 queried with: chr2:45308-46454
I remember at one time GATK GenomicsDBImport did not allow multiple intervals. It seems now this problem has been fixed. See: https://gatkforums.broadinstitute.org/gatk/discussion/13447/genomicsdbimport-intervals. Can you check your GATK version?
Yes that was version issue. It worked on GATK 4.1.7.0. Thanks!