somatic_SNVs_pipeline/tasks/bcftools_concat.wdl

task bcftools {
  
  File ref_dir
  String fasta
  File vcf_indels
  File vcf_snvs
  String sample_id
  String docker
  String cluster_config
  String disk_size

  command <<<
    set -o pipefail
    set -e
    nt=$(nproc)
    
    bcftools concat -a ${vcf_indels} ${vcf_snvs} -Oz -o ${sample_id}.strelka.vcf.gz
    bcftools norm -m -both ${sample_id}.strelka.vcf.gz -O z -o ${sample_id}.strelka.norm.vcf.gz
    bcftools view -O z -f 'PASS' ${sample_id}.strelka.norm.vcf.gz -o ${sample_id}.strelka.PASS.vcf.gz
    bcftools index -t ${sample_id}.strelka.PASS.vcf.gz

  >>>
  
  runtime {
    docker: docker
    cluster: cluster_config
    systemDisk: "cloud_ssd 40"
    dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
  }
  
  output {
    File concat_vcf = "${sample_id}.strelka.PASS.vcf.gz"
    File concat_vcf_index = "${sample_id}.strelka.PASS.vcf.gz.tbi"
  }
}