seqinr (version 4.2-36)

residuecount: Total number of residues in an ACNUC list

Description

Computes the total number of residues (nucleotides or aminoacids) in all sequences of the list of specified rank.

Usage

residuecount(lrank, socket = autosocket())

Value

A single numeric value corresponding to the total number of residues or NA in case of problem.

Arguments

lrank

the list rank on the ACNUC server

socket

an object of class sockconn connecting to a remote ACNUC database (default is a socket to the last opened database).

Author

J.R. Lobry

References

https://doua.prabi.fr/databases/acnuc.html

citation("seqinr")

See Also

choosebank, query, glr

Examples

if (FALSE) {
  ### Need internet connection
  choosebank("emblTP")
  mylist <- query("mylist", "t=CDS", virtual = TRUE)
  stopifnot(residuecount(glr("mylist")) == 1611439240)
  stopifnot(is.na(residuecount(glr("unknowlist")))) # A warning is issued
  }