coil (version 1.0.0)

coi5p_pipe: Run the entire coi5p pipeline for an input sequence.

Description

This function will take a raw DNA sequence string and run each of the coi5p methods in turn (coi5p, frame, translate, indel_check). Note that if you are not interested in all components of the output (i.e. only want sequences set in frame reading or translated), then the coi5p analysis functions can be called individually to avoid unnecessary computation.

Usage

coi5p_pipe(x, ..., name = character(), trans_table = 0,
  frame_offset = 0, indel_threshold = -363.87)

Arguments

x

a nucleotide string. Valid characters within the nucleotide string are: a,t,g,c,-,n. The nucleotide string can be input as upper case, but will be automatically converted to lower case.

...

additional arguments to be passed between methods.

name

an optional character string. Identifier for the sequence.

trans_table

The translation table to use for translating from nucleotides to amino acids. Default is 0, which indicates that censored translation should be performed. If the taxonomy of the sample is known, use the function which_trans_table() to determine the translation table to use.

frame_offset

The offset to the reading frame to be applied for translation. By default the offset is zero, so the first character in the framed sequence is considered the first nucleotide of the first codon. Passing frame_offset = 1 would make the second character in the framed sequence the the first nucleotide of the first codon.

indel_threshold

the log likelihood threshold used to assess whether or not sequences are likely to contain an indel. Default is -345.95. Values lower than this will be classified as likely to contain an indel and values higher will be classified as not likely to contain an indel.

Value

an object of class "coi5p"

See Also

coi5p

frame

translate

indel_check

which_trans_table

Examples

# NOT RUN {
dat = coi5p_pipe(example_nt_string)
#full coi5p object can then be printed
dat
#components of output coi5p object can be called individually:
dat$raw    #raw input sequence
dat$name   #name that was passed
dat$framed #sequence in common reading frame
dat$aaSeq  #sequence translated to amino acids (censored)
dat$indel_likely #whether an insertion or deletion likely exists in the sequence
dat$stop_codons #whether or not there are stop codons in the amino acid sequence

dat = coi5p_pipe(example_nt_string , trans_table = 5)
dat$aaSeq #sequence translated to amino acids using designated translation table
# }