peptides which do not have protein assignment drop out

writeFasta(file, ...)

Arguments

file

where to write

...

fasta list or single file

Value

writes a file.

Examples

#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))
if (FALSE) {
writeFasta(decoyDB, file="test.fasta")
}