Merge branch 'master' into firstbranch

modules/bowtie2/align/main.nf

@@ -29,6 +29,8 @@ process BOWTIE2_ALIGN {
         def unaligned = save_unaligned ? "--un-gz ${prefix}.unmapped.fastq.gz" : ''
         """
         INDEX=`find -L ./ -name "*.rev.1.bt2" | sed 's/.rev.1.bt2//'`
+        [ -z "\$INDEX" ] && INDEX=`find -L ./ -name "*.rev.1.bt2l" | sed 's/.rev.1.bt2l//'`
+        [ -z "\$INDEX" ] && echo "BT2 index files not found" 1>&2 && exit 1
         bowtie2 \\
             -x \$INDEX \\
             -U $reads \\
@@ -49,6 +51,8 @@ process BOWTIE2_ALIGN {
         def unaligned = save_unaligned ? "--un-conc-gz ${prefix}.unmapped.fastq.gz" : ''
         """
         INDEX=`find -L ./ -name "*.rev.1.bt2" | sed 's/.rev.1.bt2//'`
+        [ -z "\$INDEX" ] && INDEX=`find -L ./ -name "*.rev.1.bt2l" | sed 's/.rev.1.bt2l//'`
+        [ -z "\$INDEX" ] && echo "BT2 index files not found" 1>&2 && exit 1
         bowtie2 \\
             -x \$INDEX \\
             -1 ${reads[0]} \\

modules/gatk4/mergebamalignment/main.nf

@@ -43,4 +43,15 @@ process GATK4_MERGEBAMALIGNMENT {
         gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
     END_VERSIONS
     """
+
+    stub:
+    def prefix = task.ext.prefix ?: "${meta.id}"
+    """
+    touch ${prefix}.bam
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
+    END_VERSIONS
+    """
 }

modules/gatk4/mutect2/main.nf

@@ -57,4 +57,18 @@ process GATK4_MUTECT2 {
         gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
     END_VERSIONS
     """
+
+    stub:
+    def prefix = task.ext.prefix ?: "${meta.id}"
+    """
+    touch ${prefix}.vcf.gz
+    touch ${prefix}.vcf.gz.tbi
+    touch ${prefix}.vcf.gz.stats
+    touch ${prefix}.f1r2.tar.gz
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
+    END_VERSIONS
+    """
 }

modules/gatk4/revertsam/main.nf

@@ -39,4 +39,15 @@ process GATK4_REVERTSAM {
         gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
     END_VERSIONS
     """
+
+    stub:
+    def prefix = task.ext.prefix ?: "${meta.id}"
+    """
+    touch ${prefix}.reverted.bam
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
+    END_VERSIONS
+    """
 }

modules/gatk4/samtofastq/main.nf

@@ -40,4 +40,17 @@ process GATK4_SAMTOFASTQ {
         gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
     END_VERSIONS
     """
+
+    stub:
+    def prefix = task.ext.prefix ?: "${meta.id}"
+    """
+    touch ${prefix}.fastq.gz
+    touch ${prefix}_1.fastq.gz
+    touch ${prefix}_2.fastq.gz
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        gatk4: \$(echo \$(gatk --version 2>&1) | sed 's/^.*(GATK) v//; s/ .*\$//')
+    END_VERSIONS
+    """
 }

modules/happy/happy/main.nf

@@ -0,0 +1,42 @@
+def VERSION = '0.3.14'
+
+process HAPPY_HAPPY {
+    tag "$meta.id"
+    label 'process_medium'
+
+    conda (params.enable_conda ? "bioconda::hap.py=0.3.14" : null)
+    container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ?
+        'https://depot.galaxyproject.org/singularity/hap.py:0.3.14--py27h5c5a3ab_0':
+        'quay.io/biocontainers/hap.py:0.3.14--py27h5c5a3ab_0' }"
+
+    input:
+    tuple val(meta), path(truth_vcf), path(query_vcf), path(bed)
+    tuple path(fasta), path(fasta_fai)
+
+    output:
+    tuple val(meta), path('*.csv'), path('*.json')  , emit: metrics
+    path "versions.yml"                             , emit: versions
+
+    when:
+    task.ext.when == null || task.ext.when
+
+    script:
+    def args = task.ext.args ?: ''
+    def prefix = task.ext.prefix ?: "${meta.id}"
+
+    """
+    hap.py \\
+        $truth_vcf \\
+        $query_vcf \\
+        $args \\
+        --reference $fasta \\
+        --threads $task.cpus \\
+        -R $bed \\
+        -o $prefix
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        hap.py: $VERSION
+    END_VERSIONS
+    """
+}

modules/happy/happy/meta.yml

@@ -0,0 +1,67 @@
+name: "happy_happy"
+description: Hap.py is a tool to compare diploid genotypes at haplotype level. Rather than comparing VCF records row by row, hap.py will generate and match alternate sequences in a superlocus. A superlocus is a small region of the genome (sized between 1 and around 1000 bp) that contains one or more variants.
+keywords:
+  - happy
+  - benchmark
+  - haplotype
+tools:
+  - "happy":
+      description: "Haplotype VCF comparison tools"
+      homepage: "https://www.illumina.com/products/by-type/informatics-products/basespace-sequence-hub/apps/hap-py-benchmarking.html"
+      documentation: "https://github.com/Illumina/hap.py"
+      tool_dev_url: "https://github.com/Illumina/hap.py"
+      doi: ""
+      licence: "['BSD-2-clause']"
+
+input:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - truth_vcf:
+      type: file
+      description: gold standard VCF file
+      pattern: "*.{vcf,vcf.gz}"
+  - query_vcf:
+      type: file
+      description: VCF/GVCF file to query
+      pattern: "*.{vcf,vcf.gz}"
+  - bed:
+      type: file
+      description: BED file
+      pattern: "*.bed"
+  - fasta:
+      type: file
+      description: FASTA file of the reference genome
+      pattern: "*.{fa,fasta}"
+  - fasta_fai:
+      type: file
+      description: The index of the reference FASTA
+      pattern: "*.fai"
+
+output:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - summary:
+      type: file
+      description: A CSV file containing the summary of the benchmarking
+      pattern: "*.summary.csv"
+  - extended:
+      type: file
+      description: A CSV file containing extended info of the benchmarking
+      pattern: "*.extended.csv"
+  - runinfo:
+      type: file
+      description: A JSON file containing the run info
+      pattern: "*.runinfo.json"
+  - versions:
+      type: file
+      description: File containing software versions
+      pattern: "versions.yml"
+
+authors:
+  - "@nvnieuwk"

modules/happy/prepy/main.nf

