rna-seq-tophat/tasks/trimmomatic.wdl

task trimmomatic {
	File read1
	File read2
	File adapter
	String baseout
	String docker
	String cluster

	command {
  	trimmomatic PE -threads 20 -phred33 ${read1} ${read2}  -baseout ${baseout} ILLUMINACLIP:${adapter}:2:30:10:1:true HEADCROP:10 LEADING:10 TRAILING:10 SLIDINGWINDOW:4:15 MINLEN:36
	}

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

	output {
		File read_1p = "${baseout}/${baseout}_1P.fq.gz"
		File read_2p = "${baseout}/${baseout}_2P.fq.gz"
	}
}