@@ -2,8 +2,8 @@ process BISMARK_GENOME_PREPARATION {
22 tag " $genome "
33 label ' process_high'
44
5- conda " bioconda::bismark=0.23.0"
6- container " quay.io/biocontainers/bismark:0.23.0--hdfd78af_1"
5+ // conda "bioconda::bismark=0.23.0"
6+ // container "quay.io/biocontainers/bismark:0.23.0--hdfd78af_1"
77
88 input:
99 path genome
@@ -27,61 +27,6 @@ process BISMARK_GENOME_PREPARATION {
2727 """
2828}
2929
30- // process BISMARK_ALIGN {
31- // tag "$meta.id"
32- // label 'process_high'
33-
34- // conda "bioconda::bismark=0.23.0"
35- // container "quay.io/biocontainers/bismark:0.23.0--hdfd78af_1"
36-
37- // input:
38- // tuple val(meta), path(reads)
39- // path index
40-
41- // output:
42- // tuple val(meta), path("*.bam"), emit: bam
43- // tuple val(meta), path("*_report.txt"), emit: report
44- // path "versions.yml", emit: versions
45-
46- // script:
47- // def args = task.ext.args ?: ''
48- // def prefix = task.ext.prefix ?: "${meta.id}"
49- // """
50- // echo "Copying index directory..."
51- // mkdir -p copied_index
52- // cp -rL $index/* ./copied_index/
53-
54- // echo "Contents of copied index directory:"
55- // ls -lR ./copied_index
56-
57- // echo "Searching for genome file:"
58- // genome_file=\$(find ./copied_index -name "*.fa" -o -name "*.fasta" | head -n 1)
59-
60- // if [ -z "\$genome_file" ]; then
61- // echo "No genome file found in index directory"
62- // exit 1
63- // else
64- // echo "Found genome file: \$genome_file"
65- // echo "First few lines of genome file:"
66- // head -n 5 "\$genome_file"
67- // fi
68-
69- // bismark --genome ./copied_index \\
70- // $args \\
71- // -1 ${reads[0]} \\
72- // -2 ${reads[1]} \\
73- // --basename $prefix
74-
75- // mv ${prefix}_pe.bam ${prefix}.bam
76-
77- // cat <<-END_VERSIONS > versions.yml
78- // "${task.process}":
79- // bismark: \$( bismark --version | sed -e "s/Bismark Version: v//g" )
80- // END_VERSIONS
81- // """
82- // }
83-
84- // In modules/bismark.nf
8530process BISMARK_ALIGN {
8631 tag " $meta . id "
8732 label ' process_bismark_align'
@@ -108,13 +53,6 @@ process BISMARK_ALIGN {
10853 def reads_command = single_end ? " -s $reads " : " -1 ${ reads[0]} -2 ${ reads[1]} "
10954
11055 """
111- echo "Debug: Starting BISMARK_ALIGN for ${ meta.id} "
112- echo "Debug: Single-end: ${ single_end} "
113- echo "Debug: Reads: ${ reads} "
114- echo "Debug: Index: ${ index} "
115- echo "Debug: Prefix: ${ prefix} "
116- echo "Debug: CPUs: ${ task.cpus} "
117- echo "Debug: Memory: ${ task.memory} "
11856
11957 bismark \\
12058 $args \\
@@ -124,9 +62,6 @@ process BISMARK_ALIGN {
12462 -p ${ task.cpus} \\
12563 $reads_command
12664
127- echo "Debug: BISMARK_ALIGN completed for ${ meta.id} "
128-
129-
13065 cat <<-END_VERSIONS > versions.yml
13166 "${ task.process} ":
13267 bismark: \$ (bismark --version | sed -e "s/Bismark Version: v//g")
0 commit comments