Skip to contents

Read the headers from a profile HMM into a list.

Usage

multi_read_hmm_header(directory)

Arguments

directory

Path to a directory containing profile HMMs (*.hmm)

Value

A base::list().

Examples

hmms <- multi_read_hmm_header(padlocdev_example("padloc-db/hmm"))
hmms[1]
#> $PDLC00150
#> $PDLC00150$header
#> $PDLC00150$header$`HMMER3/f`
#> [1] "[3.3 | Nov 2019]"
#> 
#> $PDLC00150$header$NAME
#> [1] "HamA1_00001"
#> 
#> $PDLC00150$header$LENG
#> [1] "310"
#> 
#> $PDLC00150$header$ALPH
#> [1] "amino"
#> 
#> $PDLC00150$header$ACC
#> [1] "PDLC00150"
#> 
#> $PDLC00150$header$DESC
#> [1] "DUF1837"
#> 
#> $PDLC00150$header$MAXL
#> integer(0)
#> 
#> $PDLC00150$header$RF
#> [1] "no"
#> 
#> $PDLC00150$header$MM
#> [1] "no"
#> 
#> $PDLC00150$header$CONS
#> [1] "yes"
#> 
#> $PDLC00150$header$CS
#> [1] "no"
#> 
#> $PDLC00150$header$MAP
#> [1] "yes"
#> 
#> $PDLC00150$header$DATE
#> [1] "Thu Jan  7 16:34:27 2021"
#> 
#> $PDLC00150$header$COM
#> integer(0)
#> 
#> $PDLC00150$header$NSEQ
#> [1] "49"
#> 
#> $PDLC00150$header$EFFN
#> [1] "1.187317"
#> 
#> $PDLC00150$header$CKSUM
#> [1] "1833952886"
#> 
#> $PDLC00150$header$GA
#> integer(0)
#> 
#> $PDLC00150$header$TC
#> integer(0)
#> 
#> $PDLC00150$header$NC
#> integer(0)
#> 
#> $PDLC00150$header$`STATS LOCAL MSV`
#> [1] "-10.8865  0.70101"
#> 
#> $PDLC00150$header$`STATS LOCAL VITERBI`
#> [1] "-11.9205  0.70101"
#> 
#> $PDLC00150$header$`STATS LOCAL FORWARD`
#> [1] "-5.4252  0.70101"
#> 
#> 
#>