rna-seq-tophat/tasks/tophat2.wdl

task tophat2 {
   File gtf
   File genome_directory
   String idx_prefix
	File read1
	File read2
   String baseout
   String docker
   String cluster

   command {
      tophat2 -p 24 -o ${baseout} -G ${gtf} --library-type fr-unstranded --solexa-quals ${genome_directory}/${idx_prefix} ${read1} ${read2}
   }
   
   runtime { 
      docker: docker 
		cluster: cluster
		systemDisk: "cloud_ssd 40"
		dataDisk: "cloud_ssd 200 /cromwell_root/"
   }

   output {
      File accepted_hits = "${baseout}/${baseout}_accepted_hits.bam"
      File unmapped_bam = "${baseout}/${baseout}_unmapped.bam"
   }
}