From 02213d0324072ba2322cfc53302461c7cfabc30f Mon Sep 17 00:00:00 2001
From: Martin Aryee <martin.aryee@gmail.com>
Date: Sat, 16 Mar 2024 20:48:38 -0400
Subject: [PATCH] Sort chunks in chunk_pe_fastq.wdl

---
 chunk_pe_fastq.wdl | 33 +++++++++++++++++++++++++++++++--
 1 file changed, 31 insertions(+), 2 deletions(-)

diff --git a/chunk_pe_fastq.wdl b/chunk_pe_fastq.wdl
index 72a8433..e0fbaea 100644
--- a/chunk_pe_fastq.wdl
+++ b/chunk_pe_fastq.wdl
@@ -28,9 +28,17 @@ workflow chunk_pe_fastq {
 		}
 	}
 	
+	call sort_chunks as sort_chunks1 { input:
+		chunks = flatten(chunk_fastq1.chunks)
+	}
+
+	call sort_chunks as sort_chunks2 { input:
+		chunks = flatten(chunk_fastq2.chunks)
+	}
+	
 	output {
-		Array[File] fastq1_chunks = flatten(chunk_fastq1.chunks)
-		Array[File] fastq2_chunks = flatten(chunk_fastq2.chunks)
+		Array[File] fastq1_chunks = sort_chunks1.sorted_chunks
+		Array[File] fastq2_chunks = sort_chunks2.sorted_chunks
 	}
 }
 
@@ -58,3 +66,24 @@ task chunk_fastq_file {
         Array[File] chunks = glob("*.fastq.gz")
     }
 }
+
+task sort_chunks { 
+	input {
+		Array[String] chunks
+	}
+	
+	command {
+		echo ${sep=" " chunks} | tr " " "\n" | sort
+	}
+	
+	runtime {
+        continueOnReturnCode: false
+        docker: "ubuntu"
+        cpu: 1
+    } 
+    
+    output {
+        Array[String] sorted_chunks = read_lines(stdout())
+    }
+    
+}
\ No newline at end of file