# HG changeset patch
# User bgruening
# Date 1377975605 14400
# Node ID c650ee099c51ce9e1e42439af105188ed8dd7f0d
# Parent 2b3adbe83979b5a541c7bae4ab2d7fb0e07cacee
Uploaded
diff -r 2b3adbe83979 -r c650ee099c51 cmsearch.xml
--- a/cmsearch.xml Sat Aug 31 09:58:04 2013 -0400
+++ b/cmsearch.xml Sat Aug 31 15:00:05 2013 -0400
@@ -5,9 +5,7 @@
infernal
- ##TODO reformat the tabular-output:
- ## -"Each line consists of 18 space-delimited fields followed by a free text target sequence description, as follows:"
- ## remove line that starts with #
+ temp_tabular_output=$(mktemp);
cmsearch
--cpu 12
@@ -21,7 +19,7 @@
$nohmm
$mid
$bitscore_thresholds
- --tblout $outfile
+ --tblout $temp_tabular_output
$g
#if $A:
$A $multiple_alignment_output
@@ -29,13 +27,13 @@
#if $inclusion_thresholds_opts.inclusion_thresholds_selector == "--incE":
--incE $inclusion_thresholds_opts.incE
- #else:
+ #elif $inclusion_thresholds_opts.inclusion_thresholds_selector == "--incT":
--incT $inclusion_thresholds_opts.incT
#end if
#if $reporting_thresholds_opts.reporting_thresholds_selector == "-E":
-E $reporting_thresholds_opts.E
- #else:
+ #elif $reporting_thresholds_opts.reporting_thresholds_selector == "-T":
-T $reporting_thresholds_opts.T
#end if
@@ -46,9 +44,11 @@
#end if
$seqdb
-
+ 2>&1
+ ;
- 2>&1
+ sed -e 's/#.*$//' -e '/^$/d' -e 's/ /\t/g' -e 's/\t/ /18g' $temp_tabular_output > $outfile
+
@@ -117,9 +117,11 @@
-
+
+
+
@@ -140,15 +142,16 @@
-
-
+
+
+
@@ -172,7 +175,6 @@
-