Fixed bug - no output with <1 contig

rotary-genomics · Apr 15, 2022 · fd5acee · fd5acee
1 parent e636236
commit fd5acee
Show file tree

Hide file tree

Showing 2 changed files with 4 additions and 4 deletions.
diff --git a/config.yaml b/config.yaml
@@ -12,7 +12,7 @@ qc_short_r2: "QC_R2.fastq.gz"
 #####################
 
 ## Directory for database installation
-db_dir: "/Data/databases/nanopore-workflows"
+db_dir: "/Data/databases/toyako"
 
 ## Resources
 threads: 40

diff --git a/rules/toyako.smk b/rules/toyako.smk
@@ -659,11 +659,11 @@ checkpoint split_circular_and_linear_contigs:
 
         os.makedirs(output[0], exist_ok=True)
 
-        # Only output files if there is >1 entry
-        if circular_contigs.shape[0] > 1:
+        # Only output files if there is >=1 entry
+        if circular_contigs.shape[0] >= 1:
             circular_contigs['#seq_name'].to_csv(os.path.join(output[0], 'circular.list'), header=None, index=False)
 
-        if linear_contigs.shape[0] > 1:
+        if linear_contigs.shape[0] >= 1:
             linear_contigs['#seq_name'].to_csv(os.path.join(output[0], 'linear.list'), header=None, index=False)