GC {seqinr}R Documentation

Calculates the fractional GC content of nucleic acid sequences

Description

This calculates the fraction of G+C bases of the input nucleic acid sequence(s) It reads in nucleic acid sequences, sums the number of 'G' and 'C' bases and writes out the result as the fraction (in the interval 0.0 to 1.0) of the length of the whole sequence.

Usage

GC(seq)

Arguments

seq any nucleic acid sequence

Value

It returns the fraction of G+C as a numeric vector

Author(s)

D. Charif

References

  To have an overview of the seqinR's functionnality, please consult this vignette: Charif, D., Lobry, J.R. (2005) SeqinR: a contributed package to the R project for statistical computing devoted to biological sequences retrieval and analysis. Springer Verlag, Biological and Medical Physics/Biomedical Series, in preparation.

Examples

   s=s2c("agtctggggggccccttttaagtagatagatagctagtcgta")
   GC(s)

[Package seqinr version 1.0-3 Index]