Mercurial > repos > kkonganti > cfsan_cronology
diff 0.1.0/modules/custom/bwa_idx_mem/main.nf @ 0:c8597e9e1a97
"planemo upload"
author | kkonganti |
---|---|
date | Mon, 27 Nov 2023 12:37:44 -0500 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/0.1.0/modules/custom/bwa_idx_mem/main.nf Mon Nov 27 12:37:44 2023 -0500 @@ -0,0 +1,53 @@ +process BWA_IDX_MEM { + tag "$meta.id" + label 'process_micro' + + module (params.enable_module ? "${params.swmodulepath}${params.fs}bwa${params.fs}0.7.17" : null) + conda (params.enable_conda ? "bioconda::bwa=0.7.17 conda-forge::perl" : null) + container "${ workflow.containerEngine == 'singularity' && !task.ext.singularity_pull_docker_container ? + 'https://depot.galaxyproject.org/singularity/bwa:0.7.17--he4a0461_11' : + 'quay.io/biocontainers/bwa:0.7.17--he4a0461_11' }" + + input: + tuple val(meta), path(genome), path(reads) + + output: + tuple val(meta), path("*.sam"), emit: aligned_sam + path "versions.yml" , emit: versions + + when: + + + script: + def args = task.ext.args ?: '' + def args2 = task.ext.args2 ?: '' + def prefix = task.ext.prefix ?: "${meta.id}" + """ + + bwa index $args $genome + if [ "${params.fq_single_end}" = "false" ]; then + bwa mem \\ + $args2 \\ + -t $task.cpus \\ + -a $genome \\ + ${reads[0]} > ${prefix}.aligned_1.sam + bwa mem \\ + $args2 \\ + -t $task.cpus \\ + -a $genome \\ + ${reads[1]} > ${prefix}.aligned_2.sam + else + bwa mem \\ + $args2 \\ + -t $task.cpus \\ + -a $genome \\ + $reads > ${prefix}.aligned.sam + + fi + + cat <<-END_VERSIONS > versions.yml + "${task.process}": + bwa: \$(echo \$(bwa 2>&1) | sed 's/^.*Version: //; s/Contact:.*\$//') + END_VERSIONS + """ +} \ No newline at end of file