Can anyone give me the code for in-silico digestion of whole genome with two restriction enzyme in "R" through Biostrings package?
1 answer
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)
Log in to answer this question.
More posts like this
-
RRBS (Reduced Representation Bisulfite Sequencing)
written by sk_24 •
-
ı cant install clusterprofiler on my Ubuntu 20.04.6 LTS
written by iamsmor •
-
Any online tools to do in silico RE digestion for unmapped nanopore sequencing reads?
written by Yingzi •
-
Genome-wide In Silico Multiple Restriction Digest
written by shakuntala.mitra •
-
Software to in silico test restriction enzymes for GBS sequencing and similar approaches
written by Denis
-
In silico genome comparison tools
written by jingjin2203 •
-
With Python, how can I solve this restriction mapping problem?
written by francois •
-
Digesting contig sequence with multiple enzymes
written by darshan.wash •
-
Restriction Sites Fragments For Gbs/Ddrad
written by Bart Ulaszewski
-
In Silico Genome Digestion ?
written by Bioch'Ti •
