Merge branch 'master' into motus_downloaddb

modules/antismash/antismashlitedownloaddatabases/main.nf

@@ -27,9 +27,7 @@ process ANTISMASH_ANTISMASHLITEDOWNLOADDATABASES {
 
     output:
     path("antismash_db") , emit: database
-    path("css"), emit: css_dir
-    path("detection"), emit: detection_dir
-    path("modules"), emit: modules_dir
+    path("antismash_dir"), emit: antismash_dir
     path "versions.yml", emit: versions
 
     when:
@@ -37,11 +35,19 @@ process ANTISMASH_ANTISMASHLITEDOWNLOADDATABASES {
 
     script:
     def args = task.ext.args ?: ''
+    conda = params.enable_conda
     """
     download-antismash-databases \\
         --database-dir antismash_db \\
         $args
 
+    if [[ $conda = false ]]; \
+        then \
+            cp -r /usr/local/lib/python3.8/site-packages/antismash antismash_dir; \
+        else \
+            cp -r \$(python -c 'import antismash;print(antismash.__file__.split("/__")[0])') antismash_dir; \
+    fi
+
     cat <<-END_VERSIONS > versions.yml
     "${task.process}":
         antismash-lite: \$(antismash --version | sed 's/antiSMASH //')

modules/antismash/antismashlitedownloaddatabases/meta.yml

@@ -50,21 +50,11 @@ output:
       type: directory
       description: Download directory for antiSMASH databases
       pattern: "antismash_db"
-  - css_dir:
+  - antismash_dir:
       type: directory
       description: |
-        antismash/outputs/html/css folder which is being created during the antiSMASH database downloading step. These files are normally downloaded by download-antismash-databases itself, and must be retrieved by the user by manually running the command with conda or a standalone installation of antiSMASH. Therefore we do not recommend using this module for production pipelines, but rather require users to specify their own local copy of the antiSMASH database in pipelines.
-      pattern: "css"
-  - detection_dir:
-      type: directory
-      description: |
-        antismash/detection folder which is being created during the antiSMASH database downloading step. These files are normally downloaded by download-antismash-databases itself, and must be retrieved by the user by manually running the command with conda or a standalone installation of antiSMASH. Therefore we do not recommend using this module for production pipelines, but rather require users to specify their own local copy of the antiSMASH database in pipelines.
-      pattern: "detection"
-  - modules_dir:
-      type: directory
-      description: |
-        antismash/modules folder which is being created during the antiSMASH database downloading step. These files are normally downloaded by download-antismash-databases itself, and must be retrieved by the user by manually running the command with conda or a standalone installation of antiSMASH. Therefore we do not recommend using this module for production pipelines, but rather require users to specify their own local copy of the antiSMASH database in pipelines.
-      pattern: "modules"
+        antismash installation folder which is being modified during the antiSMASH database downloading step. The modified files are normally downloaded by download-antismash-databases itself, and must be retrieved by the user by manually running the command with conda or a standalone installation of antiSMASH. Therefore we do not recommend using this module for production pipelines, but rather require users to specify their own local copy of the antiSMASH database and installation folder in pipelines.
+      pattern: "antismash_dir"
 
 authors:
   - "@jasmezz"

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/**
@@ -1723,6 +1731,10 @@ seqwish/induce:
   - modules/seqwish/induce/**
   - tests/modules/seqwish/induce/**
 
+shigatyper:
+  - modules/shigatyper/**
+  - tests/modules/shigatyper/**
+
 shovill:
   - modules/shovill/**
   - tests/modules/shovill/**
@@ -1731,6 +1743,10 @@ sistr:
   - modules/sistr/**
   - tests/modules/sistr/**
 
+slimfastq:
+  - modules/slimfastq/**
+  - tests/modules/slimfastq/**
+
 snapaligner/index:
   - modules/snapaligner/index/**
   - tests/modules/snapaligner/index/**
@@ -1779,6 +1795,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/**
@@ -1920,6 +1940,10 @@ unzip:
   - modules/unzip/**
   - tests/modules/unzip/**
 
+vardictjava:
+  - modules/vardictjava/**
+  - tests/modules/vardictjava/**
+
 variantbam:
   - modules/variantbam/**
   - tests/modules/variantbam/**

tests/modules/antismash/antismashlitedownloaddatabases/test.yml

@@ -1,8 +1,8 @@
 - name: antismash antismashlitedownloaddatabases test_antismash_antismashlitedownloaddatabases
   command: nextflow run tests/modules/antismash/antismashlitedownloaddatabases -entry test_antismash_antismashlitedownloaddatabases -c tests/config/nextflow.config
   tags:
-    - antismash
     - antismash/antismashlitedownloaddatabases
+    - antismash
   files:
     - path: output/antismash/versions.yml
       md5sum: 24859c67023abab99de295d3675a24b6
@@ -12,6 +12,5 @@
     - path: output/antismash/antismash_db/pfam
     - path: output/antismash/antismash_db/resfam
     - path: output/antismash/antismash_db/tigrfam
-    - path: output/antismash/css
-    - path: output/antismash/detection
-    - path: output/antismash/modules
+    - path: output/antismash/antismash_dir
+    - path: output/antismash/antismash_dir/detection/hmm_detection/data/bgc_seeds.hmm

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