Skip to content

kDifferentialExpression usecase isn't working #88

@mr-eyes

Description

@mr-eyes

I followed all the steps in the usecase directory to run the snakemake but it's failing.

./downloadData.sh: 7: Syntax error: "(" unexpected

snakemake -j2 --use-conda
Building DAG of jobs...
Using shell: /usr/bin/bash
Provided cores: 2
Rules claiming more threads will be scaled down.
Job stats:
job                       count    min threads    max threads
----------------------  -------  -------------  -------------
all                           1              1              1
kDifferntialExpression        1              1              1
kmc                           2              2              2
total                         4              1              2

Select jobs to execute...

[Tue Aug 31 04:21:50 2021]
rule kmc:
    input: inputData/HBR_Rep1_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz, inputData/HBR_Rep1_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz, inputData/HBR_Rep2_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz, inputData/HBR_Rep2_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz, inputData/HBR_Rep3_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz, inputData/HBR_Rep3_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz
    output: sample.kmc_pre, sample.kmc_suf, sample.kmc.time
    log: sample.kmc.log
    jobid: 2
    wildcards: experiment=sample
    threads: 2
    resources: tmpdir=/tmp

Activating conda environment: /home/mabuelanin/dib-dev/kProcessor_wrapper/build/usecases/kDifferentialExpression/.snakemake/conda/cb6edb19c50825bf280817fbb4ac06b5
[Tue Aug 31 04:21:50 2021]
Error in rule kmc:
    jobid: 2
    output: sample.kmc_pre, sample.kmc_suf, sample.kmc.time
    log: sample.kmc.log (check log file(s) for error message)
    conda-env: /home/mabuelanin/dib-dev/kProcessor_wrapper/build/usecases/kDifferentialExpression/.snakemake/conda/cb6edb19c50825bf280817fbb4ac06b5
    shell:

        echo inputData/HBR_Rep1_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz inputData/HBR_Rep1_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz inputData/HBR_Rep2_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz inputData/HBR_Rep2_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz inputData/HBR_Rep3_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz inputData/HBR_Rep3_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz |tr -s ' ' $'
' > sample.lst
        /usr/bin/time -v -o sample.kmc.time kmc -ci1 -t2 -k31 -m2  @sample.lst sample  ./ &> sample.kmc.log

        (one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)

Removing output files of failed job kmc since they might be corrupted:
sample.kmc.time
Shutting down, this might take some time.
Exiting because a job execution failed. Look above for error message
Complete log: /home/mabuelanin/dib-dev/kProcessor_wrapper/build/usecases/kDifferentialExpression/.snakemake/log/2021-08-31T042149.582856.snakemake.log

Full Log


Building DAG of jobs...
Using shell: /usr/bin/bash
Provided cores: 10
Rules claiming more threads will be scaled down.
Conda environments: ignored
Job stats:
job                       count    min threads    max threads
----------------------  -------  -------------  -------------
all                           1              1              1
kDifferntialExpression        1              1              1
kmc                           2              4              4
total                         4              1              4

Select jobs to execute...

[Tue Aug 31 04:00:51 2021]
rule kmc:
    input: inputData/UHR_Rep1_ERCC-Mix1_Build37-ErccTranscripts-chr22.read1.fastq.gz, inputData/UHR_Rep1_ERCC-Mix1_Build37-ErccTranscripts-chr22.read2.fastq.gz, inputData/UHR_Rep2_ERCC-Mix1_Build37-ErccTranscripts-chr22.read1.fastq.gz, inputData/UHR_Rep2_ERCC-Mix1_Build37-ErccTranscripts-chr22.read2.fastq.gz, inputData/UHR_Rep3_ERCC-Mix1_Build37-ErccTranscripts-chr22.read1.fastq.gz, inputData/UHR_Rep3_ERCC-Mix1_Build37-ErccTranscripts-chr22.read2.fastq.gz
    output: control.kmc_pre, control.kmc_suf, control.kmc.time
    log: control.kmc.log
    jobid: 3
    wildcards: experiment=control
    threads: 4
    resources: tmpdir=/tmp


[Tue Aug 31 04:00:51 2021]
rule kmc:
    input: inputData/HBR_Rep1_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz, inputData/HBR_Rep1_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz, inputData/HBR_Rep2_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz, inputData/HBR_Rep2_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz, inputData/HBR_Rep3_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz, inputData/HBR_Rep3_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz
    output: sample.kmc_pre, sample.kmc_suf, sample.kmc.time
    log: sample.kmc.log
    jobid: 2
    wildcards: experiment=sample
    threads: 4
    resources: tmpdir=/tmp

[Tue Aug 31 04:00:51 2021]
Error in rule kmc:
    jobid: 3
    output: control.kmc_pre, control.kmc_suf, control.kmc.time
    log: control.kmc.log (check log file(s) for error message)
    shell:
        
	echo inputData/UHR_Rep1_ERCC-Mix1_Build37-ErccTranscripts-chr22.read1.fastq.gz inputData/UHR_Rep1_ERCC-Mix1_Build37-ErccTranscripts-chr22.read2.fastq.gz inputData/UHR_Rep2_ERCC-Mix1_Build37-ErccTranscripts-chr22.read1.fastq.gz inputData/UHR_Rep2_ERCC-Mix1_Build37-ErccTranscripts-chr22.read2.fastq.gz inputData/UHR_Rep3_ERCC-Mix1_Build37-ErccTranscripts-chr22.read1.fastq.gz inputData/UHR_Rep3_ERCC-Mix1_Build37-ErccTranscripts-chr22.read2.fastq.gz |tr -s ' ' $'
' > control.lst
	/usr/bin/time -v -o control.kmc.time kmc -ci1 -t4 -k31 -m2  @control.lst control  ./ &> control.kmc.log
        
        (one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)

Removing output files of failed job kmc since they might be corrupted:
control.kmc.time
[Tue Aug 31 04:00:51 2021]
Error in rule kmc:
    jobid: 2
    output: sample.kmc_pre, sample.kmc_suf, sample.kmc.time
    log: sample.kmc.log (check log file(s) for error message)
    shell:
        
	echo inputData/HBR_Rep1_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz inputData/HBR_Rep1_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz inputData/HBR_Rep2_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz inputData/HBR_Rep2_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz inputData/HBR_Rep3_ERCC-Mix2_Build37-ErccTranscripts-chr22.read1.fastq.gz inputData/HBR_Rep3_ERCC-Mix2_Build37-ErccTranscripts-chr22.read2.fastq.gz |tr -s ' ' $'
' > sample.lst
	/usr/bin/time -v -o sample.kmc.time kmc -ci1 -t4 -k31 -m2  @sample.lst sample  ./ &> sample.kmc.log
        
        (one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)

Removing output files of failed job kmc since they might be corrupted:
sample.kmc.time
Shutting down, this might take some time.
Exiting because a job execution failed. Look above for error message
Complete log: /home/mabuelanin/dib-dev/kProcessor_wrapper/build/usecases/kDifferentialExpression/.snakemake/log/2021-08-31T040051.177052.snakemake.log


Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions