#!/usr/bin/env cwl-runner

cwlVersion: v1.0
class: Workflow
label: "sequence (bam or fastqs) to trimmed fastqs and HISAT alignments"
requirements:
    - class: MultipleInputFeatureRequirement
    - class: SubworkflowFeatureRequirement
    - class: InlineJavascriptRequirement
    - class: StepInputExpressionRequirement
    - class: SchemaDefRequirement
      types:
          - $import: ../types/sequence_data.yml

inputs:
    unaligned:
        type: ../types/sequence_data.yml#sequence_data
    adapters:
        type: File
    adapter_trim_end:
        type: string
    adapter_min_overlap:
        type: int
    max_uncalled:
        type: int
    min_readlength:
        type: int
    read_group_id:
        type: string
    read_group_fields:
        type:
            type: array
            items: string
    reference_index:
        type: File
        secondaryFiles: [".1.ht2", ".2.ht2", ".3.ht2", ".4.ht2", ".5.ht2", ".6.ht2", ".7.ht2", ".8.ht2"]
    strand:
        type:
          - "null"
          - type: enum
            symbols: ["first", "second", "unstranded"]
    unzip_fastqs:
        type: boolean?
outputs:
    fastqs:
        type: File[]
        outputSource: trim_fastq/fastqs
    aligned_bam:
        type: File
        outputSource: hisat2_align/aligned_bam
steps:
    sequence_to_fastq:
        run: ../tools/sequence_to_fastq.cwl
        in: 
            bam:
                source: unaligned
                valueFrom: "$(self.sequence.hasOwnProperty('bam')? self.sequence.bam : null)"
            fastq1:
                source: unaligned
                valueFrom: "$(self.sequence.hasOwnProperty('fastq1')? self.sequence.fastq1 : null)"
            fastq2:
                source: unaligned
                valueFrom: "$(self.sequence.hasOwnProperty('fastq2')? self.sequence.fastq2 : null)"
            unzip_fastqs: unzip_fastqs
        out:
            [fastq1, fastq2]
    trim_fastq:
        run: ../tools/trim_fastq.cwl
        in:
            reads1: sequence_to_fastq/fastq1
            reads2: sequence_to_fastq/fastq2
            adapters: adapters
            adapter_trim_end: adapter_trim_end
            adapter_min_overlap: adapter_min_overlap
            max_uncalled: max_uncalled
            min_readlength: min_readlength
        out:
            [fastqs]
    hisat2_align:
        run: ../tools/hisat2_align.cwl
        in:
            reference_index: reference_index
            fastq1: 
                source: trim_fastq/fastqs
                valueFrom: $(self[0])
            fastq2: 
                source: trim_fastq/fastqs
                valueFrom: $(self[1])
            read_group_id: read_group_id
            read_group_fields: read_group_fields
            strand: strand
        out:
            [aligned_bam]