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