This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to run SPAdes inside docker?

Hi, Any running example in cwl script to run SPAdes inside docker? I can run SPAdes without docker but would like to know how to run with docker.

cwl docker container
• 2,873 views • link
updated by Michael R. Crusoe Michael R. Crusoe • written by jeetsahu • jeetsahu

Hello jeetsahu, can you share your CWL tool description so we can help you?

• 1 views • link
written by Michael R. Crusoe Michael R. Crusoe

1 answer

Here is an example CWL description that uses SPAdes and a Docker software container from biocontainers.pro:

https://github.com/EBI-Metagenomics/CWL-assembly/blob/38738ecc25412a1614241ed9fddd063ea0843561/cwl/assembly/spades.cwl

#!/usr/bin/env cwl-runner
cwlVersion: v1.0
class: CommandLineTool


label: "spades: de novo genomic assembler"

requirements:
  DockerRequirement:
    dockerPull: "quay.io/biocontainers/spades:3.12.0--1"
  InlineJavascriptRequirement: {}

hints:
  SoftwareRequirement:
    packages:
      spades:
        specs: [ "https://identifiers.org/rrid/RRID:SCR_000131" ]
        version: [ "3.12.0" ]

baseCommand: [ spades.py ]

arguments:
  - valueFrom: $(runtime.outdir)
    prefix: -o
#  - valueFrom: $(runtime.tmpdir)
#    prefix: --tmp-dir
  - valueFrom: $(runtime.ram/1000)
    prefix: --memory
  - valueFrom: $(runtime.cores)
    prefix: --threads

inputs:
  forward_reads:
    type: File?
#     format: edam:format_1930  # FASTQ
    inputBinding:
      prefix: "-1"
  reverse_reads:
    type: File?
#     format: edam:format_1930  # FASTQ
    inputBinding:
      prefix: "-2"
  single_reads:
    type: File?
#     format: edam:format_1930  # FASTQ
    inputBinding:
      prefix: "-s"
  interleaved_reads:
    type: File?
    #     format: edam:format_1930  # FASTQ
    inputBinding:
      prefix: "--12"

stdout: stdout.txt
stderr: stderr.txt

outputs:
  stdout: stdout
  stderr: stderr
  contigs:
    type: File
    format: edam:format_1929  # FASTA
    outputBinding:
      glob: contigs.fasta

  scaffolds:
    type: File
    format: edam:format_1929  # FASTA
    outputBinding:
      glob: scaffolds.fasta

  assembly_graph:
    type: File
    #format: edam:format_TBD  # FASTG
    outputBinding:
      glob: assembly_graph.fastg

  contigs_assembly_graph:
    type: File
    outputBinding:
      glob: contigs.paths

  scaffolds_assembly_graph:
    type: File
    outputBinding:
      glob: scaffolds.paths

  contigs_before_rr:
    label: contigs before repeat resolution
    type: File
    format: edam:format_1929  # FASTA
    outputBinding:
      glob: before_rr.fasta

  params:
    label: information about SPAdes parameters in this run
    type: File
    format: iana:text/plain
    outputBinding:
      glob: params.txt

  log:
    label: SPAdes log
    type: File
    format: iana:text/plain
    outputBinding:
      glob: spades.log

  internal_config:
    label: internal configuration file
    type: File
    # format: text/plain
    outputBinding:
      glob: dataset.info

  internal_dataset:
    label: internal YAML data set file
    type: File
    outputBinding:
      glob: input_dataset.yaml

$namespaces:
 edam: http://edamontology.org/
 iana: https://www.iana.org/assignments/media-types/
 s: http://schema.org/
$schemas:
 - http://edamontology.org/EDAM_1.16.owl
 - https://schema.org/docs/schema_org_rdfa.html

s:license: "https://www.apache.org/licenses/LICENSE-2.0"
s:copyrightHolder: "EMBL - European Bioinformatics Institute"

doc: |
  https://arxiv.org/abs/1604.03071
  http://cab.spbu.ru/files/release3.12.0/manual.html#meta
• 1 views • link
written by Michael R. Crusoe Michael R. Crusoe

Log in to answer this question.

More posts like this