# 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 @@ -