pfamscan: Query the online service PfamScan

Description Usage Arguments Details References Examples

View source: R/pfamscan.R

Description

This uses the EBI API. Details here: https://www.ebi.ac.uk/Tools/common/tools/help/index.html

Usage

1
2
3
4
5
6
7
8
9
pfamscan_retrieve(query_id, ..., api_url = "https://www.ebi.ac.uk/",
  path = "Tools/services/rest/pfamscan", errorfun = message,
  verb = TRUE, json_args = list())

pfamscan(fasta_str, email, ..., evalue = 10, asp = "false",
  api_url = "https://www.ebi.ac.uk/",
  path = "Tools/services/rest/pfamscan", maxchecktime = 120,
  wait = 1, errorfun = message, query_size = 100L, verb = TRUE,
  json_args = list())

Arguments

query_id

Character scalar. Id of the query to retrieve, e.g., "pfamscan-R20191015-22133-0197-415408-p2m"

...

Further arguments passed to httr::POST() and httr::GET().

api_url, path

URL path to the API.

errorfun

A function to either print error or message.

verb

Logical scalar. When TRUE the function will print informative messages.

json_args

Arguments passed to jsonlite::fromJSON() (see details).

fasta_str

A string representing a single sequence in fasta format (details)

email, evalue, asp

Query parameters, see here.

maxchecktime

Max wait time for the results (seconds).

wait

Time, in seconds, to wait between checks.

query_size

Maximum number of sequences per query (limit 100).

Details

The function pfamscan_retrieve can be used to retrieve the results of an specific query. This is useful if, for example, the query was correctly posted but the function ran out of time waiting for it to be completed. You can check the list of posted queries during the session using pfamscanr_queries

If fasta_str is a vector including multiple sequences to submit, the function will make either a single call to the API or multiple depending on query_size. As of now, the API limits 100 sequences per query.

The function returns the response (in JSON format) parsed using the function jsonlite::fromJSON().

References

Madeira F, Park YM, Lee J, et al. The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Research. 2019 Jul;47(W1):W636-W641. doi: 10.1093/nar/gkz268.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
## Not run: 
# This is an example extracted directly from the website
# https://www.ebi.ac.uk/Tools/pfa/pfamscan/
ans <- pfamscan(
  fasta_str = ">sp|P35858|ALS_HUMAN Insulin-like growth factor-binding protein complex acid labile subunit OS=Homo sapiens GN=IGFALS PE=1 SV=1
MALRKGGLALALLLLSWVALGPRSLEGADPGTPGEAEGPACPAACVCSYDD
DADELSVFCSSRNLTRLPDGVPGGTQALWLDGNNLSSVPPAAFQNLSSLGF
LNLQGGQLGSLEPQALLGLENLCHLHLERNQLRSLALGTFAHTPALASLGL
SNNRLSRLEDGLFEGLGSLWDLNLGWNSLAVLPDAAFRGLGSLRELVLAGN
RLAYLQPALFSGLAELRELDLSRNALRAIKANVFVQLPRLQKLYLDRNLIA
AVAPGAFLGLKALRWLDLSHNRVAGLLEDTFPGLLGLRVLRLSHNAIASLR
PRTFKDLHFLEELQLGHNRIRQLAERSFEGLGQLEVLTLDHNQLQEVKAGA
FLGLTNVAVMNLSGNCLRNLPEQVFRGLGKLHSLHLEGSCLGRIRPHTFTG
LSGLRRLFLKDNGLVGIEEQSLWGLAELLELDLTSNQLTHLPHRLFQGLGK
LEYLLLSRNRLAELPADALGPLQRAFWLDVSHNRLEALPNSLLAPLGRLRY
LSLRNNSLRTFTPQPPGLERLWLEGNPWDCGCPLKALRDFALQNPSAVPRF
VQAICEGDDCQPPAYTYNNITCASPPEVVGLDLRDLSEAHFAPC",
  email     = "your@email.com",
  httr::config(connecttimeout=60)
)

## End(Not run)

USCbiostats/pfamscanr documentation built on Dec. 31, 2020, 5:37 p.m.