Skip to content

Commit

Permalink
Added sample tools.groovy sample + portability fixes
Browse files Browse the repository at this point in the history
- avoid zcat due to OS X issue with syntax
- define type for all run_salmon stages (otherwise
  breaks with newer bpipe versions)
- remove quant type from create_salmon_index (legacy code)
  • Loading branch information
mcmero committed Sep 27, 2021
1 parent f2c784b commit aa5eab9
Show file tree
Hide file tree
Showing 3 changed files with 16 additions and 6 deletions.
2 changes: 1 addition & 1 deletion .gitignore
Original file line number Diff line number Diff line change
@@ -1,4 +1,4 @@
tools*
tools/*
ref*
*/__pycache__
allvars-case*
Expand Down
7 changes: 2 additions & 5 deletions MINTIE.groovy
Original file line number Diff line number Diff line change
Expand Up @@ -92,7 +92,7 @@ assemble = {
""", "assemble"
} else {
exec """
rlens=`zcat $input1 $input2 \
rlens=`gunzip -c $input1 $input2 \
| awk -v mrl=$min_read_length 'BEGIN {minlen = mrl; maxlen = 0} {
if (NR % 4 == 2) {
rlen = length(\$1) ;
Expand Down Expand Up @@ -138,9 +138,6 @@ assemble = {
create_salmon_index = {
def sample_name = branch.name
def salmon_index = sample_name + "/all_fasta_index"
if (type == "quant") {
salmon_index = sample_name + "/salmon_quant_index"
}
output.dir = salmon_index
def index_fasta = output.dir + "/" + sample_name + ".fasta"
produce(index_fasta, '*.bin'){
Expand Down Expand Up @@ -307,7 +304,7 @@ run { fastqCaseFormat * [ fastq_dedupe +
trim +
assemble +
create_salmon_index +
[ fastqCaseFormat * [ run_salmon ],
[ fastqCaseFormat * [ run_salmon.using(type: "case") ],
fastqControlFormat * [ run_salmon.using(type:"controls") ] ] +
create_ec_count_matrix +
run_de +
Expand Down
13 changes: 13 additions & 0 deletions tools.groovy
Original file line number Diff line number Diff line change
@@ -0,0 +1,13 @@
// Path to tools used by the MINTIE pipeline
bpipe="bpipe"
fastuniq="fastuniq"
dedupe="dedupe.sh"
trimmomatic="trimmomatic"
fasta_formatter="fasta_formatter"
samtools="samtools"
bedtools="bedtools"
soapdenovotrans="SOAPdenovo-Trans-127mer"
salmon="salmon"
gmap="gmap"
R="R"
python="python"

0 comments on commit aa5eab9

Please sign in to comment.