@@ -0,0 +1,41 @@
+def VERSION = '0.3.14'
+
+process HAPPY_PREPY {
+    tag "$meta.id"
+    label 'process_medium'
+
+    conda (params.enable_conda ? "bioconda::hap.py=0.3.14" : null)
+    container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ?
+        'https://depot.galaxyproject.org/singularity/hap.py:0.3.14--py27h5c5a3ab_0':
+        'quay.io/biocontainers/hap.py:0.3.14--py27h5c5a3ab_0' }"
+
+    input:
+    tuple val(meta), path(vcf), path(bed)
+    tuple path(fasta), path(fasta_fai)
+
+    output:
+    tuple val(meta), path('*.vcf.gz')  , emit: preprocessed_vcf
+    path "versions.yml"                , emit: versions
+
+    when:
+    task.ext.when == null || task.ext.when
+
+    script:
+    def args = task.ext.args ?: ''
+    def prefix = task.ext.prefix ?: "${meta.id}"
+
+    """
+    pre.py \\
+        $args \\
+        -R $bed \\
+        --reference $fasta \\
+        --threads $task.cpus \\
+        $vcf \\
+        ${prefix}.vcf.gz
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        pre.py: $VERSION
+    END_VERSIONS
+    """
+}

modules/happy/prepy/meta.yml

@@ -0,0 +1,55 @@
+name: "happy_prepy"
+description: Pre.py is a preprocessing tool made to preprocess VCF files for Hap.py
+keywords:
+  - happy
+  - benchmark
+  - haplotype
+tools:
+  - "happy":
+      description: "Haplotype VCF comparison tools"
+      homepage: "https://www.illumina.com/products/by-type/informatics-products/basespace-sequence-hub/apps/hap-py-benchmarking.html"
+      documentation: "https://github.com/Illumina/hap.py"
+      tool_dev_url: "https://github.com/Illumina/hap.py"
+      doi: ""
+      licence: "['BSD-2-clause']"
+
+input:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - vcf:
+      type: file
+      description: VCF file to preprocess
+      pattern: "*.{vcf,vcf.gz}"
+  - bed:
+      type: file
+      description: BED file
+      pattern: "*.bed"
+  - fasta:
+      type: file
+      description: FASTA file of the reference genome
+      pattern: "*.{fa,fasta}"
+  - fasta_fai:
+      type: file
+      description: The index of the reference FASTA
+      pattern: "*.fai"
+
+output:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - vcf:
+      type: file
+      description: A preprocessed VCF file
+      pattern: "*.vcf.gz"
+  - versions:
+      type: file
+      description: File containing software versions
+      pattern: "versions.yml"
+
+authors:
+  - "@nvnieuwk"

modules/metaphlan3/main.nf

