Intersect exome region from whole genome gvcf
0
0
Entering edit mode
3.9 years ago
MAPK ★ 2.1k

I have been trying to get a joint genotype for exome bed region from whole genome gvcf. However, this is really slow and I am not sure if this is the right approach to do it. I have more than 600 samples and this has been running for more than 10 days. Can someone please clarify?

module load java
/gatk-4.1.7.0/gatk --java-options "-Xms4G -Xmx32G -DGATK_STACKTRACE_ON_USER_EXCEPTION=true" \
 GenomicsDBImport \
--genomicsdb-workspace-path /scratch/chr19 \
--batch-size 50 \
-L /scratch/chr19.exome.bed \
--sample-name-map /scratch/AllSamplesGvcfChr19.list \
--reader-threads 16
gatk GenomicsDBImport • 974 views
ADD COMMENT
1
Entering edit mode

I got the same problem, at the end I moved back to gatk/CombineGVCF ...

ADD REPLY
0
Entering edit mode

Was it faster? Did you do GenotypeGVCFs afterwards?

ADD REPLY
1
Entering edit mode

yes and yes

ADD REPLY

Login before adding your answer.

Traffic: 1126 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6