Module genomecov added

This commit is contained in:
JoseEspinosa 2020-07-14 16:00:10 +02:00
parent f6894c7d3a
commit 9d79d64eb9
8 changed files with 85 additions and 0 deletions

View file

@ -0,0 +1,7 @@
FROM nfcore/base
LABEL authors="Jose Espinosa-Carrasco" \
description="Docker image containing all requirements for nf-core/modules/bedtools/genomecov"
COPY environment.yml /
RUN conda env create -f /environment.yml && conda clean -a
ENV PATH /opt/conda/envs/nf-core-bedtools-complementbed/bin:$PATH

View file

@ -0,0 +1,9 @@
# You can use this file to create a conda environment for this pipeline:
# conda env create -f environment.yml
name: nf-core-bedtools-genomecov
channels:
- conda-forge
- bioconda
- defaults
dependencies:
- bioconda::bedtools=2.29.2

View file

@ -0,0 +1,18 @@
process GENOMECOV {
tag {bam}
container 'quay.io/biocontainers/bedtools:2.29.2--hc088bd4_0'
input:
path (bam)
path (fasta_sizes)
val (genomecov_args)
output:
stdout()
script:
"""
bedtools genomecov -ibam ${bam} -g ${fasta_sizes} ${genomecov_args}
"""
}

View file

@ -0,0 +1,30 @@
name: bedtools genomecov
description: Returns feature coverage for a given genome in different formats
keywords:
- complement
tools:
- bedtools:
description: |
Bedtools is a software package that provides with a toolset to perform genome arithmetic operations.
homepage: https://bedtools.readthedocs.io/en/latest/index.html
documentation: https://bedtools.readthedocs.io/en/latest/index.html
doi: 10.093/bioinformatics/btq033
input:
-
- input_file:
type: file
description: Input genomic coordinates file
pattern: "*.{bam}"
- fasta_sizes:
type: file
description: Genome chromosome sizes
pattern: "*.{txt,sizes}"
output:
-
- index:
type: stdout,file
description:
pattern: "stdout,*.{bed,bedGraph}"
authors:
- "@JoseEspinosa"

View file

@ -0,0 +1 @@
chr1 1780869

View file

@ -0,0 +1,18 @@
#!/usr/bin/env nextflow
nextflow.preview.dsl = 2
params.genomecov_args = '' //'-bg'
include check_output from '../../../../tests/functions/check_process_outputs.nf' // params(params)
include GENOMECOV from '../main.nf' params(params)
// Define input channels
ch_input = Channel.fromPath('./input_data/JK2067_downsampled_s0.1.bam')
chrom_sizes = Channel.fromPath('./input_data/genome.sizes')
// Run the workflow
workflow {
GENOMECOV(ch_input, chrom_sizes, params.genomecov_args)
// .check_output()
}

View file

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