Skip to content
Open
Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
33 commits
Select commit Hold shift + click to select a range
d49aac1
added first draft for picard CollectMultipleMetrics
Mar 5, 2018
d7bcac1
- copied samtools files from deepWGBS
WiebkeSchmitt Mar 7, 2018
e21da91
ádded bwa files
WiebkeSchmitt Mar 7, 2018
be3d040
corrected bwa files
WiebkeSchmitt Mar 7, 2018
23adab9
started samtools-sort
WiebkeSchmitt Mar 7, 2018
4850c3c
´merge fix for bioconda-tool-bwa-index.cwl
WiebkeSchmitt Mar 8, 2018
d6eb571
samtools-sort now uses -o as default option
WiebkeSchmitt Mar 8, 2018
69fdc60
first draft for GAL workflow
WiebkeSchmitt Mar 8, 2018
d9d27ff
Merge branch 'master' of github.com:deepPipelines/deepChIPseq into Wi…
WiebkeSchmitt Mar 12, 2018
037c972
first draf of GAL pipeline
WiebkeSchmitt Mar 12, 2018
dc8a1f2
bwa-aln parameters are now set
WiebkeSchmitt Mar 12, 2018
67071ed
-added some parameters to tools -refined workflow
WiebkeSchmitt Mar 13, 2018
cfbad15
job file for GAL-workflow added
WiebkeSchmitt Mar 13, 2018
b04a0e5
corrected GAL workflow and adjusted job file accordingly
WiebkeSchmitt Mar 14, 2018
4ec556c
bwa-aln step and bwa-sampe step are now working, error in samtools-vi…
WiebkeSchmitt Mar 15, 2018
54095c1
- corrected prefix and reference genome in job file - corrected workf…
WiebkeSchmitt Mar 15, 2018
0239fdc
name for bwa-aln works now
WiebkeSchmitt Mar 15, 2018
00ad8bb
debugged workflow, error in CollectMultipleMetrics still needs fixing
WiebkeSchmitt Mar 16, 2018
2237a8e
workflow is now running for testdata in test folder
WiebkeSchmitt Mar 16, 2018
4a227c0
workflow now outputs .bam file from picardMarkDuplicates step
WiebkeSchmitt Mar 16, 2018
1eb0636
started dividing GAL workflow into two parts for parallelization
WiebkeSchmitt Mar 19, 2018
9ca00cb
separated parts of GAL workflow (part 1 and part 2) are now validated
WiebkeSchmitt Mar 20, 2018
50a7ba9
started parallel workflow for GAL workflow
WiebkeSchmitt Mar 20, 2018
0dd9c76
first draft of parallel workflow for GAL is now vailid
WiebkeSchmitt Mar 22, 2018
2d36f42
started scattering the GAL workflow
WiebkeSchmitt Mar 22, 2018
7f1c551
parallel workflow now valid and running, but with error in PicardMark…
WiebkeSchmitt Mar 23, 2018
a00f537
parallel workflow for GAL is now running for sample data
WiebkeSchmitt Mar 23, 2018
2438490
parallel worflow now outputs bam-File from Picard-MarkDuplicates step
WiebkeSchmitt Mar 23, 2018
3ad0f71
secondaryFiles-Option now in bwa-aln tool
WiebkeSchmitt Mar 23, 2018
8d558eb
- changed tool for bwa sampe to allow usage of secondary files - GAL-…
WiebkeSchmitt Mar 26, 2018
5ac805f
modified bwa-index tool to create secondaryFiles as output, testing y…
WiebkeSchmitt Mar 26, 2018
611f6b4
bwa-index is now running without error, but still without output
WiebkeSchmitt Mar 27, 2018
cce084e
bwa index tool is debugged further, but still not running as desired
WiebkeSchmitt Mar 28, 2018
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
12 changes: 8 additions & 4 deletions CWL/tools/bioconda-tool-bwa-aln.yml
Original file line number Diff line number Diff line change
Expand Up @@ -33,7 +33,13 @@ outputs:
inputs:

prefix:
type: string
type: File
secondaryFiles:
- ".amb"
- ".ann"
- ".bwt"
- ".pac"
- ".sa"
doc: |
reference prefix
inputBinding:
Expand Down Expand Up @@ -84,7 +90,7 @@ inputs:
prefix: "-e"

minIndelEndDist:
type: int?]
type: int?
doc: |
do not put an indel within INT bp towards the ends [5]
inputBinding:
Expand Down Expand Up @@ -139,8 +145,6 @@ inputs:
position: 1
prefix: "-q"



doc: |
Usage: bwa aln [options] <prefix> <in.fq>

Expand Down
153 changes: 82 additions & 71 deletions CWL/tools/bioconda-tool-bwa-index.cwl
Original file line number Diff line number Diff line change
@@ -1,15 +1,20 @@
#!/usr/bin/env cwl-runner

cwlVersion: "v1.0"

class: CommandLineTool

cwlVersion: "v1.0"

s:author:
- class: s:Person
s:identifier: https://orcid.org/0000-0001-6231-4417
s:email: mailto:karl.nordstroem@uni-saarland.de
s:name: Karl Nordström

requirements:
- class: InlineJavascriptRequirement
- class: InitialWorkDirRequirement
listing: input

