Mercurial > repos > nml > metaspades

changeset 8:f338e9942b27 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
"planemo upload for repository https://github.com/galaxyproject/tools-iuc/tree/master/tools/spades commit 2609367ebbce7e0c3f5b4b2d7c60a82d4bd03246"
author iuc
date Tue, 25 Jan 2022 07:57:04 +0000
parents 7c55feb65057
children e731ea24e48d
files metaspades.xml write_tsv_script.py
diffstat 2 files changed, 18 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/metaspades.xml	Mon Jan 24 18:45:46 2022 +0000
+++ b/metaspades.xml	Tue Jan 25 07:57:04 2022 +0000
@@ -136,12 +136,12 @@
             <output_collection name="out_cr" type="list" count="3">
                 <element name="ecoli_1K.fastq.gz_1.00.0_0.cor">
                     <assert_contents>
-                        <has_size value="130317" delta="1000"/>
+                        <has_size value="130317" delta="5000"/>
                     </assert_contents>
                 </element>
                 <element name="ecoli_1K.fastq.gz_2.00.0_0.cor">
                     <assert_contents>
-                        <has_size value="130317" delta="1000"/>
+                        <has_size value="130317" delta="5000"/>
                     </assert_contents>
                 </element>
                 <element name="ecoli_1K.fastq.gz__unpaired.00.0_0.cor">
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/write_tsv_script.py	Tue Jan 25 07:57:04 2022 +0000
@@ -0,0 +1,16 @@
+#!/usr/bin/env python
+
+import re
+import sys
+
+search_str = r"^>(NODE|\S+)_(\d+)(?:_|\s)length_(\d+)_cov_(\d+\.*\d*).*\$"
+
+replace_str = r"\1_\2\t\3\t\4"
+
+cmd = re.compile(search_str)
+
+sys.stdout.write("#name\tlength\tcoverage\n")
+
+for i, line in enumerate(sys.stdin):
+    if cmd.match(line):
+        sys.stdout.write(cmd.sub(replace_str, line))