This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Mutate genomes in silico

Hello, I'm looking for a tool that might introduce mutations to genomes. The main idea is to test how mutations will affect mapping (both genome to genome, and short reads to genome). Any recommendations please? Have a great weekend!

genome

You mean to mutate a fasta (genome sequence) or a fastq (reads) file?

The genome sequences, Thanks

3 answers

BamSurgeon for spiking mutations into the reads.

snp-mutator for fasta files.

Try simutator, e.g., every 2kb where a region of length 1.5kb gets 30 snps, 2 insertion and 4 deletions up to length 10bp.

mutate.sh from BBMap suite. Take a look at in-line help for additional options.

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.

Log in to answer this question.