making loose end cluster numbers build from existing

DOH-JDJ0303 · Apr 27, 2024 · d6067d9 · d6067d9
1 parent fa05f11
commit d6067d9
Showing 1 changed file with 6 additions and 0 deletions.
diff --git a/modules/local/bind-clusters.nf b/modules/local/bind-clusters.nf
@@ -17,6 +17,12 @@ process BIND_CLUSTERS {
     prefix = "${taxa}-${segment}"
     """
     header="seq,taxa,segment,cluster"
+    if [ -f "${prefix}-looseends.csv " ]
+    then
+        max_n=\$(cat *main.csv | cut -f 4 -d ',' | sort -rn)
+        mv ${prefix}-looseends.csv tmp && cat tmp | grep -v "\$header" | tr ',' '\t' | awk -v max=\${max_n} -v OFS=',' '{print $1,$2,$3,$4+max}' > ${prefix}-looseends.csv
+        rm tmp
+    fi
     echo \$header > ${prefix}-clusters.csv
     cat ${clusters} | grep -v "\$header" >> ${prefix}-clusters.csv
     """