отправка задания создания змеи в кластер из (правильной) среды conda - PullRequest
Новая
       43

отправка задания создания змеи в кластер из (правильной) среды conda

1 голос
/ 07 апреля 2020

Я пишу файл snakemake, который будет выполнять несколько операций над несколькими выборками. После того, как я проверил рабочий процесс, запущенный на моем локальном компьютере, я теперь работаю над запуском рабочего процесса в кластере.

Мои первые два правила независимы друг от друга, первое использует fastqc, а другое bwa mem

Эти два правила выглядят так (на данный момент я вызываю рабочий процесс только для одного SAMPLE = 'NIPT-PearlPPlasma-03-PPx_S3downSample'): 

rule fastQC:
    input:
        R1 = FQDIR + "{sample}_R1_001.fastq.gz",
        R2 = FQDIR + "{sample}_R2_001.fastq.gz"
    output:
        directory("fastQC/{sample}")
    conda:
        "envs/NIPTlibPrep.yaml"
    log:
        "logs/fastQC/{sample}.log" # log was giving an error when running at the command line
    shell:
        # 2> {log} at the end of the command removed
        # See wrapper at https:/snakemake-wrappers.readthedocs.io/en/stable/wrappers/fastqc.html
        "mkdir -p fastQC/{wildcards.sample} | fastqc --outdir fastQC/{wildcards.sample} -f fastq {input.R1} {input.R2}"


rule bwa_map: 
    input:
        R1 = FQDIR + "{sample}_R1_001.fastq.gz",
        R2 = FQDIR + "{sample}_R2_001.fastq.gz",
        REF = config['ref']    
    output:
        # wrap output in temp
        "aligned/{sample}.bam"
    log:
        "logs/bwa_mem/{sample}.log" 
    conda:
        "envs/NIPTlibPrep.yaml"
    shell:
        "bwa mem {input.REF} {input.R1} {input.R2} "
        "| samtools view -Sb - > {output} 2> {log}"

Но когда я звоню: 

snakemake -p -s Snakefile_v4_ngs_bngs05b --cluster qsub -j 5 --use-conda

Я получаю: 

Error in rule bwa_map:
    jobid: 10
    output: aligned/NIPT-PearlPPlasma-03-PPx_S3downSample.bam
    log: logs/bwa_mem/NIPT-PearlPPlasma-03-PPx_S3downSample.log (check log file(s) for error message)
    conda-env: /nexusb/nipt/200311_A00154_0454_AHHHKMDRXX/testMetrics/outSnakeMake_test/.snakemake/conda/38107c2c
    shell:
        bwa mem /home/ngs/data/genomes/b37/human_g1k_v37.fasta /nexusb/Novaseq/200311_A00154_0454_AHHHKMDRXX/Unaligned/NIPT-PearlPPlasma-03-PPx_S3downSample_R1_001.fastq.gz /nexusb/Novaseq/200311_A00154_0454_AHHHKMDRXX/Unaligned/NIPT-PearlPPlasma-03-PPx_S3downSample_R2_001.fastq.gz | samtools view -Sb - > aligned/NIPT-PearlPPlasma-03-PPx_S3downSample.bam 2> logs/bwa_mem/NIPT-PearlPPlasma-03-PPx_S3downSample.log
        (one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)
    cluster_jobid: Your job 381368 ("snakejob.bwa_map.10.sh") has been submitted

Error executing rule bwa_map on cluster (jobid: 10, external: Your job 381368 ("snakejob.bwa_map.10.sh") has been submitted, jobscript: /nexusb/nipt/200311_A00154_0454_AHHHKMDRXX/testMetrics/outSnakeMake_test/.snakemake/tmp.bnhr7qck/snakejob.bwa_map.10.sh). For error details see the cluster log and the log files of the involved rule(s).
[Wed Apr  8 17:21:45 2020]
Error in rule fastQC:
    jobid: 1
    output: fastQC/NIPT-PearlPPlasma-03-PPx_S3downSample
    log: logs/fastQC/NIPT-PearlPPlasma-03-PPx_S3downSample.log (check log file(s) for error message)
    conda-env: /nexusb/nipt/200311_A00154_0454_AHHHKMDRXX/testMetrics/outSnakeMake_test/.snakemake/conda/38107c2c
    shell:
        mkdir -p fastQC/NIPT-PearlPPlasma-03-PPx_S3downSample | fastqc --outdir fastQC/NIPT-PearlPPlasma-03-PPx_S3downSample -f fastq /nexusb/Novaseq/200311_A00154_0454_AHHHKMDRXX/Unaligned/NIPT-PearlPPlasma-03-PPx_S3downSample_R1_001.fastq.gz /nexusb/Novaseq/200311_A00154_0454_AHHHKMDRXX/Unaligned/NIPT-PearlPPlasma-03-PPx_S3downSample_R2_001.fastq.gz
        (one of the commands exited with non-zero exit code; note that snakemake uses bash strict mode!)
    cluster_jobid: Your job 381369 ("snakejob.fastQC.1.sh") has been submitted

Error executing rule fastQC on cluster (jobid: 1, external: Your job 381369 ("snakejob.fastQC.1.sh") has been submitted, jobscript: /nexusb/nipt/200311_A00154_0454_AHHHKMDRXX/testMetrics/outSnakeMake_test/.snakemake/tmp.bnhr7qck/snakejob.fastQC.1.sh). For error details see the cluster log and the log files of the involved rule(s).
Shutting down, this might take some time.
Exiting because a job execution failed. Look above for error message

При первом выполнении рабочего процесса я заметил, что среда была создана в .snakemake/conda (относительно Snakefile). Когда я вызываю скрипт во второй раз, не изменяя директивы conda, snakemake использует тот же conda-env.

Описание моей среды выглядит так: 

channels:
  - bioconda
  - conda-forge
dependencies:
  - bwa=0.7.17
  - samtools=1.9
  - picard=2.22.1
  - mosdepth=0.2.6
  - python=3.7.6
  - pandas=1.0.3
  - fastqc=0.11.9

и он сохраняется в envs/NIPTlibPrep.yaml (относительно Snakefile)

Тот факт, что рабочий процесс завершается локально, но не может быть запущен на кластере, я нахожу действительно озадачивающим, особенно учитывая тот факт, что среда с правильными зависимостями была успешно создан.

1 Ответ

0 голосов
/ 07 апреля 2020

Взгляните на: https://snakemake.readthedocs.io/en/stable/snakefiles/deployment.html#integrated -package-management

Вам нужно будет создать среду conda yaml для правила sort_bam, примерно так: 

channels:
  - bioconda
dependencies:
  - samtools

И тогда в вашем правиле вам нужно обратиться к этому файлу под атрибутом conda: 

rule sort_bam:
    input:
        "aligned/{sample}.bam"
    output:
        protected("aligned/{sample}.sorted.bam")
    params:
        THREADS = config['sort_threads']
    conda:
        "samtools.yaml"
    shell:
        "samtools sort -T aligned/{wildcards.sample} "
        "-O bam {input} > {output}"

Затем вы можете позвонить snakemake -np -s Snakefile_v4_ngs_bngs05b --cluster qsub -j 5 --use-conda, и Snakemake позаботится о отдых.

...