Description

Convert HMMER’s –tblout/–domtblout output to a clean, properly-delimited TSV

Input

name:type
description
pattern

meta

:map

Groovy Map containing sample information e.g. [ id:'sample1' ]

labels

:list

A list of strings, one per entry in files, in the same order. Each becomes the value of the profile column for every row taken from the corresponding file – e.g. the name of the HMM profile that particular tblout/domtblout came from. Not derived from the input filenames, so the module doesn’t depend on any particular producer’s naming convention.

files

:file

One or more HMMER –tblout or –domtblout files (gzipped or not), all of the same kind as given in format. Column names are hardcoded to the hmmsearch/hmmscan/phmmer layout (19 tblout / 23 domtblout columns); other tools with a different column count for the same format, e.g. nhmmer’s 16-column tblout, are rejected with an error rather than silently mislabeled.

*.{tbl,domtbl}{,.gz}

format

:string

Which HMMER output format files contains – either tblout or domtblout

^(tblout|domtblout)$

Output

name:type
description
pattern

tsv

meta

:map

Groovy Map containing sample information e.g. [ id:'sample1' ]

*.tsv.gz

:file

All rows from every input file, concatenated, as a single gzipped TSV with a profile column (from labels) first. Column names and count depend on format: 20 columns (profile + HMMER’s 19) for tblout, 24 (profile + HMMER’s 23) for domtblout.

*.{tblout,domtblout}.tsv.gz

versions_gawk

${task.process}

:string

The name of the process

gawk

:string

The name of the tool

awk -Wversion | sed '1!d; s/.*Awk //; s/,.*//'

:eval

The expression to obtain the version of the tool

Topics

name:type
description
pattern

versions

${task.process}

:string

The name of the process

gawk

:string

The name of the tool

awk -Wversion | sed '1!d; s/.*Awk //; s/,.*//'

:eval

The expression to obtain the version of the tool