diff --git a/gatk.wdl b/gatk.wdl index 997266c838c7d35afd4d1afd77573a0456d07043..ab1c215b390e4d67a1e4122aa1d6d556aa4e0389 100644 --- a/gatk.wdl +++ b/gatk.wdl @@ -194,3 +194,28 @@ task CombineGVCFs { File output_gvcf_index = output_basename + ".vcf.gz.tbi" } } + +task SplitNCigarReads { + String? preCommand + + File input_bam + File ref_fasta + File ref_fasta_index + File ref_dict + 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" + } +}