12.15. De Novo Sequence Assembly - Clara Genomics Analysis

This is a reference pipeline using Clara Genomic Analysis tools to assemble genome with Clara Deploy SDK.

These tools exploit the abilities of GPU to accelerate gene sequencing.

12.15.1. Pipeline Definition

Copy
Copied!

            
            api-version: 0.4.0
name: denovo-gpu
parameters:
  DOCKER_IMAGE: claraparabricks/cga_cuda10
  DOCKER_TAG: v0.5.2

  MAPPER_KMER_SIZE: 15
  MAPPER_WINDOW_SIZE: 5
  MAPPER_ADDITIONAL_PARAMS: ''

  RACON_LOOPS: 5
  RACON_THREADS: 5
  RACON_POLISH_BATCH_SIZE: 4
  RACON_ADDITIONAL_PARAMS: ''
  JOB_ID: '.'

operators:
- name: mapper
  description: CUDA mapper
  container:
    image: ${{DOCKER_IMAGE}}
    tag: ${{DOCKER_TAG}}
    command: ["pef", "cudamapper",
              "-i", "/input",
              "-o", "/mapperOutput/${{JOB_ID}}",
              "-p", "-k${{MAPPER_KMER_SIZE}}-w${{MAPPER_WINDOW_SIZE}}"]
  requests:
    gpu: 1
    memory: 16384
  input:
  - path: /input/
  output:
  - path: /mapperOutput

- name: miniasm
  description: Miniasm
  container:
    image: ${{DOCKER_IMAGE}}
    tag: ${{DOCKER_TAG}}
    command: ["pef", "miniasm",
              "-i", "/input",
              "-l", "/mapperOutput/${{JOB_ID}}/overlaps.paf",
              "-o", "/asmOutput/${{JOB_ID}}"]
  input:
  - path: /input/
  - from: mapper
    path: /mapperOutput
  output:
  - path: /asmOutput
  requests:
    memory: 16384

- name: racon
  description: Polish Assembly using racon
  container:
    image: ${{DOCKER_IMAGE}}
    tag: ${{DOCKER_TAG}}
    command: ["pef", "racon",
              "-i", "/input",
              "-a", "/asmOutput/${{JOB_ID}}/reads.fa",
              "-o", "/raconOutput/${{JOB_ID}}",
              "-t", "${{RACON_THREADS}}",
              "-l", "${{RACON_LOOPS}}"]
  requests:
    cpu: 4
    gpu: 1
    memory: 16384
  input:
  - path: /input/
  - from: miniasm
    path: /asmOutput
  - from: mapper
    path: /mapperOutput
  output:
  - path: /raconOutput/

12.15.2. Executing the Pipeline

Please refer to the Run Reference Pipelines using Local Input Files in the How to run a Reference Pipeline section to learn how to register a pipeline and execute the pipeline using local input files.

12.15.3. Data Input

Input requires a folder containing the following files:

sample.fasta - Input fasta sample file for all-to-all mapping

jobConfig(optional) - A file containing param and value in shell script style. A sample (sample_job_config.sh) is provided. Following is content of a jobConfig file with default values.

Copy
Copied!

            
            KMER_SIZE=15         # length of kmer to use for minimizers
WINDOW_SIZE=5        # length of window to use for minimizers
INDEX_SIZE=10000     # length of batch size used for query

RACON_LOOPS=5        # Number of polishing loops
RACON_THREADS=15     # number of threads
POLISH_BATCH_SIZE=6  # number of batches for CUDA accelerated polishing

12.15.4. Data Output

Assembled and polished sequence