Skip to content
Snippets Groups Projects
Unverified Commit 7619538a authored by Laurent Modolo's avatar Laurent Modolo
Browse files

SNP_calling.nf: fix collect problem

parent c9fc9da5
No related branches found
No related tags found
No related merge requests found
...@@ -139,19 +139,18 @@ sambamba sort -t ${task.cpus} --tmpdir=./tmp -o ${file_id}_sorted.bam ${bam} ...@@ -139,19 +139,18 @@ sambamba sort -t ${task.cpus} --tmpdir=./tmp -o ${file_id}_sorted.bam ${bam}
sorted_bam_files.into { sorted_bam_files.into {
sorted_bam_file_norm; sorted_bam_file_norm;
sorted_bam_file_tumor sorted_bam_file_tumor;
} }
collect_sorted_bam_file_norm = sorted_bam_file_norm collect_sorted_bam_file_norm = sorted_bam_file_norm
.filter{ normal_sample.contains(it[0]) } .filter{ normal_sample.contains(it[0]) }
.map { it -> it[1]} .map { it -> it[1]}
.collect() .buffer( size: normal_sample.size())
.map { it -> ["normal_sample", it]} .map { it -> ["normal_sample", it]}
collect_sorted_bam_file_tumor = sorted_bam_file_tumor collect_sorted_bam_file_tumor = sorted_bam_file_tumor
.filter{ tumor_sample.contains(it[0]) } .filter{ tumor_sample.contains(it[0]) }
.map { it -> it[1]} .map { it -> it[1]}
.collect() .buffer( size: tumor_sample.size())
.map { it -> ["tumor_sample", it]} .map { it -> ["tumor_sample", it]}
collect_sorted_bam_file = Channel.create() collect_sorted_bam_file = Channel.create()
...@@ -246,7 +245,7 @@ samtools faidx ${fasta} ...@@ -246,7 +245,7 @@ samtools faidx ${fasta}
haplotypecaller_bam_files_norm = haplo_bam_files_norm haplotypecaller_bam_files_norm = haplo_bam_files_norm
.filter{ "normal_sample" == it[0] } .filter{ "normal_sample" == it[0] }
haplotypecaller_bam_files_tumor = haplo_bam_files_tumor haplotypecaller_bam_files_tumor = haplo_bam_files_tumor
.filter{ "tumor_sample" == it[0] } .filter{ "tumor_sample" == it[0] }
indexed_bam_files.into { indexed_bam_files.into {
index_bam_files_norm; index_bam_files_norm;
...@@ -257,6 +256,7 @@ indexed_bam_files_norm = index_bam_files_norm ...@@ -257,6 +256,7 @@ indexed_bam_files_norm = index_bam_files_norm
indexed_bam_files_tumor = index_bam_files_tumor indexed_bam_files_tumor = index_bam_files_tumor
.filter{ "tumor_sample" == it[0] } .filter{ "tumor_sample" == it[0] }
/*
process HaplotypeCaller { process HaplotypeCaller {
tag "$file_id" tag "$file_id"
cpus 4 cpus 4
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment