Skip to content

An R client for EMBL-EBI's PfamScan API

License

Unknown, MIT licenses found

Licenses found

Unknown
LICENSE
MIT
LICENSE.md
Notifications You must be signed in to change notification settings

USCbiostats/pfamscanr

Repository files navigation

pfamscanr

Travis build status AppVeyor build status

The goal of pfamscanr is to …

Installation

You can download the development version from GitHub with:

# install.packages("devtools")
devtools::install_github("USCbiostats/pfamscanr")

Example

This is an example extracted directly from the website https://www.ebi.ac.uk/Tools/pfa/pfamscan/

library(pfamscanr)
ans <- pfamscan(
  fasta_str = ">sp|P35858|ALS_HUMAN OS=Homo sapiens GN=IGFALS PE=1 SV=1
  MALRKGGLALALLLLSWVALGPRSLEGADPGTPGEAEGPACPAACVCSYDDDADEL
  SVFCSSRNLTRLPDGVPGGTQALWLDGNNLSSVPPAAFQNLSSLGFLNLQGGQLGS
  LEPQALLGLENLCHLHLERNQLRSLALGTFAHTPALASLGLSNNRLSRLEDGLFEG
  LGSLWDLNLGWNSLAVLPDAAFRGLGSLRELVLAGNRLAYLQPALFSGLAELRELD
  LSRNALRAIKANVFVQLPRLQKLYLDRNLIAAVAPGAFLGLKALRWLDLSHNRVAG
  LLEDTFPGLLGLRVLRLSHNAIASLRPRTFKDLHFLEELQLGHNRIRQLAERSFEG
  LGQLEVLTLDHNQLQEVKAGAFLGLTNVAVMNLSGNCLRNLPEQVFRGLGKLHSLH
  LEGSCLGRIRPHTFTGLSGLRRLFLKDNGLVGIEEQSLWGLAELLELDLTSNQLTH
  LPHRLFQGLGKLEYLLLSRNRLAELPADALGPLQRAFWLDVSHNRLEALPNSLLAP
  LGRLRYLSLRNNSLRTFTPQPPGLERLWLEGNPWDCGCPLKALRDFALQNPSAVPR
  FVQAICEGDDCQPPAYTYNNITCASPPEVVGLDLRDLSEAHFAPC",
  email = "your@email.com",
  httr::config(connecttimeout=60)
)
ans
#>   model_length
#> 1           28
#> 2           61
#> 3           61
#> 4           61
#> 5           61
#> 6           61
#>                                                                                                                                                                                                                                                                                                    align
#> 1                                                                                                             #HMM       aCPrpChCs......alvVnCsergLtavPrdlP, #MATCH     aCP+ C+Cs      +l+V+Cs+r+Lt++P+++P, #PP        7*******************************99, #SEQ       ACPAACVCSydddadELSVFCSSRNLTRLPDGVP
#> 2                 #HMM       sLdLsnNrltslddeaFkglsnLkvLdLsnNllttlspgafsgLpsLrsLdLsgNrL, #MATCH      L L+ N+l+s+++ aF++ls+L  L+L   +l +l+p+a+ gL++L +L+L++N+L, #PP        6789***************************************************98, #SEQ       ALWLDGNNLSSVPPAAFQNLSSLGFLNLQGGQLGSLEPQALLGLENLCHLHLERNQL
#> 3 #HMM       pnLtsLdLsnNrltslddeaFkglsnLkvLdLsnNllttlspgafsgLpsLrsLdLsgNrL, #MATCH     p L sL LsnNrl+ l+d+ F+gl +L  L+L+ N+l +l   af+gL sLr+L L gNrL, #PP        6799*******************************************************98, #SEQ       PALASLGLSNNRLSRLEDGLFEGLGSLWDLNLGWNSLAVLPDAAFRGLGSLRELVLAGNRL
#> 4         #HMM       nLtsLdLsnNrltslddeaFkglsnLkvLdLsnNllttlspgafsgLpsLrsLdLsgNr, #MATCH     +L++LdLs+N l +++  +F +l++L+ L+L +Nl+  + pgaf gL++Lr+LdLs+Nr, #PP        59********************************************************8, #SEQ       ELRELDLSRNALRAIKANVFVQLPRLQKLYLDRNLIAAVAPGAFLGLKALRWLDLSHNR
#> 5         #HMM       LtsLdLsnNrltslddeaFkglsnLkvLdLsnNllttlspgafsgLpsLrsLdLsgNrL, #MATCH     L++L Ls+N + sl++ +Fk+l+ L++L+L++N+++ l +++f+gL +L+ L+L++N+L, #PP        99*******************************************************98, #SEQ       LRVLRLSHNAIASLRPRTFKDLHFLEELQLGHNRIRQLAERSFEGLGQLEVLTLDHNQL
#> 6     #HMM       nLtsLdLsnNrltslddeaFkglsnLkvLdLsnNllttlspgafsgLpsLrsLdLsgNrL, #MATCH     +L+sL+L  + l  +++++F+gls+L+ L L++N l  +++++++gL+ L +LdL++N+L, #PP        69********************************************************98, #SEQ       KLHSLHLEGSCLGRIRPHTFTGLSGLRRLFLKDNGLVGIEEQSLWGLAELLELDLTSNQL
#>   env.to env.from  name        acc sig  evalue
#> 1     73       40 LRRNT PF01462.18   1 1.5e-09
#> 2    134       74 LRR_8  PF13855.6   1 4.6e-10
#> 3    206      146 LRR_8  PF13855.6   1 5.1e-14
#> 4    278      219 LRR_8  PF13855.6   1 1.5e-16
#> 5    350      290 LRR_8  PF13855.6   1 4.2e-16
#> 6    446      386 LRR_8  PF13855.6   1   5e-13
#>                                    desc hmm.to hmm.from act_site   type bits
#> 1 Leucine rich repeat N-terminal domain     28        1       NA Family 37.3
#> 2                   Leucine rich repeat     61        5       NA Repeat 39.1
#> 3                   Leucine rich repeat     61        1       NA Repeat 51.8
#> 4                   Leucine rich repeat     60        2       NA Repeat 59.9
#> 5                   Leucine rich repeat     61        3       NA Repeat 58.4
#> 6                   Leucine rich repeat     61        2       NA Repeat 48.6
#>      clan seq.to seq.from            seq.name
#> 1 No_clan     73       40 sp|P35858|ALS_HUMAN
#> 2  CL0022    134       78 sp|P35858|ALS_HUMAN
#> 3  CL0022    206      146 sp|P35858|ALS_HUMAN
#> 4  CL0022    277      219 sp|P35858|ALS_HUMAN
#> 5  CL0022    350      292 sp|P35858|ALS_HUMAN
#> 6  CL0022    446      387 sp|P35858|ALS_HUMAN

About

An R client for EMBL-EBI's PfamScan API

Topics

Resources

License

Unknown, MIT licenses found

Licenses found

Unknown
LICENSE
MIT
LICENSE.md

Code of conduct

Stars

Watchers

Forks

Releases

No releases published

Packages

No packages published

Languages