diff --git a/nextflow.config b/nextflow.config
index 2eb4cb1ec03a5f078d5be9eee4e5a1b2f98cf8f7..1a935aa991270f04f369e2b3d212498d2da076f4 100644
--- a/nextflow.config
+++ b/nextflow.config
@@ -1,6 +1,6 @@
 process.executor='slurm'
 process.queue='128GB,256GB,256GBv1'
 process.clusterOptions = '--hold'
-trace.enabled = true
+trace.enabled = false
 trace.file = 'pipeline_trace.txt'
 trace.field = 'task_id,native_id,process,name,status,exit,submit,start,complete,duration,realtime'
diff --git a/workflow/process_scripts b/workflow/process_scripts
index 0bdf56d5bdb9252ce418fcf7d0c7bdbcee9ff4a1..291e1c352873ce14881ab8b1822fa7d3183b946b 160000
--- a/workflow/process_scripts
+++ b/workflow/process_scripts
@@ -1 +1 @@
-Subproject commit 0bdf56d5bdb9252ce418fcf7d0c7bdbcee9ff4a1
+Subproject commit 291e1c352873ce14881ab8b1822fa7d3183b946b
diff --git a/workflow/scripts/dea.R b/workflow/scripts/dea.R
index 058dd1817f6d448baf3a9c84fdd22f8fc197ad3b..42a57ff6c8986ad375655d7510edf1b1c98da82c 100755
--- a/workflow/scripts/dea.R
+++ b/workflow/scripts/dea.R
@@ -94,6 +94,7 @@ for (i in 1:a) {
       filt.out <- na.omit(output[output$fdr < 0.05,])
       if (nrow(filt.out) > 2) {
       	 subset <- logcpm[row.names(logcpm) %in% filt.out$symbol,]
+	 subset <- subset[!apply(subset, 1, function(x) {any(x == 0)}),]
       	 gnames <- filt.out[c('ensembl','symbol')]
       	 s <- merge(gnames,subset,by.x="ensembl",by.y="row.names",all.x=FALSE,all.y=TRUE,sort=FALSE)
       	 STREE <- hclust(dist(t(subset)))
@@ -148,6 +149,7 @@ for (i in 1:a) {
       filt.out <- na.omit(output[output$fdr < 0.05,])
       if (nrow(filt.out) > 2) {
       subset <- logcpm[row.names(logcpm) %in% filt.out$symbol,]
+      subset <- subset[!apply(subset, 1, function(x) {any(x == 0)}),]
       gnames <- filt.out[c('ensembl','symbol')]
       s <- merge(gnames,subset,by.x="ensembl",by.y="row.names",all.x=FALSE,all.y=TRUE,sort=FALSE)
       STREE <- hclust(dist(t(subset)))