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

Hi! I want to simulate the human genome. I found a suitable BBMap program. The 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/ in=reference/gh19.fasta out = reference/hg19_with_mut.fasta id = 0.99 prefix = bbmap

I can use the program BBMap , 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

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

Do you need to simulate a genome at this point? There are plenty available in databases. You can run 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

Written by Brian Bushnell
Last modified August 6, 2019

Description:  Creates a mutant version of a genome.

Usage: 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.


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