Skip to content
Snippets Groups Projects
Unverified Commit 6e52362d authored by Nicolas Servant's avatar Nicolas Servant Committed by GitHub
Browse files

Merge pull request #47 from nservant/dev

Fix typo issue in dev
parents a1948117 ef4fb113
No related branches found
No related tags found
No related merge requests found
...@@ -47,7 +47,7 @@ def helpMessage() { ...@@ -47,7 +47,7 @@ def helpMessage() {
Contacts calling Contacts calling
--min_restriction_fragment_size Minimum size of restriction fragments to consider. Default: None --min_restriction_fragment_size Minimum size of restriction fragments to consider. Default: None
--max_restriction_framgnet_size Maximum size of restriction fragmants to consider. Default: None --max_restriction_fragment_size Maximum size of restriction fragments to consider. Default: None
--min_insert_size Minimum insert size of mapped reads to consider. Default: None --min_insert_size Minimum insert size of mapped reads to consider. Default: None
--max_insert_size Maximum insert size of mapped reads to consider. Default: None --max_insert_size Maximum insert size of mapped reads to consider. Default: None
--saveInteractionBAM Save BAM file with interaction tags (dangling-end, self-circle, etc.). Default: False --saveInteractionBAM Save BAM file with interaction tags (dangling-end, self-circle, etc.). Default: False
...@@ -505,10 +505,9 @@ if (!params.dnase){ ...@@ -505,10 +505,9 @@ if (!params.dnase){
set val(oname), file("${prefix}.mapstat") into all_mapstat set val(oname), file("${prefix}.mapstat") into all_mapstat
script: script:
sample = prefix.toString() - ~/(_R1$|_R2$|_val_1$|_val_2$|_1$|_2$)/ sample = prefix.toString() - ~/(_R1|_R2|_val_1|_val_2|_1$|_2)/
tag = prefix.toString() =~/_R1$|_val_1$|_1$/ ? "R1" : "R2" tag = prefix.toString() =~/_R1|_val_1|_1/ ? "R1" : "R2"
oname = prefix.toString() - ~/(\.[0-9]+)$/ oname = prefix.toString() - ~/(\.[0-9]+)$/
""" """
samtools merge -@ ${task.cpus} \\ samtools merge -@ ${task.cpus} \\
-f ${prefix}_bwt2merged.bam \\ -f ${prefix}_bwt2merged.bam \\
...@@ -546,10 +545,9 @@ if (!params.dnase){ ...@@ -546,10 +545,9 @@ if (!params.dnase){
set val(oname), file("${prefix}.mapstat") into all_mapstat set val(oname), file("${prefix}.mapstat") into all_mapstat
script: script:
sample = prefix.toString() - ~/(_R1$|_R2$|_val_1$|_val_2$|_1$|_2$)/ sample = prefix.toString() - ~/(_R1|_R2|_val_1|_val_2|_1|_2)/
tag = prefix.toString() =~/_R1$|_val_1$|_1$/ ? "R1" : "R2" tag = prefix.toString() =~/_R1|_val_1|_1/ ? "R1" : "R2"
oname = prefix.toString() - ~/(\.[0-9]+)$/ oname = prefix.toString() - ~/(\.[0-9]+)$/
""" """
echo "## ${prefix}" > ${prefix}.mapstat echo "## ${prefix}" > ${prefix}.mapstat
echo -n "total_${tag}\t" >> ${prefix}.mapstat echo -n "total_${tag}\t" >> ${prefix}.mapstat
...@@ -563,6 +561,7 @@ if (!params.dnase){ ...@@ -563,6 +561,7 @@ if (!params.dnase){
} }
} }
println(bwt2_merged_bam)
process combine_mapped_files{ process combine_mapped_files{
tag "$sample = $r1_prefix + $r2_prefix" tag "$sample = $r1_prefix + $r2_prefix"
...@@ -721,16 +720,15 @@ process merge_sample { ...@@ -721,16 +720,15 @@ process merge_sample {
output: output:
file("mstats/") into all_mstats file("mstats/") into all_mstats
script: script:
sample = prefix.toString() - ~/(_R1$|_R2$|_val_1$|_val_2$|_1$|_2$)/ sample = prefix.toString() - ~/(_R1|_R2|_val_1|_val_2|_1|_2)/
if ( (fstat =~ /.mapstat/) ){ ext = "mmapstat" } if ( (fstat =~ /.mapstat/) ){ ext = "mmapstat" }
if ( (fstat =~ /.pairstat/) ){ ext = "mpairstat" } if ( (fstat =~ /.pairstat/) ){ ext = "mpairstat" }
if ( (fstat =~ /.RSstat/) ){ ext = "mRSstat" } if ( (fstat =~ /.RSstat/) ){ ext = "mRSstat" }
"""
""" mkdir -p mstats/${sample}
mkdir -p mstats/${sample} merge_statfiles.py -f ${fstat} > mstats/${sample}/${prefix}.${ext}
merge_statfiles.py -f ${fstat} > mstats/${sample}/${prefix}.${ext} """
"""
} }
......
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