NGS reads simulation
Entering edit mode
6.9 years ago
wangyi2412 ▴ 230

Hi, everyone! 

      I am studying the performance of my algorithm, where I need simulation. I looked up the ones used by 1000genomes, but people there said it was outdated, and suggested finding a new one.

      I overheard some software called ART, but cannot find it on web. 

      I also read some similar paper which used simulation, but all of them did not point out what existing software or data sets they used.

      Beside the software or datasets to do the simulations in box, I want to know more about the details of and the principle behind the simulation.

       The settings are like this:

1 First construct the diploid of a human(only consider SNPs/indels, not including other type of variations)

2 generate templates with Gussian distributed length and coming with equal prob from the 4 strand of DNA(+/- strand of two homologous chromosomes, with the error rate similar with that of the sequencing machine like illumina hiseq 2000

3 get 100 bp reads from each template.

         The key is how to construct the diploid of a human so that it best resemble a "typical" person in a population in study. Anyone has any idea? Randomly select of a bp to be different from the ref with the prob. of the mutation rate, say 1%? But the mutation rate should be different on different regions, so how to simulate this scenario? Or to the aim of the study, as long as the simulation is not for study depending on the distribution of the variations, this could be omitted

       Thank you very much!






sequencing simulation • 9.5k views
Entering edit mode
6.9 years ago
Felix Francis ▴ 570

Here is a list of genetic simulation resources. I believe there are several ones that would suit your needs.



Software Resource Brief Description and Homepage
Aladyn Tools to investigate how demographic parameters, populations genetics and abiotic conditions affect the rate of adaptation
ALF A Simulation Framework for Genome Evolution
ART ART is a set of simulation tools to generate synthetic next-generation sequencing data by mimicking real sequencing process with empirical error models or quality profiles.
Bayesian Serial SimCoal Bayesian Serial SimCoal, (BayeSSC) is a modification of SIMCOAL 1.0, a program written by Laurent Excoffier, John Novembre, and Stefan Schneider.
BaySICS An integral platform with a graphical interface for statistical inference based on approximate Bayesian computation.
BEERS BEERS was designed to benchmark RNA-Seq alignment algorithms and also algorithms that aim to reconstruct different isoforms and alternate splicing from RNA-Seq data
BOTTLENECK Bottleneck is a program for detecting recent effective population size reductions from allele data frequencies
BottleSim BottleSim is a computer simulation program for simulating the process of population bottlenecks
CASS Protein Sequence Simulation
CDPOP CDPOP is a landscape genetics tool for simulating the emergence of spatial genetic structure in populations resulting from specified landscape processes governing organism movement behavior.
Classical Genetics Simulator Web-based simulation software
CoaSim CoaSim is a tool for simulating the coalescent process with recombination and geneconversion under various demographic models.
cosi The cosi package is written in C and is available as a tar file.
CS-PSeq-Gen A program to simulate the evolution of protein sequences under the constraints of the information of a particular reconstructed phylogeny
DAWG An application designed to simulate the evolution of recombinant DNA sequences in continuous time
Easypop EASYPOP is an individual based model intended to simulate datasets under a very broad range of conditions
EggLib EggLib is a C++/Python library and program package for evolutionary genetics and genomics.
EpiSIM EpiSIM: simulation of multiple epistasis, linkage disequilibrium patterns and haplotype blocks for genome-wide interaction analysis
EvolSimulator A simulation test bed for hypotheses of genome evolution
EvolveAGene A realistic coding sequence simulation program that separates mutation from selection and allows the user to set selection conditions
fastsimcoal A continuous-­‐time coalescent simulator of genomic diversity under arbitrarily complex evolutionary scenarios
FastSLINK Simulation of Marker and Phenotype Data in Pedigrees
FFPopSim C++/Python library for population genetics.
FLUX SIMULATOR The Flux Simulator aims at providing a deterministic in silico reproduction of the experimental pipelines for RNA-Seq, employing a minimal set of parameters.
forqs Forward-in-time simulation of Recombination, Quantitative Traits, and Selection
ForSim ForSim: A Forward Evolutionary Computer Simulation
ForwSim The program given below is based on the algorithm described in Padhukasahasram et al. 2008 to simulate genetic drift in a standard Wright-Fisher process.
FPG Forward Population Genetic simulation
FREGENE FREGENE is a C++ program that simulates sequence-like data over large genomic regions in large diploid populations.
FIGG FIGG is a genome simulation tool that uses known or theorized variation frequency, per a given fragment size and grouped by GC content across a genome to model new genomes in FASTA format while tracking applied mutations for use in analysis
fwdpp A C++ template library for implementing efficient forward simulations.
GAMETES Genetic Architecture Model Emulator for Testing and Evaluating Software: Simulates complex SNP models with pure, strict epistatic interactions with n-loci.
GASP Genometric Analysis Simulation Program. A software tool for testing and investigating methods in statistical genetics by generating samples of family data based on user specified models.
GCTA Genome-wide Complex Trait Analysis
GemSIM Next generation sequencing read simulator
GeneArtisan Simulation of Markers in Case-Control Study Designs
GENOME A rapid coalescent-based whole genome simulator
GenomePop2 GenomePop2 is a specialization of the program GenomePop just to manage SNPs under more flexible and useful settings. If you need models with more than 2 alleles please use the GenomePop program version.
GenomeSimla GenomeSIMLA is currently under development- however, we have a beta release that we are asking to be tested
GENS2 Simulates interactions among two genetic and one environmental factor and also allows for epistatic interactions.
GWAsimulator A rapid whole genome simulation program
HAP-SAMPLE An association simulator for candidate regions or genome scans
HAPGEN A simulator for the simulation of case control datasets at SNP markers
HapSim A simulation tool for generating haplotype data with pre-specified allele frequencies and LD coefficients
HAPSIMU A program that simulates heterogeneous populations with various known and controllable structures under the continuous migration model or the discrete model
IBDsim IBDSim is a computer package for the simulation of genotypic data under general isolation by distance models.
indel-Seq-Gen A biological sequence simulation program that simulates highly divergent DNA sequences and protein superfamilies
Indelible A powerful and flexible simulator of biological evolution
invertFREGENE InvertFREGENE is a forward-in-time simulator of inversions in population genetic data
kernalPop A spatially explicit population genetic simulation engine
MaCS Markovian Coalescent Simulator
Mason A package for the simulation of nucleotide data.
mbs modifying Hudson's ms software to generate samples of DNA sequences with a biallelic site under selection
Mendel's Accountant Mendel's Accountant (MENDEL) is an advanced numerical simulation program for modeling genetic change over time and was developed collaboratively by Sanford, Baumgardner, Brewer, Gibson and ReMine
MetaSim A tool to generate collections of synthetic reads that reflect the diverse taxonomical composition of typical metagenome data sets
mlcoalsim Multilocus Coalescent Simulations
ms The purpose of this program is to allow one to investigate the statistical properties of such samples, to evaluate estimators or statistical tests, and generally to aid in the interpretation of polymorphism data sets.
msHOT The purpose of this program is to allow one to investigate the statistical properties of such samples, to evaluate estimators or statistical tests, and generally to aid in the interpretation of polymorphism data sets.
msms A coalescent Simlation tool with selection.
MySSP A program for the simulation of DNA sequence evolution across a phylogenetic tree
Nemo A forward-time, individual-based, genetically explicit, and stochastic simulation program designed to study the evolution of genetic markers, life history traits, and phenotypic traits in a flexible (meta-)population framework.
NetRecodon Coalescent simulation of coding DNA sequences with recombination (inter and intracodon), migration and demography
PEDAGOG Software for simulating eco-evolutionary population dynamics
phenosim A tool to add phenotypes to simulated genotypes
PhyloSim An R package for the Monte Carlo simulation of sequence evolution
pIRS Profile-based Illumina pair-end reads simulator
ProteinEvolver Simulation of protein evolution along phylogenies under structure-based substitution models
QMSim QTL and Marker Simulator
quantiNEMO An individual-based program for the analysis of quantitative traits with explicit genetic architecture potentially under selection in a structured population
RECOAL Simulates new haplotype data from a reference population of haplotypes.
Recodon Coalescent simulation of coding DNA sequences with recombination, migration and demography
rlsim A package for simulating RNA-seq library preparation with parameter estimation
Rmetasim Rmetasim is a front-end for the metasim engine that is implemented as a package that runs in the statistical computing environment R
RNA Seq Simulator RSS takes SAM alignment files from RNA-Seq data and simulates over dispersed, multiple replica, differential, non-stranded RNA-Seq datasets.
Rose Random model of sequence evolution
scrm A coalescent simulator optimized for long sequences and large samples.
SelSim SelSim is a program for Monte Carlo simulation of DNA polymorphism data for a recom- bining region within which a single bi-allelic site has experienced natural selection
Seq-Gen An application for the Monte Carlo simulation of molecular sequence evolution along phylogenetic trees.
SEQPower Statistical power analysis for sequence-based association studies
SeqSIMLA SeqSIMLA can simulate sequence data with user-specified disease and quantitative trait models. Family or unrelated case-control data can be simulated.
Serial NetEvolve A flexible utility for generating serially-sampled sequences along a tree or recombinant network
SFS_CODE SFS_CODE can perform forward population genetic simulations under a general Wright-Fisher model with arbitrary migration, demographic, selective, and mutational effects.
SIBSIM Quantitative phenotype simulation in extended pedigrees
SimAdapt A spatially explicit, individual-based, forward-time, landscape-genetic simulation model combined with a landscape cellular automaton.
SIMCOAL2 A coalescent program for the simulation of complex recombination patterns over large genomic regions under various demographic models
SimCopy An R package simulating the evolution of copy number profiles along a tree.
SIMLA SIMLA is a SIMuLAtion program that generates data sets of families for use in Linkage and Association studies.
SimPed A Simulation Program to Generate Haplotype and Genotype Data for Pedigree Structures
Simprot A program to simulate protein evolution by substitution, insertion and deletion
SimRare Rare variant simulation and analysis tool
simuGWAS A forward-time simulator that simulates realistic samples for genome-wide association studies.
simuPOP simuPOP is a general-purpose individual-based forward-time population genetics simulation environment.
SISSI A software tool to generate data of related sequences along a given phylogeny, taking into account user defined system of neighbourhoods and instantaneous rate matrices.
SMARTPOP Simulating Mating Alliance as a Reproductive Tactic for Populations
SNPsim Coalescent simulation of hotspot recombination
SPIP SPIP simulates the transmission of genes from parents to offspring in a population having demographic structure defined by the user
Splatche Spatial and Temporal Coalescences in Heterogeneous Environment
srv Simulator of Rare Varaints (srv) is a simulator for the simulation of the introduction and evolution of (rare) genetic variants.
SUP SLINK/FastSLINK utility program
TreesimJ A flexible, forward-time population genetic simulator
Vortex VORTEX is an individual-based simulation model for population viability analysis (PVA).
Wessim Whole Exome Sequencing SIMulator
Entering edit mode

Whoa! That's one gigantic list! Goes to show the richness in just any subdomain of bioinformatics.

Entering edit mode

5 years later, this is still a very impressive list. I have nothing to add but a (unfortunately not so recent) paper that might be of use as additional reference

Originally posted by @Joseph Hughes

Entering edit mode

BBMap's RandomReads: Generates single-ended or paired Illumina reads, or PacBio reads, from a genome. Also has a metagenome mode.

Entering edit mode
6.9 years ago
rtliu ★ 2.1k

Try read simulators in  including ART, wgsim etc


Login before adding your answer.

Traffic: 2078 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6