battenberg_only/tasks/battenberg.wdl

task battenberg {
    String sample_id
    File tumor_bam
    File tumor_bam_bai
    String normal_id
    File normal_bam
    File normal_bam_bai

    String docker
    String cluster
    String disk_size

    command <<<

        set -o pipefail
        set -e
        cp ${normal_bam} ./${normal_id}.bam
        cp ${tumor_bam} ./${sample_id}.bam
        
        cp ${normal_bam_bai} ./${normal_id}.bam.bai
        cp ${tumor_bam_bai} ./${sample_id}.bam.bai
      
        cp /opt/battenberg_reference/battenberg_wgs.R ./
        /opt/conda/envs/batternberg/bin/Rscript ./battenberg_wgs.R -t ${sample_id} -n ${normal_id} --nb ./${normal_id}.bam --tb  ./${sample_id}.bam --cpu 4 --sex female -o ./
        mkdir ${sample_id}_chr_png
        cp ${sample_id}_subclones_chr*.png ./${sample_id}_chr_png/
        tar zcvf ./${sample_id}_chr_png.tgz ./${sample_id}_chr_png/

    >>>

    runtime {
        docker : docker
        cluster: cluster
        systemDisk: "cloud_ssd 100"
        dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
    }

    output {
        File subclones = "${sample_id}_subclones.txt"
        File rho = "${sample_id}_rho_and_psi.txt"
        File average_png = "${sample_id}_BattenbergProfile_average.png"
        File subclones_png = "${sample_id}_BattenbergProfile_subclones.png"
        File distance_png = "${sample_id}_distance.png"
        File tumor_png = "${sample_id}.tumour.png"
        File germline_png = "${sample_id}.germline.png"
        File coverage_png = "${sample_id}_coverage.png"
        File alleleratio_png = "${sample_id}_alleleratio.png"
        File chr_png = "${sample_id}_chr_png.tgz"
    }

}