rnaseq-engineering/tasks/hisat2.wdl

task hisat2 {
  File idx
  File read_1P
  File read_2P
  String sample_name=sub(basename(read_1P),"\\_R1.\\S+$", "")
  String idx_prefix
  String docker
  String cluster
  String disk_size
   

  command {
    nt=$(nproc)
    hisat2 -t -p $nt -x ${idx}/${idx_prefix} -1 ${read_1P} -2 ${read_2P} -S ${sample_name}.sam  --un-conc-gz ${sample_name}_un.fq.gz
  }
   
  runtime { 
    docker: docker
    cluster: cluster
    systemDisk: "cloud_ssd 40"
    dataDisk: "cloud_ssd " + disk_size + " /cromwell_root/"
  }

  output {
    File sam=sample_name + ".sam"
    File unmapread_1p=sample_name + "_un.fq.1.gz"
    File unmapread_2p=sample_name + "_un.fq.2.gz"
  }
}