retrieve_seqs: function to automatically retrieve sequences from a databank...

Description Usage Arguments Value Author(s) See Also Examples

View source: R/functions.R

Description

this is a function to retrieve sequences of DNA, RNA or protein from databases structured under ACNUC and located on the web, the retrieval is limited to accession numbers as selection criteria, remote access to ACNUC databases works by opening a socket connection on a port (for example on port number 5558 at pbil.univ-lyon1.fr) and by communicating on this socket following the protocol described here

Usage

1
retrieve_seqs(seqnames, acnucdb)

Arguments

seqnames

a character vector specifying the accession numbers of the sequences to be retrieved

acnucdb

a character string specifying the name of the ACNUC database to be searched. Use 'choosebank' without arguments to see a list of available databases

Value

this function returns a named list with the sequences matching the sequence names of class SeqAcnucWeb retrieved by 'query' based on the provided accession numbers the name of each list item is composed of the accession number and the length of the sequence

Author(s)

gerardo esteban antonicelli

See Also

'print_alignment' 'clean_alignment' 'load_alignment' 'make_tree' 'max_parsimony'

Examples

1
2
3
4
5
6
7
8
seqs <- retrieve_seqs(c('P06747', 'P0C569', 'O56773', 'Q5VKP1'), 'swissprot')
## Not run: viruses_annot <- read.csv(system.file('extdata',
                                               'CoV_genomes_annot.csv',
                                               package='firstPackage'))
## End(Not run)
## Not run: accessions <- viruses_annot$accession
## Not run: seqs <- retrieve_seqs(accessions, 'genbank')
## Not run: str(seqs)

geantonicelli/firstPackage documentation built on Aug. 24, 2020, 3:14 a.m.