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 6 -G ${gtf} --library-type fr-unstranded --solexa-quals -o ${baseout} ${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"
   }
}