diff --git a/scripts/paris/ranks_tsv_to_h5an.py b/scripts/paris/ranks_tsv_to_h5an.py
new file mode 100644
index 0000000000000000000000000000000000000000..8f45829cb1fbf4f15af2464881943fcc376f841a
--- /dev/null
+++ b/scripts/paris/ranks_tsv_to_h5an.py
@@ -0,0 +1,14 @@
+import anndata as ad
+import signatures
+
+if __name__ == "__main__":
+    import sys
+
+    assert(len(sys.argv) == 4)
+
+    ranks = signatures.load_ranks_csv(sys.argv[1])
+    adata = ad.read(sys.argv[2])
+    adata.layers["ranks"] = ranks
+    adata.write(sys.argv[3], compressed="gzip")
+
+    print(adata, out=sys.stderr, flush=True)
diff --git a/scripts/paris/signatures.py b/scripts/paris/signatures.py
index 16a51da3492c74a68ed747745e0a018a5244e177..8bab33548e4cdc14f8f8930a006ef2f6f286336a 100644
--- a/scripts/paris/signatures.py
+++ b/scripts/paris/signatures.py
@@ -272,6 +272,7 @@ def load_ranks_csv(rankfile):
             print("\r", n, end=" ", file=sys.stderr, flush=True)
             ranks_row = numpy.array([float(r) for r in row[1:]])
             ranks_l.append(ranks_row)
+
     assert(len(ranks_l) > 0)
     assert(len(ranks_l) == len(genes))
     ngenes = len(genes)