Skip to content
Snippets Groups Projects
Commit a744be62 authored by nservant's avatar nservant
Browse files

[MODIF] add gc content for compartments

parent a6faf8c7
No related branches found
No related tags found
No related merge requests found
...@@ -978,10 +978,10 @@ process compartment_calling { ...@@ -978,10 +978,10 @@ process compartment_calling {
script: script:
""" """
#cooltools genome binnify ${chrsize} ${res} > genome_bins.txt cooltools genome binnify ${chrsize} ${res} > genome_bins.txt
#cooltools genome gc genome_bins.txt ${fasta} > genome_gc.txt cooltools genome gc genome_bins.txt ${fasta} > genome_gc.txt
cooltools call-compartments --contact-type cis -o ${sample}_compartments ${cool} 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
""" """
} }
......
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