writeFasta | R Documentation |
peptides which do not have protein assignment drop out
writeFasta(file, ...)
file |
where to write |
... |
fasta list or single file |
writes a file.
#example how to create a protein db with decoy sequences
library(seqinr)
#library(prozor)
file = system.file("extdata/fgcz_contaminants2021_20210929.fasta.gz",package = "prozor")
fasta = readPeptideFasta(file = file)
revfasta <- reverseSeq(fasta)
decoyDB <- c(fasta,revfasta)
stopifnot(length(decoyDB) == 2 * length(fasta))
## Not run:
writeFasta(decoyDB, file="test.fasta")
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.