Add the star-index example for the debugging blog post. blog-examples
authorWard Vandewege <ward@curii.com>
Tue, 30 Nov 2021 16:56:31 +0000 (11:56 -0500)
committerWard Vandewege <ward@jhvc.com>
Tue, 30 Nov 2021 17:35:50 +0000 (12:35 -0500)
refs #18178

Arvados-DCO-1.1-Signed-off-by: Ward Vandewege <ward@curii.com>

blog-examples/debugging/fastqc/fastqc_blog.cwl [moved from blog-examples/debugging/fastqc_example/fastqc_blog.cwl with 100% similarity]
blog-examples/debugging/fastqc/fastqc_blog.yml [moved from blog-examples/debugging/fastqc_example/fastqc_blog.yml with 100% similarity]
blog-examples/debugging/fastqc/fastqc_blog_broken.cwl [moved from blog-examples/debugging/fastqc_example/fastqc_blog_broken.cwl with 100% similarity]
blog-examples/debugging/fastqc/fastqc_blog_mismatch.yml [moved from blog-examples/debugging/fastqc_example/fastqc_blog_mismatch.yml with 100% similarity]
blog-examples/debugging/star-index/.gitignore [new file with mode: 0644]
blog-examples/debugging/star-index/README.md [new file with mode: 0644]
blog-examples/debugging/star-index/STAR-index-broken.cwl [new file with mode: 0755]
blog-examples/debugging/star-index/STAR-index.cwl [new file with mode: 0755]
blog-examples/debugging/star-index/STAR-index.yml [new file with mode: 0644]

diff --git a/blog-examples/debugging/star-index/.gitignore b/blog-examples/debugging/star-index/.gitignore
new file mode 100644 (file)
index 0000000..681307f
--- /dev/null
@@ -0,0 +1 @@
+rnaseq
diff --git a/blog-examples/debugging/star-index/README.md b/blog-examples/debugging/star-index/README.md
new file mode 100644 (file)
index 0000000..70872ac
--- /dev/null
@@ -0,0 +1,16 @@
+
+To test this:
+
+  mkdir rnaseq
+  cd rnaseq
+  wget --mirror --no-parent --no-host --cut-dirs=1 https://download.jutro.arvadosapi.com/c=9178fe1b80a08a422dbe02adfd439764+925/
+  cd ..
+
+Then:
+
+   arvados-cwl-runner STAR-index-broken.cwl STAR-index.yml
+
+or:
+
+   arvados-cwl-runner STAR-index.cwl STAR-index.yml
+
diff --git a/blog-examples/debugging/star-index/STAR-index-broken.cwl b/blog-examples/debugging/star-index/STAR-index-broken.cwl
new file mode 100755 (executable)
index 0000000..10b3fa1
--- /dev/null
@@ -0,0 +1,58 @@
+#!/usr/bin/env cwl-runner
+cwlVersion: v1.0
+class: CommandLineTool
+
+requirements:
+  DockerRequirement:
+    dockerPull: "quay.io/biocontainers/star:2.7.5c--0"
+
+inputs:
+
+  InputFiles:
+    format: http://edamontology.org/format_1930
+    type: File[]
+    inputBinding:
+      prefix: "--genomeFastaFiles"
+
+  IndexName:
+    type: string
+    inputBinding:
+      prefix: "--genomeDir"
+      valueFrom: ./$(self)
+
+#Optional Inputs
+
+  Gtf:
+    type: File?
+    inputBinding:
+      prefix: "--sjdbGTFfile"
+
+  Overhang:
+    type: int?
+    inputBinding:
+      prefix: "--sjdbOverhang"
+
+  Junctions:
+    type: File?
+    inputBinding:
+      prefix: "--sjdbFileChrStartEnd"
+
+  GenomeSize:
+    type: int?
+    inputBinding:
+      prefix: "--genomeSAindexNbases"
+
+  GenomeBits:
+    type: int?
+    inputBinding:
+      prefix: "--genomeChrBinNbits"
+
+baseCommand: [STAR, --runMode, genomeGenerate]
+
+arguments: [--runThreadN, $(runtime.cores)]
+
+outputs:
+  indexes:
+    type: Directory
+    outputBinding:
+      glob: ./$(inputs.IndexName)/
diff --git a/blog-examples/debugging/star-index/STAR-index.cwl b/blog-examples/debugging/star-index/STAR-index.cwl
new file mode 100755 (executable)
index 0000000..74ee7db
--- /dev/null
@@ -0,0 +1,61 @@
+#!/usr/bin/env cwl-runner
+cwlVersion: v1.0
+class: CommandLineTool
+
+requirements:
+  DockerRequirement:
+    dockerPull: "quay.io/biocontainers/star:2.7.5c--0"
+  ResourceRequirement:
+    coresMin: 2
+    ramMin: 12000
+
+inputs:
+
+  InputFiles:
+    format: http://edamontology.org/format_1930
+    type: File[]
+    inputBinding:
+      prefix: "--genomeFastaFiles"
+
+  IndexName:
+    type: string
+    inputBinding:
+      prefix: "--genomeDir"
+      valueFrom: ./$(self)
+
+#Optional Inputs
+
+  Gtf:
+    type: File?
+    inputBinding:
+      prefix: "--sjdbGTFfile"
+
+  Overhang:
+    type: int?
+    inputBinding:
+      prefix: "--sjdbOverhang"
+
+  Junctions:
+    type: File?
+    inputBinding:
+      prefix: "--sjdbFileChrStartEnd"
+
+  GenomeSize:
+    type: int?
+    inputBinding:
+      prefix: "--genomeSAindexNbases"
+
+  GenomeBits:
+    type: int?
+    inputBinding:
+      prefix: "--genomeChrBinNbits"
+
+baseCommand: [STAR, --runMode, genomeGenerate]
+
+arguments: [--runThreadN, $(runtime.cores)]
+
+outputs:
+  indexes:
+    type: Directory
+    outputBinding:
+      glob: ./$(inputs.IndexName)/
diff --git a/blog-examples/debugging/star-index/STAR-index.yml b/blog-examples/debugging/star-index/STAR-index.yml
new file mode 100644 (file)
index 0000000..189faf2
--- /dev/null
@@ -0,0 +1,9 @@
+InputFiles:
+  - class: File
+    location: rnaseq/reference_data/chr1.fa
+    format: http://edamontology.org/format_1930
+IndexName: 'hg19-chr1-STAR-index'
+Gtf:
+  class: File
+  location: rnaseq/reference_data/chr1-hg19_genes.gtf
+Overhang: 99