############# READ BEFORE RUN SCRIPT ###########################

• "Biostrings", "tidyr", "dplyr", "stringr" packages must be installed.

Biostrings package can report a error in installing process, if it is the case run the following code:

if (!require("BiocManager", quietly = TRUE)) install.packages("BiocManager")

BiocManager::install("Biostrings")

for install the other packages run this code:

install.packages("tidyr") install.packages("dplyr") install.packages("stringr")

############# RUNNING THE SCRIPT ###########################

For execute this code, your table must be in csv format "your_tbl.csv".

your_tbl must have the following columns "ID", "sequencep", "sequencen".

sequencep = columns containning the aminoacid sequence with "p" before the modified aminoacid. sequencen = aminoacid sequence.

you have to rename your FASTA file to "human.FASTA", or rename the code in script.

IMPORTANT: you must have "your_tbl", "human/mouse/rat.fasta" in the same folder in R project.