#!/usr/bin/env cwl-runner
label: "Non-Coding Bacterial Genes"
cwlVersion: v1.0
class: Workflow

requirements:
  - class: SubworkflowFeatureRequirement

inputs:
  asn_cache: Directory
  seqids: File
  16s_blastdb_dir: Directory
  23s_blastdb_dir: Directory
  model_path: File
  rfam_amendments: File
  rfam_stockholm: File

outputs:
  # asncache:
  #   type: Directory
  #   outputSource: generate_16s/asncache
  annotations_5s:
    type: File
    outputSource: generate_5s/annots
  annotations_16s:
    type: File
    outputSource: generate_16s/annotations
  annotations_23s:
    type: File
    outputSource: generate_23s/annotations
    
steps:
  generate_5s:
    run: wf_gcmsearch.cwl
    in:
      asn_cache: asn_cache
      seqids: seqids
      model_path: model_path
      rfam_amendments: rfam_amendments
      rfam_stockholm: rfam_stockholm
    out: [ annots ]

  generate_16s:
    run: wf_blastn.cwl
    in:
      asn_cache: asn_cache
      seqids: seqids
      blastdb_dir: 16s_blastdb_dir
      blastdb:
        default: blastdb
      product_name: 
        default: "16S ribosomal RNA"
      outname:
        default: annotations_16s.asn
    out: [asncache, annotations]

  generate_23s:
    run: wf_blastn.cwl
    in:
      asn_cache: asn_cache
      seqids: seqids
      blastdb_dir: 23s_blastdb_dir
      blastdb:
        default: Ribosom23S
      product_name: 
        default: "23S ribosomal RNA"
      outname:
        default: annotations_23s.asn
    out: [asncache, annotations]