diff --git a/gatk.wdl b/gatk.wdl
index 997266c838c7d35afd4d1afd77573a0456d07043..2a2b2b0237e306356c435a10c1d893042576d4cc 100644
--- a/gatk.wdl
+++ b/gatk.wdl
@@ -194,3 +194,26 @@ task CombineGVCFs {
         File output_gvcf_index = output_basename + ".vcf.gz.tbi"
     }
 }
+
+task SplitNCigarReads {
+    String? preCommand
+
+    File input_bam
+    File ref_fasta
+    String output_bam
+    String gatk_jar
+
+    command {
+        set -e -o pipefail
+        ${preCommand}
+        java -Xms4G -jar ${gatk_jar} \
+        -I ${input_bam} \
+        -R ${ref_fasta} \
+        -O ${output_bam} # might have to be -o depending on GATK version
+    }
+
+    output {
+        File bam = output_bam
+        File bam_index = output_bam + ".bai"
+    }
+}