diff --git a/somaticseq.wdl b/somaticseq.wdl index 4a8b860d71f6e8750f482b91419fd619a21968d1..349c7d2f4b9042167990a1b26d7b7cbb18c140b3 100644 --- a/somaticseq.wdl +++ b/somaticseq.wdl @@ -280,9 +280,14 @@ task ModifyStrelka { } command { + set -e -o pipefail + ~{installDir}/modify_Strelka.py \ -infile ~{strelkaVCF} \ - -outfile ~{outputVCFName} + -outfile "modified_strelka.vcf" + + first_FORMAT_line_num=$(grep -n -m 1 '##FORMAT' "modified_strelka.vcf" | cut -d : -f 1) + sed "$first_FORMAT_line_num"'i##FORMAT=<ID=GT,Number=1,Type=String,Description="Genotype">' "modified_strelka.vcf" > ~{outputVCFName} } output {