Bedtools genomecov adapted to the current module standard

This commit is contained in:
JoseEspinosa 2020-07-15 18:15:32 +02:00
parent 2cbc988d5f
commit 0b123f9e13
6 changed files with 22 additions and 8 deletions

View file

@ -4,4 +4,4 @@ LABEL authors="Jose Espinosa-Carrasco" \
COPY environment.yml /
RUN conda env create -f /environment.yml && conda clean -a
ENV PATH /opt/conda/envs/nf-core-bedtools-complementbed/bin:$PATH
ENV PATH /opt/conda/envs/nf-core-bedtools-genomecov/bin:$PATH

View file

@ -1,18 +1,33 @@
def MODULE = "bedtools_genomecov"
params.publish_dir = MODULE
params.publish_results = "default"
process BEDTOOLS_GENOMECOV {
tag {bam}
publishDir "${params.out_dir}/${params.publish_dir}",
mode: params.publish_dir_mode,
saveAs: { filename ->
if (params.publish_results == "none") null
else filename }
//container "docker.pkg.github.com/nf-core/$MODULE"
container 'quay.io/biocontainers/bedtools:2.29.2--hc088bd4_0'
conda "${moduleDir}/environment.yml"
input:
path (bam)
path (chrom_sizes)
val (bedtools_genomecov_args)
output:
stdout()
path "${bam}.bed", emit: coverage
path "*.version.txt", emit: version
script:
"""
bedtools genomecov -ibam ${bam} -g ${chrom_sizes} ${bedtools_genomecov_args}
bedtools genomecov -ibam ${bam} -g ${chrom_sizes} ${bedtools_genomecov_args} > ${bam}.bed
bedtools --version | sed -n "s/.*\\(v.*\$\\)/\\1/p" > bedtools.version.txt
"""
}

View file

@ -0,0 +1 @@
../../../../../tests/data/bam/JK2067_downsampled_s0.1.bam

View file

@ -1 +0,0 @@
chr1 1780869

View file

@ -0,0 +1 @@
../../../../../tests/data/bed/genome.sizes

View file

@ -2,9 +2,11 @@
nextflow.preview.dsl = 2
params.out_dir = "test_output"
params.fastqc_args = ''
params.publish_dir_mode = "copy"
params.bedtools_genomecov_args = '' //'-bg'
include check_output from '../../../../tests/functions/check_process_outputs.nf' // params(params)
include BEDTOOLS_GENOMECOV from '../main.nf' params(params)
// Define input channels
@ -14,5 +16,4 @@ chrom_sizes = Channel.fromPath('./input_data/genome.sizes')
// Run the workflow
workflow {
BEDTOOLS_GENOMECOV(ch_input, chrom_sizes, params.bedtools_genomecov_args)
// .check_output()
}

View file

@ -1,2 +0,0 @@
docker.enabled = true
params.outdir = './results'