kofamscan: TSV output; fix merge_kofamscan_hits.py

metashot · Nov 10, 2022 · 6adfed8 · 6adfed8
1 parent 63cd12c
commit 6adfed8
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 3 deletions.
diff --git a/bin/merge_kofamscan_hits.py b/bin/merge_kofamscan_hits.py
@@ -17,8 +17,7 @@
 for input_fn in args.input:
     input_bn = os.path.basename(input_fn)
     input_id = input_bn.split(args.suffix)[0]
-    input_df = pd.read_fwf(input_fn, skiprows=[1],
-        colspecs = [(0, 1), (2, 22), (22, 29), (54, 1024)])
+    input_df = pd.read_csv(input_fn, delimiter='\t', skiprows=[1])
     input_df = input_df[input_df['#'] == '*']
     input_df = input_df.groupby(["KO", "KO definition"]) \
         [["gene name"]].nunique()
@@ -27,3 +26,10 @@
 
 compl_df.fillna(0, inplace=True)
 compl_df.to_csv('kofamscan_hits.tsv', sep='\t', float_format='%.d')
+Footer
+© 2022 GitHub, Inc.
+Footer navigation
+Terms
+Privacy
+Security
+
diff --git a/modules/kofamscan.nf b/modules/kofamscan.nf
@@ -15,7 +15,7 @@ process kofamscan {
     script:
     """
     exec_annotation \
-        -f detail \
+        -f detail-tsv \
         --cpu ${task.cpus} \
         --tmp-dir tmp \
         -o ${id}.hits.txt \