Add site-specific sequence information.

This function adds a new "site_sequence" column to the variable information table. It is the protein sequence around the glycosylation site. If the head and tail amino acids of the peptide sequence are insufficient, fill with "X".

This function requires the following columns in the variable information tibble:

• "protein": The protein uniprot accession.

• "protein_site": The site on the protein sequence.

Usage

add_site_seq(exp, fasta, n_aa = 7)

Arguments

exp

A glyexp::experiment() object.

fasta

A character string specifying the path to the FASTA file containing protein sequences.

n_aa

The number of amino acids to the left and right of the glycosylation site. For example, if n_aa = 5, the resulting sequence will contain 11 amino acids.

Value

A glyexp::experiment() object with the new "site_sequence" column.