Sfoglia il codice sorgente

add annotate and faidx

master
YaqingLiu 4 anni fa
parent
commit
eece15f9b1
5 ha cambiato i file con 15 aggiunte e 14 eliminazioni
  1. +1
    -1
      defaults
  2. +2
    -2
      inputs
  3. +3
    -3
      tasks/access.wdl
  4. +5
    -4
      tasks/batch.wdl
  5. +4
    -4
      workflow.wdl

+ 1
- 1
defaults Vedi File

{ {
"fasta": "GRCh38.d1.vd1.fa", "fasta": "GRCh38.d1.vd1.fa",
"ref_dir": "oss://pgx-reference-data/GRCh38.d1.vd1/",
"faidx": "oss://pgx-reference-data/GRCh38.d1.vd1/",
"docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/cnvkit:0.9.7", "docker": "registry.cn-shanghai.aliyuncs.com/pgx-docker-registry/cnvkit:0.9.7",
"disk_size": "400", "disk_size": "400",
"cluster_config": "OnDemand bcs.a2.7xlarge img-ubuntu-vpc", "cluster_config": "OnDemand bcs.a2.7xlarge img-ubuntu-vpc",

+ 2
- 2
inputs Vedi File

"{{ project_name }}.tumor_bam": {{ tumor_bam | tojson }}, "{{ project_name }}.tumor_bam": {{ tumor_bam | tojson }},
"{{ project_name }}.tumor_bai": {{ tumor_bai | tojson }}, "{{ project_name }}.tumor_bai": {{ tumor_bai | tojson }},
"{{ project_name }}.fasta": "{{ fasta }}", "{{ project_name }}.fasta": "{{ fasta }}",
"{{ project_name }}.ref_dir": "{{ ref_dir }}",
"{{ project_name }}.ref_dir": "{{ ref_flat }}",
"{{ project_name }}.faidx": "{{ faidx }}",
"{{ project_name }}.faidx": "{{ ref_flat }}",
"{{ project_name }}.docker": "{{ docker }}", "{{ project_name }}.docker": "{{ docker }}",
"{{ project_name }}.bed": "{{ bed }}", "{{ project_name }}.bed": "{{ bed }}",
"{{ project_name }}.disk_size": "{{ disk_size }}", "{{ project_name }}.disk_size": "{{ disk_size }}",

+ 3
- 3
tasks/access.wdl Vedi File

task access { task access {
File ref_dir
String fasta
File faidx
File fasta
String docker String docker
String cluster_config String cluster_config
String disk_size String disk_size


command <<< command <<<
cnvkit.py access -s 5000 -o access-5kb-mappable.hg38.bed ${ref_dir}/${fasta}
cnvkit.py access -s 5000 -o access-5kb-mappable.hg38.bed ${fasta}
>>> >>>


runtime { runtime {

+ 5
- 4
tasks/batch.wdl Vedi File

Array[File] normal_bam Array[File] normal_bam
Array[File] normal_bai Array[File] normal_bai
File bed File bed
File ref_dir
String fasta
File faidx
File fasta
File ref_flat File ref_flat
File access_bed File access_bed
String docker String docker
mkdir -p /cromwell_root/tmp/cnvkit mkdir -p /cromwell_root/tmp/cnvkit
cp ${sep=' ' normal_bai} /cromwell_root/tmp/cnvkit cp ${sep=' ' normal_bai} /cromwell_root/tmp/cnvkit
cp ${sep=' ' tumor_bai} /cromwell_root/tmp/cnvkit cp ${sep=' ' tumor_bai} /cromwell_root/tmp/cnvkit
cp ${faidx} /cromwell_root/tmp/cnvkit
cd /cromwell_root/tmp/cnvkit cd /cromwell_root/tmp/cnvkit
mkdir results mkdir results
cnvkit.py batch faidx --no-rebuild${sep=' ' tumor_bam} --normal ${sep=' ' normal_bam} \
cnvkit.py batch ${sep=' ' tumor_bam} --normal ${sep=' ' normal_bam} \
--targets ${bed} --annotate ${ref_flat} \ --targets ${bed} --annotate ${ref_flat} \
--fasta ${ref_dir}/${fasta} --access ${access_bed} \
--fasta ${fasta} --access ${access_bed} \
--output-reference ${sample_id}.reference.cnn --output-dir /cromwell_root/tmp/cnvkit/results/ \ --output-reference ${sample_id}.reference.cnn --output-dir /cromwell_root/tmp/cnvkit/results/ \
--diagram --scatter --diagram --scatter
>>> >>>

+ 4
- 4
workflow.wdl Vedi File

Array[File] normal_bam Array[File] normal_bam
Array[File] normal_bai Array[File] normal_bai
File bed File bed
File ref_dir
String fasta
File faidx
File fasta
File ref_flat File ref_flat
String docker String docker
String cluster_config String cluster_config
call access.access as access { call access.access as access {
input: input:
fasta = fasta, fasta = fasta,
ref_dir = ref_dir,
faidx = faidx,
docker = docker, docker = docker,
cluster_config = cluster_config, cluster_config = cluster_config,
disk_size = disk_size disk_size = disk_size
input: input:
sample_id = sample_id, sample_id = sample_id,
fasta = fasta, fasta = fasta,
ref_dir = ref_dir,
faidx = faidx,
ref_flat = ref_flat, ref_flat = ref_flat,
tumor_bam = tumor_bam, tumor_bam = tumor_bam,
tumor_bai = tumor_bai, tumor_bai = tumor_bai,

Loading…
Annulla
Salva