diff --git a/snakemodules/analysis.smk b/snakemodules/analysis.smk
index 0806bef26814e2216bef2dffa994692ae4f2852b..babecada1119c31ec7f1f76f6e6ac4fabf0dcfb3 100644
--- a/snakemodules/analysis.smk
+++ b/snakemodules/analysis.smk
@@ -482,7 +482,8 @@ rule coverage_analysis:
         fasta = RESULTS_DIR + "/supercontigs.fa"
     output:
         RESULTS_DIR + "/heatmap_coverage.pdf",
-        RESULTS_DIR + "/heatmap_coverage.txt"
+        RESULTS_DIR + "/heatmap_coverage.txt",
+        RESULTS_DIR + "/contigs_sorted_by_length.txt"
     conda:
         os.path.join(WORKFLOW_PATH,"snakemodules/envs/notebooks.yml")
     benchmark:
diff --git a/snakemodules/notebooks/contig_coverage_heatmap.py.ipynb b/snakemodules/notebooks/contig_coverage_heatmap.py.ipynb
index d40567d9a37f0999210b9e842810f9078f543064..7e637e3a9277bf06e65c47c51cd5206604f8a2c8 100644
--- a/snakemodules/notebooks/contig_coverage_heatmap.py.ipynb
+++ b/snakemodules/notebooks/contig_coverage_heatmap.py.ipynb
@@ -215,7 +215,8 @@
    "metadata": {},
    "outputs": [],
    "source": [
-    "coverage_df_cleaned.to_csv(snakemake.output[1], sep='\\t', index=False)"
+    "coverage_df_cleaned.to_csv(snakemake.output[1], sep='\\t')\n",
+    "df_sort_contig.to_csv(snakemake.output[2],sep='\\t', index=False)"
    ]
   },
   {