add filtering before reformating cluster

genepi · Jul 31, 2024 · 4226acc · 4226acc
1 parent 8392a21
commit 4226acc
Showing 1 changed file with 8 additions and 1 deletion.
diff --git a/lib/workflows/umi-pipeline.nf b/lib/workflows/umi-pipeline.nf
@@ -55,7 +55,14 @@ workflow UMI_PIPELINE {
 
     main:
         CLUSTER( detected_umis_ch, raw )
-        REFORMAT_FILTER_CLUSTER( CLUSTER.out.cluster_fastas, raw, umi_parse_clusters )
+
+        CLUSTER.out.cluster_fastas
+        .transpose( by: 1 )
+        .filter{ sample, fasta -> fasta.countFasta() >= params.min_reads_per_cluster }
+        .groupTuple()
+        .set{ filtered_clusters }
+
+        REFORMAT_FILTER_CLUSTER( filtered_clusters, raw, umi_parse_clusters )
 }
 
 //////////////////