diff --git a/main.nf b/main.nf index 8f6fffbc5760af09831c5d2f74247f2d198e55c6..f7c977b0e69a6f394d34b7408d149ca7ee58da4e 100644 --- a/main.nf +++ b/main.nf @@ -978,10 +978,10 @@ process compartment_calling { script: """ - #cooltools genome binnify ${chrsize} ${res} > genome_bins.txt - #cooltools genome gc genome_bins.txt ${fasta} > genome_gc.txt + cooltools genome binnify ${chrsize} ${res} > genome_bins.txt + cooltools genome gc genome_bins.txt ${fasta} > genome_gc.txt cooltools call-compartments --contact-type cis -o ${sample}_compartments ${cool} - #awk -F"\t" 'NR>1{OFS="\t"; if(\$6==""){\$6=0}; print \$1,\$2,\$3,\$6}' ${sample}_compartments.cis.vecs.tsv | sort -k1,1 -k2,2n > ${sample}_compartments.cis.E1.bedgraph + awk -F"\t" 'NR>1{OFS="\t"; if(\$6==""){\$6=0}; print \$1,\$2,\$3,\$6}' ${sample}_compartments.cis.vecs.tsv | sort -k1,1 -k2,2n > ${sample}_compartments.cis.E1.bedgraph """ }