DRAGEN gVCF files and joint calling GATK-GenomicsDB
REQUIRED for all errors and issues: a) GATK version used:
4.5
b) Exact command used:
gatk GenomicsDBImport [GVCFS] --geomicsdb-workspace-path test1 --intervals chr20
gatk GenotypeGVCFs -R references/hs38DH.fa -V gendb://test1 -O test1.vcf
c) Entire program log: The tools ran accordingly, but the VCF output after joint had no variants.
My question is: Does gVCF files from DRAGEN can be used for the joint calling using the GenomicsDB workflow?
• 1,025 views
•
link
updated
by
DBScan •
•
written
by
brunomiwa •
•
written
by
brunomiwa •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
GenomicsDBImport completes successfully, but SelectVariants/GenotypeGVCFs hangs at “Starting traver…
written by chiyong4783 •
-
Error while doing joint genotype calling using GATK
written by ngseq •
-
GATK GenomicsDBImport update database error using snakemake
written by Peter Chung
-
GenotypeGVCF too many genotypes from pooled samples
written by Vic •
-
GATK GenomicsDBImport does not support GVCFs with MNPs
written by Rodolfo Adrián •
-
GATK GenomicsDBImport for consolidate GVCF
written by mscpythonstudy •
-
How to call denovo variants from multiple trios cohort?
written by DareDevil
-
Dealing with GATK illegal character
written by MAPK
-
Error in GATK GenomicsDBImport
written by raf.marcondes •
-
Some questions about call variants for mtDNA using GATK HaplotypeCaller
written by MatthewP
I think I did it once, but why would you use GATK with DRAGEN gVCF files? If you have a DRAGEN, the iterative gVCF Genotyper does a much faster and easier job.