diff --git a/CHANGELOG.md b/CHANGELOG.md
index b7a8741f309cd7fa4f12fc5043f7817e94594411..836af4ac2fc7b81d157813a97fcaea4c2a013269 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -12,6 +12,7 @@ that users understand how the changes affect the new version.
 version 5.0.0-dev
 ---------------------------
 + Complete `parameter_meta` for tasks missing the outputs.
++ DeepVariant: Add an optional input for the gvcf index.
 + Samtools: `Sort` task now has `threads` in runtime instead of `1`.
 + Picard: Add parameter_meta to `SortSam`.
 + pbmm2: Add parameter_meta for `sample`.
diff --git a/deepvariant.wdl b/deepvariant.wdl
index f71a1c88521ed264aa976e7c6e396a774e1dc72d..618200aa8df9c3b4be9700d323d48eaf300f5eb6 100644
--- a/deepvariant.wdl
+++ b/deepvariant.wdl
@@ -28,11 +28,11 @@ task RunDeepVariant {
         File inputBamIndex
         String modelType
         String outputVcf
-
         String? postprocessVariantsExtraArgs
         File? customizedModel
         Int? numShards
         String? outputGVcf
+        String? outputGVcfIndex
         File? regions
         String? sampleName
         Boolean? VCFStatsReport = true
@@ -44,6 +44,7 @@ task RunDeepVariant {
 
     command {
         set -e
+
         /opt/deepvariant/bin/run_deepvariant \
         --ref ~{referenceFasta} \
         --reads ~{inputBam} \
@@ -66,10 +67,10 @@ task RunDeepVariant {
 
     output {
         File outputVCF = outputVcf
-        File outputVCFIndex = outputVCF + ".tbi"
+        File outputVCFIndex = outputVcf + ".tbi"
         Array[File] outputVCFStatsReport = glob("*.visual_report.html")
         File? outputGVCF = outputGVcf
-        File? outputGVCFIndex = outputGVcf + ".tbi"
+        File? outputGVCFIndex = outputGVcfIndex
     }
 
     parameter_meta {
@@ -84,6 +85,7 @@ task RunDeepVariant {
         customizedModel: {description: "A path to a model checkpoint to load for the `call_variants` step. If not set, the default for each --model_type will be used.", category: "advanced"}
         numShards: {description: "Number of shards for make_examples step.", category: "common"}
         outputGVcf: {description: "Path where we should write gVCF file.", category: "common"}
+        outputGVcfIndex: {description: "Path to where the gVCF index file will be written. This is needed as a workaround, set it to outputGVcf+.tbi.", category: "common"}
         regions: {description: "List of regions we want to process, in BED/BEDPE format.", category: "advanced"}
         sampleName: {description: "Sample name to use instead of the sample name from the input reads BAM (SM tag in the header).", category: "common"}
         VCFStatsReport: {description: "Output a visual report (HTML) of statistics about the output VCF.", category: "common"}