rnaseq-engineering/tasks/qualimap.wdl

task qualimap {
	File bam
	File bai
	File gtf
	String bamname = basename(bam, ".sorted.bam")
	String docker
	String cluster
	String disk_size

	command <<<
		set -o pipefail
		set -e
		nt=$(nproc)
		/opt/qualimap/qualimap rnaseq -bam ${bam} -outformat HTML -outdir ${bamname}_rnaseq -gtf ${gtf} -pe --java-mem-size=10G
		tar -zcvf ${bamname}_rnaseq_qualimap.zip ${bamname}_rnaseq
		/opt/qualimap/qualimap bamqc -bam ${bam} -outformat PDF:HTML -nt $nt -outdir ${bamname}_bamqc --java-mem-size=32G 
		tar -zcvf ${bamname}_bamqc_qualimap.zip ${bamname}_bamqc
	>>>

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

	output {
		File rnaseq_zip = "${bamname}_rnaseq_qualimap.zip"
		File bamqc_zip = "${bamname}_bamqc_qualimap.zip"
	}
}