Question: How to use BBmap (mutate.sh) to simulate the human genome?
0
gravatar for kamanovae
3 months ago by
kamanovae0
kamanovae0 wrote:

Hi! I want to simulate the human genome. I found a suitable BBMap mutate.sh) program. The mutate.sh only outputs reads that contain mutations. But I need to get a fast output file that would contain reads with mutation and without mutations. I want to try to maintain the coverage of the reference genome. How can i do this?

Now I use the command to run:

1.bbmap/mutate.sh in=reference/gh19.fasta out = reference/hg19_with_mut.fasta id = 0.99 prefix = bbmap

I can use the program BBMap randomreads.sh) , but it looks much more difficult to run and I'm afraid to unconsciously get an undesirable result

bbmap simulation human genome • 139 views
ADD COMMENTlink modified 11 weeks ago by Biostar ♦♦ 20 • written 3 months ago by kamanovae0
1

One does not really simulate a genome but you simulate reads using the reference for that genome. randomreads.sh is not difficult to run.

Do you need to simulate a genome at this point? There are plenty available in databases. You can run mutate.sh on one of them to introduce the mutations, if you don't want to simulate a new dataset.

ADD REPLYlink modified 3 months ago • written 3 months ago by genomax87k

At this stage, I need to introduce mutations in the reference genome, and then I plan to use the NanoSim program to simulate nanopore reads. My final goal is nanopore reads with introduced mutations. Which program is better for the first step?

ADD REPLYlink written 3 months ago by kamanovae0

Then use mutate.sh.

Written by Brian Bushnell
Last modified August 6, 2019

Description:  Creates a mutant version of a genome.

Usage:  mutate.sh in=<input file> out=<output file> id=<identity>

I/O parameters:
in=<file>       Input genome.
out=<file>      Output mutant genome.
vcf=<file>      Output VCF file showing variations added.

Take a look at processing parameters to control the mutations.

ADD REPLYlink written 3 months ago by genomax87k
Please log in to add an answer.

Help
Access

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