@@ -23,7 +23,7 @@ process METAPHLAN3 {
     script:
     def args = task.ext.args ?: ''
     def prefix = task.ext.prefix ?: "${meta.id}"
-    def input_type  = ("$input".endsWith(".fastq.gz")) ? "--input_type fastq" :  ("$input".contains(".fasta")) ? "--input_type fasta" : ("$input".endsWith(".bowtie2out.txt")) ? "--input_type bowtie2out" : "--input_type sam"
+    def input_type  = ("$input".endsWith(".fastq.gz") || "$input".endsWith(".fq.gz")) ? "--input_type fastq" :  ("$input".contains(".fasta")) ? "--input_type fasta" : ("$input".endsWith(".bowtie2out.txt")) ? "--input_type bowtie2out" : "--input_type sam"
     def input_data  = ("$input_type".contains("fastq")) && !meta.single_end ? "${input[0]},${input[1]}" : "$input"
     def bowtie2_out = "$input_type" == "--input_type bowtie2out" || "$input_type" == "--input_type sam" ? '' : "--bowtie2out ${prefix}.bowtie2out.txt"
 

modules/samtools/view/main.nf

@@ -41,4 +41,16 @@ process SAMTOOLS_VIEW {
         samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
     END_VERSIONS
     """
+
+    stub:
+    def prefix = task.ext.prefix ?: "${meta.id}"
+    """
+    touch ${prefix}.bam
+    touch ${prefix}.cram
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        samtools: \$(echo \$(samtools --version 2>&1) | sed 's/^.*samtools //; s/Using.*\$//')
+    END_VERSIONS
+    """
 }

modules/shigatyper/main.nf

@@ -0,0 +1,64 @@
+process SHIGATYPER {
+    tag "$meta.id"
+    label 'process_low'
+
+    conda (params.enable_conda ? "bioconda::shigatyper=2.0.1" : null)
+    container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ?
+        'https://depot.galaxyproject.org/singularity/shigatyper%3A2.0.1--pyhdfd78af_0':
+        'quay.io/biocontainers/shigatyper:2.0.1--pyhdfd78af_0' }"
+
+    input:
+    tuple val(meta), path(reads)
+
+    output:
+    tuple val(meta), path("${prefix}.tsv")     , emit: tsv
+    tuple val(meta), path("${prefix}-hits.tsv"), optional: true, emit: hits
+    path "versions.yml"                        , emit: versions
+
+    when:
+    task.ext.when == null || task.ext.when
+
+    script:
+    def args = task.ext.args ?: ''
+    prefix = task.ext.prefix ?: "${meta.id}"
+
+    if (meta.is_ont) {
+        """
+        shigatyper \\
+            $args \\
+            --SE $reads \\
+            --ont \\
+            --name $prefix
+
+        cat <<-END_VERSIONS > versions.yml
+        "${task.process}":
+            shigatyper: \$(echo \$(shigatyper --version 2>&1) | sed 's/^.*ShigaTyper //' )
+        END_VERSIONS
+        """
+    } else if (meta.single_end) {
+        """
+        shigatyper \\
+            $args \\
+            --SE $reads \\
+            --name $prefix
+
+        cat <<-END_VERSIONS > versions.yml
+        "${task.process}":
+            shigatyper: \$(echo \$(shigatyper --version 2>&1) | sed 's/^.*ShigaTyper //' )
+        END_VERSIONS
+        """
+    } else {
+        """
+        shigatyper \\
+            $args \\
+            --R1 ${reads[0]} \\
+            --R2 ${reads[1]} \\
+            --name $prefix
+
+        cat <<-END_VERSIONS > versions.yml
+        "${task.process}":
+            shigatyper: \$(echo \$(shigatyper --version 2>&1) | sed 's/^.*ShigaTyper //' )
+        END_VERSIONS
+        """
+    }
+}

modules/shigatyper/meta.yml

@@ -0,0 +1,47 @@
+name: "shigatyper"
+description: Determine Shigella serotype from Illumina or Oxford Nanopore reads
+keywords:
+  - fastq
+  - shigella
+  - serotype
+tools:
+  - "shigatyper":
+      description: "Typing tool for Shigella spp. from WGS Illumina sequencing"
+      homepage: "https://github.com/CFSAN-Biostatistics/shigatyper"
+      documentation: "https://github.com/CFSAN-Biostatistics/shigatyper"
+      tool_dev_url: "https://github.com/CFSAN-Biostatistics/shigatyper"
+      doi: "10.1128/AEM.00165-19"
+      licence: "['Public Domain']"
+
+input:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false, is_ont:false ]
+  - reads:
+      type: file
+      description: Illumina or Nanopore FASTQ file
+      pattern: "*.fastq.gz"
+
+output:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - versions:
+      type: file
+      description: File containing software versions
+      pattern: "versions.yml"
+  - tsv:
+      type: file
+      description: A TSV formatted file with ShigaTyper results
+      pattern: "*.tsv"
+  - hits:
+      type: file
+      description: A TSV formatted file with individual gene hits
+      pattern: "*-hits.tsv"
+
+authors:
+  - "@rpetit3"

modules/slimfastq/main.nf

@@ -0,0 +1,52 @@
+def VERSION = '2.04'
+
+process SLIMFASTQ {
+    tag "$meta.id"
+    label 'process_low'
+
+    conda (params.enable_conda ? "bioconda::slimfastq=2.04" : null)
+    container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ?
+        'https://depot.galaxyproject.org/singularity/slimfastq:2.04--h87f3376_2':
+        'quay.io/biocontainers/slimfastq:2.04--h87f3376_2' }"
+
+    input:
+    tuple val(meta), path(fastq)
+
+    output:
+    tuple val(meta), path("*.sfq"), emit: sfq
+    path "versions.yml"           , emit: versions
+
+    when:
+    task.ext.when == null || task.ext.when
+
+    script:
+    def args = task.ext.args ?: ''
+    def prefix = task.ext.prefix ?: "${meta.id}"
+    if (meta.single_end) {
+        """
+        gzip -d -c '${fastq}' | slimfastq \\
+            $args \\
+            -f '${prefix}.sfq'
+
+        cat <<-END_VERSIONS > versions.yml
+        "${task.process}":
+            slimfastq: ${VERSION}
+        END_VERSIONS
+        """
+    } else {
+        """
+        gzip -d -c '${fastq[0]}' | slimfastq \\
+            $args \\
+            -f '${prefix}_1.sfq'
+
+        gzip -d -c '${fastq[1]}' | slimfastq \\
+            $args \\
+            -f '${prefix}_2.sfq'
+
+        cat <<-END_VERSIONS > versions.yml
+        "${task.process}":
+            slimfastq: ${VERSION}
+        END_VERSIONS
+        """
+    }
+}

modules/slimfastq/meta.yml

@@ -0,0 +1,41 @@
+name: "slimfastq"
+description: Fast, efficient, lossless compression of FASTQ files.
+keywords:
+  - FASTQ
+  - compression
+  - lossless
+tools:
+  - "slimfastq":
+      description: "slimfastq efficiently compresses/decompresses FASTQ files"
+      homepage: "https://github.com/Infinidat/slimfastq"
+      tool_dev_url: "https://github.com/Infinidat/slimfastq"
+      licence: "['BSD-3-clause']"
+
+input:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - fastq:
+      type: file
+      description: Either a single-end FASTQ file or paired-end files.
+      pattern: "*.{fq.gz,fastq.gz}"
+
+output:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - versions:
+      type: file
+      description: File containing software versions
+      pattern: "versions.yml"
+  - sfq:
+      type: file
+      description: Either one or two sequence files in slimfastq compressed format.
+      pattern: "*.{sfq}"
+
+authors:
+  - "@Midnighter"

modules/srst2/srst2/main.nf

@@ -0,0 +1,47 @@
+process SRST2_SRST2 {
+    tag "${meta.id}"
+    label 'process_low'
+
+    conda (params.enable_conda ? "bioconda::srst2=0.2.0" : null)
+    container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ?
+        'https://depot.galaxyproject.org/singularity/srst2%3A0.2.0--py27_2':
+        'quay.io/biocontainers/srst2:0.2.0--py27_2'}"
+
+    input:
+    tuple val(meta), path(fastq_s), path(db)
+
+    output:
+    tuple val(meta), path("*_genes_*_results.txt")               , optional:true, emit: gene_results
+    tuple val(meta), path("*_fullgenes_*_results.txt")           , optional:true, emit: fullgene_results
+    tuple val(meta), path("*_mlst_*_results.txt")                , optional:true, emit: mlst_results
+    tuple val(meta), path("*.pileup")                            ,                emit: pileup
+    tuple val(meta), path("*.sorted.bam")                        ,                emit: sorted_bam
+    path "versions.yml"                                          ,                emit: versions
+
+    when:
+    task.ext.when == null || task.ext.when
+
+    script:
+    def args = task.ext.args ?: ""
+    def prefix = task.ext.prefix ?: "${meta.id}"
+    def read_s = meta.single_end ? "--input_se ${fastq_s}" : "--input_pe ${fastq_s[0]} ${fastq_s[1]}"
+    if (meta.db=="gene") {
+        database = "--gene_db ${db}"
+    } else if (meta.db=="mlst") {
+        database = "--mlst_db ${db}"
+    } else {
+        error "Please set meta.db to either \"gene\" or \"mlst\""
+    }
+    """
+    srst2 \\
+        ${read_s} \\
+        --threads $task.cpus \\
+        --output ${prefix} \\
+        ${database} \\
+        $args
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        srst2: \$(echo \$(srst2 --version 2>&1) | sed 's/srst2 //' ))
+    END_VERSIONS
+    """
+}

modules/srst2/srst2/meta.yml

@@ -0,0 +1,72 @@
+name: srst2_srst2
+description: |
+  Short Read Sequence Typing for Bacterial Pathogens is a program designed to take Illumina sequence data,
+  a MLST database and/or a database of gene sequences (e.g. resistance genes, virulence genes, etc)
+  and report the presence of STs and/or reference genes.
+keywords:
+  - mlst
+  - typing
+  - illumina
+tools:
+  - srst2:
+      description: "Short Read Sequence Typing for Bacterial Pathogens"
+      homepage: "http://katholt.github.io/srst2/"
+      documentation: "https://github.com/katholt/srst2/blob/master/README.md"
+      tool_dev_url: "https://github.com/katholt/srst2"
+      doi: "10.1186/s13073-014-0090-6"
+      licence: ["BSD"]
+
+input:
+  - meta:
+      type: map0.2.0-4
+      description: |
+        Groovy Map containing sample information
+        id: should be the identification number or sample name
+        single_end: should be true for single end data and false for paired in data
+        db: should be either 'gene' to use the --gene_db option or "mlst" to use the --mlst_db option
+        e.g. [ id:'sample', single_end:false , db:'gene']
+  - fasta:
+      type: file
+      description: |
+        gzipped fasta file. If files are NOT in
+        MiSeq format sample_S1_L001_R1_001.fastq.gz uses --forward and --reverse parameters; otherwise
+        default is _1, i.e. expect forward reads as sample_1.fastq.gz).
+      pattern: "*.fastq.gz"
+  - db:
+      type: file
+      description: Database in FASTA format
+      pattern: "*.fasta"
+
+output:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'sample', single_end:false ]
+  - versions:
+      type: file
+      description: File containing software versions
+      pattern: "versions.yml"
+  - txt:
+      type: file
+      description: A detailed report, with one row per gene per sample described here github.com/katholt/srst2#gene-typing
+      pattern: "*_fullgenes_*_results.txt"
+  - txt:
+      type: file
+      description: A tabulated summary report of samples x genes.
+      pattern: "*_genes_*_results.txt"
+  - txt:
+      type: file
+      description: A tabulated summary report of mlst subtyping.
+      pattern: "*_mlst_*_results.txt"
+  - bam:
+      type: file
+      description: Sorted BAM file
+      pattern: "*.sorted.bam"
+  - pileup:
+      type: file
+      description: SAMtools pileup file
+      pattern: "*.pileup"
+
+authors:
+  - "@jvhagey"

modules/vardictjava/main.nf

@@ -0,0 +1,49 @@
+def VERSION = '1.8.3'
+
+process VARDICTJAVA {
+    tag "$meta.id"
+    label 'process_medium'
+
+    conda (params.enable_conda ? "bioconda::vardict-java=1.8.3" : null)
+    container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ?
+        'https://depot.galaxyproject.org/singularity/vardict-java:1.8.3--hdfd78af_0':
+        'quay.io/biocontainers/vardict-java:1.8.3--hdfd78af_0' }"
+
+    input:
+    tuple val(meta), path(bam), path(bai), path(bed)
+    tuple path(fasta), path(fasta_fai)
+
+    output:
+    tuple val(meta), path("*.vcf.gz"), emit: vcf
+    path "versions.yml"              , emit: versions
+
+    when:
+    task.ext.when == null || task.ext.when
+
+    script:
+    def args = task.ext.args ?: ''
+    def args2 = task.ext.args2 ?: ''
+    def prefix = task.ext.prefix ?: "${meta.id}"
+
+    """
+    vardict-java \\
+        $args \\
+        -c 1 -S 2 -E 3 \\
+        -b $bam \\
+        -th $task.cpus \\
+        -N $prefix \\
+        -G $fasta \\
+        $bed \\
+        | teststrandbias.R \\
+        | var2vcf_valid.pl \\
+            $args2 \\
+            -N $prefix \\
+        | gzip -c > ${prefix}.vcf.gz
+
+    cat <<-END_VERSIONS > versions.yml
+    "${task.process}":
+        vardict-java: $VERSION
+        var2vcf_valid.pl: \$(echo \$(var2vcf_valid.pl -h | sed -n 2p | awk '{ print \$2 }'))
+    END_VERSIONS
+    """
+}

modules/vardictjava/meta.yml

@@ -0,0 +1,60 @@
+name: "vardictjava"
+
+description: The Java port of the VarDict variant caller
+keywords:
+  - variant calling
+  - VarDict
+  - AstraZeneca
+tools:
+  - "vardictjava":
+      description: "Java port of the VarDict variant discovery program"
+      homepage: "https://github.com/AstraZeneca-NGS/VarDictJava"
+      documentation: "https://github.com/AstraZeneca-NGS/VarDictJava"
+      tool_dev_url: "https://github.com/AstraZeneca-NGS/VarDictJava"
+      doi: "10.1093/nar/gkw227 "
+      licence: "['MIT']"
+
+input:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - bam:
+      type: file
+      description: BAM/SAM file
+      pattern: "*.{bam,sam}"
+  - bai:
+      type: file
+      description: Index of the BAM file
+      pattern: "*.bai"
+  - fasta:
+      type: file
+      description: FASTA of the reference genome
+      pattern: "*.{fa,fasta}"
+  - fasta_fai:
+      type: file
+      description: The index of the FASTA of the reference genome
+      pattern: "*.fai"
+  - bed:
+      type: file
+      description: BED with the regions of interest
+      pattern: "*.bed"
+
+output:
+  - meta:
+      type: map
+      description: |
+        Groovy Map containing sample information
+        e.g. [ id:'test', single_end:false ]
+  - versions:
+      type: file
+      description: File containing software versions
+      pattern: "versions.yml"
+  - vcf:
+      type: file
+      description: VCF file output
+      pattern: "*.vcf.gz"
+
+authors:
+  - "@nvnieuwk"

tests/config/pytest_modules.yml

@@ -891,6 +891,14 @@ hamronization/summarize:
   - modules/hamronization/summarize/**
   - tests/modules/hamronization/summarize/**
 
+happy/happy:
+  - modules/happy/happy/**
+  - tests/modules/happy/happy/**
+
+happy/prepy:
+  - modules/happy/prepy/**
+  - tests/modules/happy/prepy/**
+
 hicap:
   - modules/hicap/**
   - tests/modules/hicap/**
@@ -1719,6 +1727,10 @@ seqwish/induce:
   - modules/seqwish/induce/**
   - tests/modules/seqwish/induce/**
 
+shigatyper:
+  - modules/shigatyper/**
+  - tests/modules/shigatyper/**
+
 shovill:
   - modules/shovill/**
   - tests/modules/shovill/**
@@ -1727,6 +1739,10 @@ sistr:
   - modules/sistr/**
   - tests/modules/sistr/**
 
+slimfastq:
+  - modules/slimfastq/**
+  - tests/modules/slimfastq/**
+
 snapaligner/index:
   - modules/snapaligner/index/**
   - tests/modules/snapaligner/index/**
@@ -1775,6 +1791,10 @@ sratools/prefetch:
   - modules/sratools/prefetch/**
   - tests/modules/sratools/prefetch/**
 
+srst2/srst2:
+  - modules/srst2/srst2/**
+  - tests/modules/srst2/srst2/**
+
 ssuissero:
   - modules/ssuissero/**
   - tests/modules/ssuissero/**
@@ -1916,6 +1936,10 @@ unzip:
   - modules/unzip/**
   - tests/modules/unzip/**
 
+vardictjava:
+  - modules/vardictjava/**
+  - tests/modules/vardictjava/**
+
 variantbam:
   - modules/variantbam/**
   - tests/modules/variantbam/**

tests/modules/bowtie2/align/main.nf

@@ -32,4 +32,4 @@ workflow test_bowtie2_align_paired_end {
 
     BOWTIE2_BUILD ( fasta )
     BOWTIE2_ALIGN ( input, BOWTIE2_BUILD.out.index, save_unaligned )
-}
+}

tests/modules/bowtie2/align/nextflow.config

@@ -1,5 +1,16 @@
+params {
+    force_large_index = false
+}
+
 process {
 
     publishDir = { "${params.outdir}/${task.process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()}" }
-
+}
+
+if (params.force_large_index) {
+    process {
+        withName: BOWTIE2_BUILD {
+            ext.args = '--large-index'
+        }
+    }
 }

tests/modules/bowtie2/align/test.yml

@@ -39,3 +39,45 @@
       md5sum: 52be6950579598a990570fbcf5372184
     - path: ./output/bowtie2/bowtie2/genome.rev.2.bt2
       md5sum: e3b4ef343dea4dd571642010a7d09597
+
+- name: bowtie2 align single-end large-index
+  command: nextflow run ./tests/modules/bowtie2/align -entry test_bowtie2_align_single_end -c ./tests/config/nextflow.config -c ./tests/modules/bowtie2/align/nextflow.config --force_large_index
+  tags:
+    - bowtie2
+    - bowtie2/align
+  files:
+    - path: ./output/bowtie2/test.bam
+    - path: ./output/bowtie2/test.bowtie2.log
+    - path: ./output/bowtie2/bowtie2/genome.3.bt2l
+      md5sum: 8952b3e0b1ce9a7a5916f2e147180853
+    - path: ./output/bowtie2/bowtie2/genome.2.bt2l
+      md5sum: 22c284084784a0720989595e0c9461fd
+    - path: ./output/bowtie2/bowtie2/genome.1.bt2l
+      md5sum: 07d811cd4e350d56267183d2ac7023a5
+    - path: ./output/bowtie2/bowtie2/genome.4.bt2l
+      md5sum: c25be5f8b0378abf7a58c8a880b87626
+    - path: ./output/bowtie2/bowtie2/genome.rev.1.bt2l
+      md5sum: fda48e35925fb24d1c0785f021981e25
+    - path: ./output/bowtie2/bowtie2/genome.rev.2.bt2l
+      md5sum: 802c26d32b970e1b105032b7ce7348b4
+
+- name: bowtie2 align paired-end large-index
+  command: nextflow run ./tests/modules/bowtie2/align -entry test_bowtie2_align_paired_end -c ./tests/config/nextflow.config -c ./tests/modules/bowtie2/align/nextflow.config --force_large_index
+  tags:
+    - bowtie2
+    - bowtie2/align
+  files:
+    - path: ./output/bowtie2/test.bam
+    - path: ./output/bowtie2/test.bowtie2.log
+    - path: ./output/bowtie2/bowtie2/genome.3.bt2l
+      md5sum: 8952b3e0b1ce9a7a5916f2e147180853
+    - path: ./output/bowtie2/bowtie2/genome.2.bt2l
+      md5sum: 22c284084784a0720989595e0c9461fd
+    - path: ./output/bowtie2/bowtie2/genome.1.bt2l
+      md5sum: 07d811cd4e350d56267183d2ac7023a5
+    - path: ./output/bowtie2/bowtie2/genome.4.bt2l
+      md5sum: c25be5f8b0378abf7a58c8a880b87626
+    - path: ./output/bowtie2/bowtie2/genome.rev.1.bt2l
+      md5sum: fda48e35925fb24d1c0785f021981e25
+    - path: ./output/bowtie2/bowtie2/genome.rev.2.bt2l
+      md5sum: 802c26d32b970e1b105032b7ce7348b4

tests/modules/gatk4/mergebamalignment/main.nf

@@ -14,3 +14,14 @@ workflow test_gatk4_mergebamalignment {
 
     GATK4_MERGEBAMALIGNMENT ( input, fasta, dict )
 }
+
+workflow test_gatk4_mergebamalignment_stubs {
+     input    = [ [ id:'test' ], // meta map
+                 "test_foo.bam",
+                 "test_bar.bam"
+               ]
+    fasta    = "genome.fasta"
+    dict     = "genome.fasta.dict"
+
+    GATK4_MERGEBAMALIGNMENT ( input, fasta, dict )
+}

tests/modules/gatk4/mergebamalignment/test.yml

@@ -7,3 +7,12 @@
     - path: output/gatk4/test.bam
       md5sum: e6f1b343700b7ccb94e81ae127433988
     - path: output/gatk4/versions.yml
+
+- name: gatk4 mergebamalignment test_gatk4_mergebamalignment_stubs
+  command: nextflow run ./tests/modules/gatk4/mergebamalignment -entry test_gatk4_mergebamalignment -c ./tests/config/nextflow.config -c ./tests/modules/gatk4/mergebamalignment/nextflow.config -stub-run
+  tags:
+    - gatk4
+    - gatk4/mergebamalignment
+  files:
+    - path: output/gatk4/test.bam
+    - path: output/gatk4/versions.yml

tests/modules/gatk4/mutect2/main.nf

@@ -118,3 +118,25 @@ workflow test_gatk4_mutect2_mitochondria {
 
     GATK4_MUTECT2_MITO ( input, fasta, fai, dict, [], [], [], [] )
 }
+
+workflow test_gatk4_mutect2_tumor_normal_pair_f1r2_stubs {
+    input = [ [ id:'test', normal_id:'normal', tumor_id:'tumour' ], // meta map
+              [ "foo_paired.bam",
+                "foo_paired2.bam"
+                ],
+              [ "foo_paired.bam.bai",
+                "foo_paired2.bam.bai"
+                ],
+              []
+            ]
+
+    fasta = "genome.fasta"
+    fai = "genome.fasta.fai"
+    dict = "genome.fasta.dict"
+    germline_resource = "genome_gnomAD.r2.1.1.vcf.gz"
+    germline_resource_tbi = "genome_gnomAD.r2.1.1.vcf.gz.tbi"
+    panel_of_normals = "genome_mills_and_1000G.indels.hg38.vcf.gz"
+    panel_of_normals_tbi = "genome_mills_and_1000G.indels.hg38.vcf.gz.tbi"
+
+    GATK4_MUTECT2_F1R2 ( input, fasta, fai, dict, germline_resource, germline_resource_tbi, panel_of_normals, panel_of_normals_tbi )
+}

tests/modules/gatk4/mutect2/test.yml

@@ -69,3 +69,15 @@
       md5sum: fc6ea14ca2da346babe78161beea28c9
     - path: output/gatk4/test.vcf.gz.tbi
     - path: output/gatk4/versions.yml
+
+- name: gatk4 mutect2 test_gatk4_mutect2_tumor_normal_pair_f1r2_stubs
+  command: nextflow run ./tests/modules/gatk4/mutect2 -entry test_gatk4_mutect2_tumor_normal_pair_f1r2 -c ./tests/config/nextflow.config -c ./tests/modules/gatk4/mutect2/nextflow.config -stub-run
+  tags:
+    - gatk4
+    - gatk4/mutect2
+  files:
+    - path: output/gatk4/test.f1r2.tar.gz
+    - path: output/gatk4/test.vcf.gz
+    - path: output/gatk4/test.vcf.gz.stats
+    - path: output/gatk4/test.vcf.gz.tbi
+    - path: output/gatk4/versions.yml

tests/modules/gatk4/revertsam/main.nf

@@ -11,3 +11,11 @@ workflow test_gatk4_revertsam {
 
     GATK4_REVERTSAM ( input )
 }
+
+workflow test_gatk4_revertsam_stubs {
+    input = [ [ id:'test' ], // meta map
+              "foo_paired_end.bam"
+            ]
+
+    GATK4_REVERTSAM ( input )
+}

tests/modules/gatk4/revertsam/test.yml

@@ -7,3 +7,12 @@
     - path: output/gatk4/test.reverted.bam
       md5sum: f783a88deb45c3a2c20ca12cbe1c5652
     - path: output/gatk4/versions.yml
+
+- name: gatk4 revertsam test_gatk4_revertsam_stubs
+  command: nextflow run ./tests/modules/gatk4/revertsam -entry test_gatk4_revertsam -c ./tests/config/nextflow.config -c ./tests/modules/gatk4/revertsam/nextflow.config -stub-run
+  tags:
+    - gatk4
+    - gatk4/revertsam
+  files:
+    - path: output/gatk4/test.reverted.bam
+    - path: output/gatk4/versions.yml

tests/modules/gatk4/samtofastq/main.nf

@@ -19,3 +19,11 @@ workflow test_gatk4_samtofastq_paired_end {
 
     GATK4_SAMTOFASTQ ( input )
 }
+
+workflow test_gatk4_samtofastq_paired_end_stubs {
+    input = [ [ id:'test', single_end: false ], // meta map
+              [ "foo_paired_end.bam" ]
+            ]
+
+    GATK4_SAMTOFASTQ ( input )
+}

tests/modules/gatk4/samtofastq/test.yml

@@ -19,3 +19,13 @@
     - path: output/gatk4/test_2.fastq.gz
       md5sum: 613bf64c023609e1c62ad6ce9e4be8d7
     - path: output/gatk4/versions.yml
+
+- name: gatk4 samtofastq test_gatk4_samtofastq_paired_end_stubs
+  command: nextflow run ./tests/modules/gatk4/samtofastq -entry test_gatk4_samtofastq_paired_end -c ./tests/config/nextflow.config -c ./tests/modules/gatk4/samtofastq/nextflow.config -stub-run
+  tags:
+    - gatk4
+    - gatk4/samtofastq
+  files:
+    - path: output/gatk4/test_1.fastq.gz
+    - path: output/gatk4/test_2.fastq.gz
+    - path: output/gatk4/versions.yml

tests/modules/happy/happy/main.nf

@@ -0,0 +1,39 @@
+#!/usr/bin/env nextflow
+
+nextflow.enable.dsl = 2
+
+include { HAPPY_HAPPY } from '../../../../modules/happy/happy/main.nf'
+
+workflow test_happy_vcf {
+    
+    input = [
+        [ id:'test' ], // meta map
+        file(params.test_data['homo_sapiens']['illumina']['test_rnaseq_vcf'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['illumina']['test_genome21_indels_vcf_gz'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_bed'], checkIfExists: true)        
+    ]
+
+    fasta = Channel.value([
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta_fai'], checkIfExists: true)
+    ])
+
+    HAPPY_HAPPY ( input, fasta )
+}
+
+workflow test_happy_gvcf {
+    
+    input = [
+        [ id:'test' ], // meta map
+        file(params.test_data['homo_sapiens']['illumina']['test_rnaseq_vcf'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['illumina']['test_genome_vcf'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_bed'], checkIfExists: true)        
+    ]
+
+    fasta = Channel.value([
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta_fai'], checkIfExists: true)
+    ])
+
+    HAPPY_HAPPY ( input, fasta )
+}

tests/modules/happy/happy/nextflow.config

@@ -0,0 +1,5 @@
+process {
+
+    publishDir = { "${params.outdir}/${task.process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()}" }
+    
+}

tests/modules/happy/happy/test.yml

@@ -0,0 +1,27 @@
+- name: happy happy test_happy_vcf
+  command: nextflow run tests/modules/happy/happy -entry test_happy_vcf -c tests/config/nextflow.config
+  tags:
+    - happy
+    - happy/happy
+  files:
+    - path: output/happy/test.extended.csv
+      md5sum: ef79c7c789ef4f146ca2e50dafaf22b3
+    - path: output/happy/test.runinfo.json
+    - path: output/happy/test.summary.csv
+      md5sum: f8aa5d36d3c48dede2f607fd565894ad
+    - path: output/happy/versions.yml
+      md5sum: 82243bf6dbdc71aa63211ee2a89f47f2
+
+- name: happy happy test_happy_gvcf
+  command: nextflow run tests/modules/happy/happy -entry test_happy_gvcf -c tests/config/nextflow.config
+  tags:
+    - happy
+    - happy/happy
+  files:
+    - path: output/happy/test.extended.csv
+      md5sum: 3d5c21b67a259a3f6dcb088d55b86cd3
+    - path: output/happy/test.runinfo.json
+    - path: output/happy/test.summary.csv
+      md5sum: 03044e9bb5a0c6f0947b7e910fc8a558
+    - path: output/happy/versions.yml
+      md5sum: 551fa216952d6f5de78e6e453b92aaab

tests/modules/happy/prepy/main.nf

@@ -0,0 +1,37 @@
+#!/usr/bin/env nextflow
+
+nextflow.enable.dsl = 2
+
+include { HAPPY_PREPY } from '../../../../modules/happy/prepy/main.nf'
+
+workflow test_happy_prepy_vcf {
+    
+    input = [
+        [ id:'test' ], // meta map
+        file(params.test_data['homo_sapiens']['illumina']['test_genome21_indels_vcf_gz'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_bed'], checkIfExists: true)        
+    ]
+
+    fasta = Channel.value([
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta_fai'], checkIfExists: true)
+    ])
+
+    HAPPY_PREPY ( input, fasta )
+}
+
+workflow test_happy_prepy_gvcf {
+    
+    input = [
+        [ id:'test' ], // meta map
+        file(params.test_data['homo_sapiens']['illumina']['test_genome_vcf'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_bed'], checkIfExists: true)        
+    ]
+
+    fasta = Channel.value([
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta_fai'], checkIfExists: true)
+    ])
+
+    HAPPY_PREPY ( input, fasta )
+}

tests/modules/happy/prepy/nextflow.config

@@ -0,0 +1,5 @@
+process {
+
+    publishDir = { "${params.outdir}/${task.process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()}" }
+    
+}

tests/modules/happy/prepy/test.yml

@@ -0,0 +1,19 @@
+- name: happy prepy test_happy_prepy_vcf
+  command: nextflow run tests/modules/happy/prepy -entry test_happy_prepy_vcf -c tests/config/nextflow.config
+  tags:
+    - happy/prepy
+    - happy
+  files:
+    - path: output/happy/test.vcf.gz
+    - path: output/happy/versions.yml
+      md5sum: 814d20f1f29f23a3d21012748a5d6393
+
+- name: happy prepy test_happy_prepy_gvcf
+  command: nextflow run tests/modules/happy/prepy -entry test_happy_prepy_gvcf -c tests/config/nextflow.config
+  tags:
+    - happy/prepy
+    - happy
+  files:
+    - path: output/happy/test.vcf.gz
+    - path: output/happy/versions.yml
+      md5sum: 970a54de46e68ef6d5228a26eaa4c8e7

tests/modules/samtools/view/main.nf

@@ -22,3 +22,12 @@ workflow test_samtools_view_cram {
 
     SAMTOOLS_VIEW ( input, fasta )
 }
+
+workflow test_samtools_view_stubs {
+    input = [ [ id:'test', single_end:false ], // meta map
+                "foo_paired_end.bam",
+                []
+            ]
+
+    SAMTOOLS_VIEW ( input, [] )
+}

tests/modules/samtools/view/test.yml

@@ -14,3 +14,11 @@
     - samtools
   files:
     - path: output/samtools/test.cram
+
+- name: samtools view test_samtools_view_stubs
+  command: nextflow run ./tests/modules/samtools/view -entry test_samtools_view -c ./tests/config/nextflow.config -c ./tests/modules/samtools/view/nextflow.config -stub-run
+  tags:
+    - samtools/view
+    - samtools
+  files:
+    - path: output/samtools/test.bam

tests/modules/shigatyper/main.nf

@@ -0,0 +1,36 @@
+#!/usr/bin/env nextflow
+
+nextflow.enable.dsl = 2
+
+include { SHIGATYPER } from '../../../modules/shigatyper/main.nf'
+
+workflow test_shigatyper_pe {
+    
+    input = [ 
+        [ id:'test', single_end:false, is_ont:false ], // meta map
+        [ file(params.test_data['sarscov2']['illumina']['test_1_fastq_gz'], checkIfExists: true),
+          file(params.test_data['sarscov2']['illumina']['test_2_fastq_gz'], checkIfExists: true) ]
+    ]
+
+    SHIGATYPER ( input )
+}
+
+workflow test_shigatyper_se {
+    
+    input = [ 
+        [ id:'test', single_end:true, is_ont:false ], // meta map
+        [ file(params.test_data['sarscov2']['illumina']['test_1_fastq_gz'], checkIfExists: true) ]
+    ]
+
+    SHIGATYPER ( input )
+}
+
+workflow test_shigatyper_ont {
+    
+    input = [ 
+        [ id:'test', single_end:true, is_ont:true ], // meta map
+        [ file(params.test_data['sarscov2']['nanopore']['test_fastq_gz'], checkIfExists: true) ]
+    ]
+
+    SHIGATYPER ( input )
+}

tests/modules/shigatyper/nextflow.config

@@ -0,0 +1,5 @@
+process {
+
+    publishDir = { "${params.outdir}/${task.process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()}" }
+    
+}

tests/modules/shigatyper/test.yml

@@ -0,0 +1,29 @@
+- name: shigatyper test_shigatyper_pe
+  command: nextflow run tests/modules/shigatyper -entry test_shigatyper_pe -c tests/config/nextflow.config -c tests/modules/shigatyper/nextflow.config
+  tags:
+    - shigatyper
+  files:
+    - path: output/shigatyper/test.tsv
+      md5sum: 4f7d38c956993800546b9acb9881d717
+    - path: output/shigatyper/versions.yml
+      md5sum: d8ca45ed88dfba9bc570c01e4b49773b
+
+- name: shigatyper test_shigatyper_se
+  command: nextflow run tests/modules/shigatyper -entry test_shigatyper_se -c tests/config/nextflow.config -c tests/modules/shigatyper/nextflow.config
+  tags:
+    - shigatyper
+  files:
+    - path: output/shigatyper/test.tsv
+      md5sum: 4f7d38c956993800546b9acb9881d717
+    - path: output/shigatyper/versions.yml
+      md5sum: 8bbf165da5a5df3b7771a33aad197eec
+
+- name: shigatyper test_shigatyper_ont
+  command: nextflow run tests/modules/shigatyper -entry test_shigatyper_ont -c tests/config/nextflow.config -c tests/modules/shigatyper/nextflow.config
+  tags:
+    - shigatyper
+  files:
+    - path: output/shigatyper/test.tsv
+      md5sum: 4f7d38c956993800546b9acb9881d717
+    - path: output/shigatyper/versions.yml
+      md5sum: 0da333e1178e9e7e84a9116ad5a5ff71

tests/modules/slimfastq/main.nf

@@ -0,0 +1,46 @@
+#!/usr/bin/env nextflow
+
+nextflow.enable.dsl = 2
+
+include { SLIMFASTQ } from '../../../modules/slimfastq/main.nf'
+
+workflow test_slimfastq_single_end {
+
+    input = [
+        [ id:'test', single_end:true ], // meta map
+        file(params.test_data['sarscov2']['illumina']['test_1_fastq_gz'], checkIfExists: true)
+    ]
+
+    SLIMFASTQ ( input )
+}
+
+workflow test_slimfastq_paired_end {
+
+    input = [
+        [ id:'test', single_end:false ], // meta map
+        [ file(params.test_data['sarscov2']['illumina']['test_1_fastq_gz'], checkIfExists: true),
+         file(params.test_data['sarscov2']['illumina']['test_2_fastq_gz'], checkIfExists: true)]
+    ]
+
+    SLIMFASTQ ( input )
+}
+
+workflow test_slimfastq_nanopore {
+
+    input = [
+        [ id:'test', single_end:true ], // meta map
+        file(params.test_data['sarscov2']['nanopore']['test_fastq_gz'], checkIfExists: true)
+    ]
+
+    SLIMFASTQ ( input )
+}
+
+workflow test_slimfastq_pacbio {
+
+    input = [
+        [ id:'test', single_end:true ], // meta map
+        file(params.test_data['homo_sapiens']['pacbio']['ccs_fq_gz'], checkIfExists: true)
+    ]
+
+    SLIMFASTQ ( input )
+}

tests/modules/slimfastq/nextflow.config

@@ -0,0 +1,5 @@
+process {
+
+    publishDir = { "${params.outdir}/${task.process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()}" }
+    
+}

tests/modules/slimfastq/test.yml

@@ -0,0 +1,41 @@
+- name: slimfastq test_slimfastq_single_end
+  command: nextflow run tests/modules/slimfastq -entry test_slimfastq_single_end -c tests/config/nextflow.config
+  tags:
+    - slimfastq
+  files:
+    - path: output/slimfastq/test.sfq
+      md5sum: 6a942eeca6c99ee9a9a0cedab5d246f1
+    - path: output/slimfastq/versions.yml
+      md5sum: f52351f5c9e6259af02745c8eae5c780
+
+- name: slimfastq test_slimfastq_paired_end
+  command: nextflow run tests/modules/slimfastq -entry test_slimfastq_paired_end -c tests/config/nextflow.config
+  tags:
+    - slimfastq
+  files:
+    - path: output/slimfastq/test_1.sfq
+      md5sum: 6a942eeca6c99ee9a9a0cedab5d246f1
+    - path: output/slimfastq/test_2.sfq
+      md5sum: 0d2c60b52a39f7c2cb7843e848d90afd
+    - path: output/slimfastq/versions.yml
+      md5sum: 6239853705877651a4851c4cb6d62da4
+
+- name: slimfastq test_slimfastq_nanopore
+  command: nextflow run tests/modules/slimfastq -entry test_slimfastq_nanopore -c tests/config/nextflow.config
+  tags:
+    - slimfastq
+  files:
+    - path: output/slimfastq/test.sfq
+      md5sum: e17f14d64d3a75356b03ff2f9e8881f7
+    - path: output/slimfastq/versions.yml
+      md5sum: 33153f1103482a2bd35cb2f4c337c5e8
+
+- name: slimfastq test_slimfastq_pacbio
+  command: nextflow run tests/modules/slimfastq -entry test_slimfastq_pacbio -c tests/config/nextflow.config
+  tags:
+    - slimfastq
+  files:
+    - path: output/slimfastq/test.sfq
+      md5sum: 9e8389e47e6ddf8c25e92412dd628339
+    - path: output/slimfastq/versions.yml
+      md5sum: 1982789c3d5c7de37c0a9351e4ae63f7

tests/modules/srst2/srst2/main.nf

@@ -0,0 +1,53 @@
+#!/usr/bin/env nextflow
+
+nextflow.enable.dsl = 2
+
+include { SRST2_SRST2 } from '../../../../modules/srst2/srst2/main.nf'
+
+workflow test_srst2_srst2_exit {
+
+    input = [
+        [ id:'test', single_end:false, db:"test"], // meta map
+        [ file(params.test_data['bacteroides_fragilis']['illumina']['test1_1_fastq_gz'], checkIfExists: true),
+          file(params.test_data['bacteroides_fragilis']['illumina']['test1_2_fastq_gz'], checkIfExists: true) ],
+//        [("")]
+        file('https://raw.githubusercontent.com/nf-core/test-datasets/modules/data/delete_me/srst2/resFinder_20180221_srst2.fasta')
+    ]
+
+    SRST2_SRST2(input)
+}
+
+workflow test_srst2_srst2_mlst {
+
+    input = [
+        [ id:'test', single_end:false, db:"mlst"], // meta map
+        [ file("https://raw.githubusercontent.com/nf-core/test-datasets/modules/data/delete_me/srst2/SRR9067271_1.fastq.gz", checkIfExists: true),
+          file("https://raw.githubusercontent.com/nf-core/test-datasets/modules/data/delete_me/srst2/SRR9067271_2.fastq.gz", checkIfExists: true) ],
+        file('https://raw.githubusercontent.com/nf-core/test-datasets/modules/data/delete_me/srst2/MLST_DB.fas')
+    ]
+
+    SRST2_SRST2(input)
+}
+
+workflow test_srst2_srst2_paired_end {
+
+    input = [
+        [ id:'test', single_end:false, db:"gene"], // meta map
+        [ file(params.test_data['bacteroides_fragilis']['illumina']['test1_1_fastq_gz'], checkIfExists: true),
+          file(params.test_data['bacteroides_fragilis']['illumina']['test1_2_fastq_gz'], checkIfExists: true) ],
+        file('https://raw.githubusercontent.com/nf-core/test-datasets/modules/data/delete_me/srst2/resFinder_20180221_srst2.fasta')  // Change to params.test_data syntax after the data is included in tests/config/test_data.config
+    ]
+
+    SRST2_SRST2(input)
+}
+
+workflow test_srst2_srst2_single_end {
+
+    input = [
+        [ id:'test', single_end:true, db:"gene" ], // meta map
+        file(params.test_data['bacteroides_fragilis']['illumina']['test1_1_fastq_gz'], checkIfExists: true),
+        file('https://raw.githubusercontent.com/nf-core/test-datasets/modules/data/delete_me/srst2/resFinder_20180221_srst2.fasta')  // Change to params.test_data syntax after the data is included in tests/config/test_data.config
+    ]
+
+    SRST2_SRST2(input)
+}

tests/modules/srst2/srst2/nextflow.config

@@ -0,0 +1,5 @@
+process {
+
+    publishDir = { "${params.outdir}/${task.process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()}" }
+    
+}

tests/modules/srst2/srst2/test.yml

@@ -0,0 +1,51 @@
+- name: srst2 srst2 test_srst2_srst2_exit #Testing pipeline exit when not meta.db
+  command: nextflow run tests/modules/srst2/srst2 -entry test_srst2_srst2_exit -c tests/config/nextflow.config
+  tags:
+    - srst2/srst2
+    - srst2
+  exit_code: 1
+
+- name: srst2 srst2 test_srst2_srst2_mlst
+  command: nextflow run tests/modules/srst2/srst2 -entry test_srst2_srst2_mlst -c tests/config/nextflow.config
+  tags:
+    - srst2/srst2
+    - srst2
+  files:
+    - path: output/srst2/test__SRR9067271.MLST_DB.pileup
+      contains:
+        - "dnaJ-1	2	C	17	.........,.......	FFFFFFFFFFFFFFFFF"
+    - path: output/srst2/test__SRR9067271.MLST_DB.sorted.bam
+    - path: output/srst2/test__mlst__MLST_DB__results.txt
+      md5sum: ec1b1f69933401d67c57f64cad11a098
+    - path: output/srst2/versions.yml
+      md5sum: a0c256a2fd3636069710b8ef22ee5ea7
+
+- name: srst2 srst2 test_srst2_srst2_paired_end
+  command: nextflow run tests/modules/srst2/srst2 -entry test_srst2_srst2_paired_end -c tests/config/nextflow.config
+  tags:
+    - srst2/srst2
+    - srst2
+  files:
+    - path: output/srst2/test__genes__resFinder_20180221_srst2__results.txt
+      md5sum: 099aa6cacec5524b311f606debdfb3a9
+    - path: output/srst2/test__test1.resFinder_20180221_srst2.pileup
+      md5sum: 64b512ff495b828c456405ec7b676ad1
+    - path: output/srst2/test__test1.resFinder_20180221_srst2.sorted.bam
+    - path: output/srst2/versions.yml
+      md5sum: b446a70f1a2b4f60757829bcd744a214
+
+- name: srst2 srst2 test_srst2_srst2_single_end
+  command: nextflow run tests/modules/srst2/srst2 -entry test_srst2_srst2_single_end -c tests/config/nextflow.config
+  tags:
+    - srst2/srst2
+    - srst2
+  files:
+    - path: output/srst2/test__fullgenes__resFinder_20180221_srst2__results.txt
+      md5sum: d0762ef8c38afd0e0a34cce52ed1a3db
+    - path: output/srst2/test__genes__resFinder_20180221_srst2__results.txt
+      md5sum: b8850c6644406d8b131e471ecc3f9013
+    - path: output/srst2/test__test1_1.resFinder_20180221_srst2.pileup
+      md5sum: 5f6279dc8124aa762a9dfe3d7a871277
+    - path: output/srst2/test__test1_1.resFinder_20180221_srst2.sorted.bam
+    - path: output/srst2/versions.yml
+      md5sum: 790fe00493c6634d17801a930073218b

tests/modules/vardictjava/main.nf

@@ -0,0 +1,22 @@
+#!/usr/bin/env nextflow
+
+nextflow.enable.dsl = 2
+
+include { VARDICTJAVA } from '../../../modules/vardictjava/main.nf'
+
+workflow test_vardictjava {
+    
+    bam_input_ch = Channel.value([
+        [ id:'test' ], // meta map
+        file(params.test_data['homo_sapiens']['illumina']['test_paired_end_sorted_bam'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['illumina']['test_paired_end_sorted_bam_bai'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_bed'], checkIfExists: true)
+    ])
+
+    reference = Channel.value([
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta'], checkIfExists: true),
+        file(params.test_data['homo_sapiens']['genome']['genome_fasta_fai'], checkIfExists: true)
+    ])
+
+    VARDICTJAVA ( bam_input_ch, reference )
+}

tests/modules/vardictjava/nextflow.config

@@ -0,0 +1,5 @@
+process {
+
+    publishDir = { "${params.outdir}/${task.process.tokenize(':')[-1].tokenize('_')[0].toLowerCase()}" }
+    
+}

tests/modules/vardictjava/test.yml

@@ -0,0 +1,9 @@
+- name: vardictjava test_vardictjava
+  command: nextflow run tests/modules/vardictjava -entry test_vardictjava -c tests/config/nextflow.config
+  tags:
+    - vardictjava
+  files:
+    - path: output/vardictjava/test.vcf.gz
+      md5sum: 3f1f227afc532bddeb58f16fd3013fc8
+    - path: output/vardictjava/versions.yml
+      md5sum: 9b62c431a4f2680412b61c7071bdb1cd