From 58bdde43189de11e5ae6e05730e27525d014f1e7 Mon Sep 17 00:00:00 2001
From: "Luis M. Rodriguez-R" <lmrodriguezr@gmail.com>
Date: Thu, 29 Aug 2024 19:58:38 +0200
Subject: [PATCH] Reduce waste by inactivating empty assemblies

---
 lib/miga/version.rb |  2 +-
 scripts/cds.bash    | 10 ++++++++--
 2 files changed, 9 insertions(+), 3 deletions(-)

diff --git a/lib/miga/version.rb b/lib/miga/version.rb
index 3b14318..93d629b 100644
--- a/lib/miga/version.rb
+++ b/lib/miga/version.rb
@@ -12,7 +12,7 @@ module MiGA
   # - String indicating release status:
   #   - rc* release candidate, not released as gem
   #   - [0-9]+ stable release, released as gem
-  VERSION = [1.3, 20, 2].freeze
+  VERSION = [1.3, 20, 3].freeze
 
   ##
   # Nickname for the current major.minor version.
diff --git a/scripts/cds.bash b/scripts/cds.bash
index f726cf0..0d76da2 100755
--- a/scripts/cds.bash
+++ b/scripts/cds.bash
@@ -16,6 +16,14 @@ if [[ -e "../05.assembly/$DATASET.LargeContigs.fna.gz" \
   miga add_result -P "$PROJECT" -D "$DATASET" -r assembly -f
 fi
 
+# Check if the input assembly is empty
+ASM_LEN=$(grep -v '^>' "../05.assembly/${DATASET}.LargeContigs.fna" \
+            | wc -lc | awk '{ print $2-$1 }')
+if [[ "$ASM_LEN" -lt 1 ]] ; then
+  miga edit -P "$PROJECT" -D "$DATASET" --inactivate "Empty assembly"
+  exit 0
+fi
+
 # Run Prodigal
 TYPE=$(miga ls -P "$PROJECT" -D "$DATASET" -m type | cut -f 2)
 case "$TYPE" in
@@ -27,8 +35,6 @@ case "$TYPE" in
     P_LEN=0
     BEST_CT=0
     PROCEDURE=single
-    ASM_LEN=$(grep -v '^>' "../05.assembly/${DATASET}.LargeContigs.fna" \
-      | wc -lc | awk '{ print $2-$1 }')
     [[ "$ASM_LEN" -lt 2000 ]] && PROCEDURE=meta
     echo "# Codon table selection:" > "${DATASET}.ct.t"
     for ct in 11 4 ; do