Question: Simulate Paired-End Bisulfite Read Data
gravatar for Chris Miller
8.6 years ago by
Chris Miller21k
Washington University in St. Louis, MO
Chris Miller21k wrote:

I'm looking for a tool that can simulate bisulfite paired-end reads.

RMAP includes a some nice tools that allow you to either simulate paired end reads or bisulfite reads, but not both. It's also got parameters for the rate of CpG methylation, the rate of bisulfite conversion, and mutation rate. Does anyone know of something comparable that will generate this data for bisulfite-treated paired-end reads?

Just trying to avoid reinventing the wheel here.

ADD COMMENTlink modified 2.9 years ago by vimal.biochem0 • written 8.6 years ago by Chris Miller21k
gravatar for brentp
8.6 years ago by
Salt Lake City, UT
brentp23k wrote:

Good question. I looked a while ago; it's not a simple task since methylation has a high spatial autocorrelation and depends on gc-content and position relative to genes (not to mention repetitive elements). Maybe using something from the sequence read archive is the way to go here.

Another option is just to use a normal read simulator--your RMAP or dwgsim from dnaa tools and simulate non-bs-treated reads. Certainly, you should expect those reads to map correctly.

From there, if you know something about the patterns methylation in your genome of choice, you could manually convert some C's to T's (and G's to A's) in your reference with a probability based on what you see in the genome.

And if you build such a tool, please share.

ADD COMMENTlink written 8.6 years ago by brentp23k

Yeah, I may dig into the RMAP code, see how they did it, and hack something together here in the next few days. If I do, I'll certainly post back here.

ADD REPLYlink written 8.6 years ago by Chris Miller21k
gravatar for vimal.biochem
2.9 years ago by
vimal.biochem0 wrote:

Hi Chris, As I can see this is an old post, I am not sure if such simulator is available. Although simulators like Sherman BSSim DNemulator etc ( are present but I am looking for a tool with C methylation (CG and CH context) and the rate of bisulfite conversion both. As far as i know above-mentioned tools don't address bisulfite conversion. Please suggest me if I am not aware of something already available.

ADD COMMENTlink written 2.9 years ago by vimal.biochem0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2379 users visited in the last hour