From: Ward Vandewege <ward@curii.com>
Date: Tue, 30 Nov 2021 16:56:31 +0000 (-0500)
Subject: Add the star-index example for the debugging blog post.
X-Git-Url: https://git.arvados.org/arvados-tutorial.git/commitdiff_plain/4292ebf16822fde6b023ca1828bcae97b6d812ce

Add the star-index example for the debugging blog post.

refs #18178

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

diff --git a/blog-examples/debugging/fastqc_example/fastqc_blog.cwl b/blog-examples/debugging/fastqc/fastqc_blog.cwl
similarity index 100%
rename from blog-examples/debugging/fastqc_example/fastqc_blog.cwl
rename to blog-examples/debugging/fastqc/fastqc_blog.cwl
diff --git a/blog-examples/debugging/fastqc_example/fastqc_blog.yml b/blog-examples/debugging/fastqc/fastqc_blog.yml
similarity index 100%
rename from blog-examples/debugging/fastqc_example/fastqc_blog.yml
rename to blog-examples/debugging/fastqc/fastqc_blog.yml
diff --git a/blog-examples/debugging/fastqc_example/fastqc_blog_broken.cwl b/blog-examples/debugging/fastqc/fastqc_blog_broken.cwl
similarity index 100%
rename from blog-examples/debugging/fastqc_example/fastqc_blog_broken.cwl
rename to blog-examples/debugging/fastqc/fastqc_blog_broken.cwl
diff --git a/blog-examples/debugging/fastqc_example/fastqc_blog_mismatch.yml b/blog-examples/debugging/fastqc/fastqc_blog_mismatch.yml
similarity index 100%
rename from blog-examples/debugging/fastqc_example/fastqc_blog_mismatch.yml
rename to blog-examples/debugging/fastqc/fastqc_blog_mismatch.yml
diff --git a/blog-examples/debugging/star-index/.gitignore b/blog-examples/debugging/star-index/.gitignore
new file mode 100644
index 0000000..681307f
--- /dev/null
+++ b/blog-examples/debugging/star-index/.gitignore
@@ -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
index 0000000..70872ac
--- /dev/null
+++ b/blog-examples/debugging/star-index/README.md
@@ -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
index 0000000..10b3fa1
--- /dev/null
+++ b/blog-examples/debugging/star-index/STAR-index-broken.cwl
@@ -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
index 0000000..74ee7db
--- /dev/null
+++ b/blog-examples/debugging/star-index/STAR-index.cwl
@@ -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
index 0000000..189faf2
--- /dev/null
+++ b/blog-examples/debugging/star-index/STAR-index.yml
@@ -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