hints:
- class: ResourceRequirement
coresMin: 1
Expand All @@ -18,113 +23,121 @@ hints:
- class: DockerRequirement
dockerPull: "quay.io/biocontainers/bwa:0.6.2--1"

requirements:
- class: InlineJavascriptRequirement

inputs:

input:
type: File
inputBinding:
position: 3

algorithm:
type: string?
description: |
doc: |
BWT construction algorithm: bwtsw or is (Default: auto)
inputBinding:
position: 2
prefix: "-a"

outPrefix:
type: string?
description: |
doc: |
Prefix of the index (Default: same as fasta name)
inputBinding:
position: 2
prefix: "-p"

blockSize:
type: int?
description: |
doc: |
Block size for the bwtsw algorithm (effective with -a bwtsw) (Default: 10000000)
inputBinding:
position: 2
prefix: "-b"

altSuffix:
type: boolean?
description: |
doc: |
Index files named as <in.fasta>.64.* instead of <in.fasta>.*
inputBinding:
position: 2
prefix: "-6"

outputs:
- id: output
type: { type: array, items: File }
outputBinding:
glob:
- ${
if (inputs.p) {
return inputs.p + ".amb"
} else {
if (inputs._6 == true) {
return inputs.input.path + ".64.amb"
} else {
return inputs.input.path + ".amb"
}
}
}
- ${
if (inputs.p) {
return inputs.p + ".ann"
} else {
if (inputs._6 == true) {
return inputs.input.path + ".64.ann"
} else {
return inputs.input.path + ".ann"
}
}
}
- ${
if (inputs.p) {
return inputs.p + ".bwt"
} else {
if (inputs._6 == true) {
return inputs.input.path + ".64.bwt"
} else {
return inputs.input.path + ".bwt"
}
}
}
- ${
if (inputs.p) {
return inputs.p + ".pac"
} else {
if (inputs._6 == true) {
return inputs.input.path + ".64.pac"
} else {
return inputs.input.path + ".pac"
}
}
}
- ${
if (inputs.p) {
return inputs.p + ".sa"
} else {
if (inputs._6 == true) {
return inputs.input.path + ".64.sa"
} else {
return inputs.input.path + ".sa"
}
}
}
# - id: output
# index:
# type: File
# outputBinding:
# glob: $( inputs.input )
# secondaryFiles:
# - ".amb"
# - ".ann"
# - ".bwt"
# - ".pac"
# - ".sa"
# type: { type: array, items: File }
# outputBinding:
# glob:
# - ${
# if (inputs.p) {
# return inputs.p + ".amb"
# } else {
# if (inputs._6 == true) {
# return inputs.input.path + ".64.amb"
# } else {
# return inputs.input.path + ".amb"
# }
# }
# }
# - ${
# if (inputs.p) {
# return inputs.p + ".ann"
# } else {
# if (inputs._6 == true) {
# return inputs.input.path + ".64.ann"
# } else {
# return inputs.input.path + ".ann"
# }
# }
# }
# - ${
# if (inputs.p) {
# return inputs.p + ".bwt"
# } else {
# if (inputs._6 == true) {
# return inputs.input.path + ".64.bwt"
# } else {
# return inputs.input.path + ".bwt"
# }
# }
# }
# - ${
# if (inputs.p) {
# return inputs.p + ".pac"
# } else {
# if (inputs._6 == true) {
# return inputs.input.path + ".64.pac"
# } else {
# return inputs.input.path + ".pac"
# }
# }
# }
# - ${
# if (inputs.p) {
# return inputs.p + ".sa"
# } else {
# if (inputs._6 == true) {
# return inputs.input.path + ".64.sa"
# } else {
# return inputs.input.path + ".sa"
# }
# }
# }

baseCommand:
- bwa
- index

description: |
doc: |
Usage: bwa index [options] <in.fasta>

Options: -a STR BWT construction algorithm: bwtsw or is [auto]
Expand All @@ -135,8 +148,6 @@ description: |
Warning: `-a bwtsw' does not work for short genomes, while `-a is' and
`-a div' do not work not for long genomes.



$namespaces:
s: https://schema.org/
edam: http://edamontology.org/
Expand Down
24 changes: 22 additions & 2 deletions CWL/tools/bioconda-tool-bwa-sampe.yml
Original file line number Diff line number Diff line change
Expand Up @@ -35,7 +35,13 @@ outputs:
inputs:

prefix:
type: string
type: File
secondaryFiles:
- ".amb"
- ".ann"
- ".bwt"
- ".pac"
- ".sa"
doc: |
reference prefix
inputBinding:
Expand Down Expand Up @@ -71,6 +77,14 @@ inputs:
inputBinding:
position: 14

maximumInsertSize:
type: int?
inputBinding:
position: 5
prefix: "-a"
doc: |
maximum nsert size [500]

noSW:
type: boolean?
doc: |
Expand Down Expand Up @@ -103,7 +117,13 @@ inputs:
position: 5
prefix: "-N"


readGroupHeaderLine:
type: string?
inputBinding:
position: 5
prefix: "-r"
doc: |
read group header line such as '@RG\tID:foo\tSM:bar' [null]

doc: |
Usage: bwa sampe [options] <prefix> <in1.sai> <in2.sai> <in1.fq> <in2.fq>
Expand Down
Loading