Can anyone give me the code for in-silico digestion of whole genome with two restriction enzyme in "R" through Biostrings package?

Question

[Bhowmik]

Can anyone give me the code for in-silico digestion of whole genome with two restriction enzyme in "R" through Biostrings package?

Answer 1

[Haci]

DigestDNA() from the DECIPHER package does this, from the documentation:

# Restriction digest of Yeast Chr. 1 with EcoRI and EcoRV
data(yeastSEQCHR1)
sites <- RESTRICTION_ENZYMES[c("EcoRI", "EcoRV")]
seqs <- DigestDNA(sites, yeastSEQCHR1)
seqs[[1]]

pos <- DigestDNA(sites, yeastSEQCHR1, type="positions")
str(pos)