command <<< | command <<< | ||||
mv ${snv_gzvcf} ./${sample}.snv.vcf.gz | mv ${snv_gzvcf} ./${sample}.snv.vcf.gz | ||||
mv ${indel_gzvcf} ./${sample.indel.vcf.gz} | |||||
mv ${indel_gzvcf} ./${sample}.indel.vcf.gz | |||||
rtg index -f vcf ${sample}.snv.vcf.gz | rtg index -f vcf ${sample}.snv.vcf.gz | ||||
rtg index -f vcf ${sample}.indel.vcf.gz | rtg index -f vcf ${sample}.indel.vcf.gz | ||||
>>> | >>> | ||||
dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" | dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/" | ||||
} | } | ||||
output { | output { | ||||
File snv_gzvcf_2 = "${sample}.snv.vcf.gz" | |||||
File indel_gzvcf_2 = "${sample}.indel.vcf.gz" | |||||
File snv_gzvcf_index = "${sample}.snv.vcf.gz.tbi" | File snv_gzvcf_index = "${sample}.snv.vcf.gz.tbi" | ||||
File indel_gzvcf_index = "${sample}.indel.vcf.gz.tbi" | File indel_gzvcf_index = "${sample}.indel.vcf.gz.tbi" | ||||
} | } |
task merge { | task merge { | ||||
File snv_gzvcf | |||||
File indel_gzvcf | |||||
File snv_gzvcf_2 | |||||
File indel_gzvcf_2 | |||||
File snv_gzvcf_index | File snv_gzvcf_index | ||||
File indel_gzvcf_index | File indel_gzvcf_index | ||||
String sample | String sample | ||||
String disk_size | String disk_size | ||||
command <<< | command <<< | ||||
bcftools concat -a ${snv_gzvcf} ${indel_gzvcf} > ${sample}.vcf | |||||
bcftools concat -a ${snv_gzvcf_2} ${indel_gzvcf_2} > ${sample}.vcf | |||||
>>> | >>> | ||||
runtime { | runtime { |
call merge.merge as merge { | call merge.merge as merge { | ||||
input: | input: | ||||
snv_gzvcf_2=index.snv_gzvcf_2, | |||||
indel_gzvcf_2=index.indel_gzvcf_2, | |||||
snv_gzvcf=snv_gzvcf, | snv_gzvcf=snv_gzvcf, | ||||
indel_gzvcf=indel_gzvcf, | indel_gzvcf=indel_gzvcf, | ||||
snv_gzvcf_index=index.snv_gzvcf_index, | snv_gzvcf_index=index.snv_gzvcf_index, |