| Long Description || Fastsimcoal is a program to generate the neutral genomic molecular diversity in current or ancient samples drawn from a population with a complex demographic history. Fastsimcoal is a completely rewritten version of simcoal2 (Laval and Excoffier, 2004), a coalescent simulation program implementing a generation by generation approach while fastsimcoal is based on a much faster continuous time approximation. Despite a completely new coalescent engine, fastsimcoal uses exactly the same input files as simcoal2, and it produces very similar output files.
Fastsimcoal typically generates many replicates of random outcome of molecular diversity under a user-‐defined evolutionary scenario. The evolutionary scenario is defined in an input parameter file (extension .par) and the output diversity is written in arlequin project files (extension .arp) that can then be processed with arlequin or arlsumstat (Excoffier and Lischer, 2010) to get distributions of various summary statistics. Additional options of fastsimcoal can be specified on the command line (type "fastsimcoal -‐h" for help on command line options).
Fastsimcoal can handle very complex evolutionary scenarios including an arbitrary migration matrix between samples, historical events allowing for population resize, population fusion and fission, admixture events, changes in migration matrix, or changes in population growth rates. The time of sampling can be specified independently for each sample, allowing for serial sampling in the same or in different populations.
Different markers, such as DNA sequences, SNP, STR (microsatellite) or multi-‐locus allelic data can be generated under a variety of mutation models (e.g. finite-‐ and infinite-‐site models for DNA sequences, stepwise or generalized stepwise mutation model for STRs data, infinite-‐allele model for standard multi-‐allelic data).
fastsimcoal can simulate data in genomic regions with arbitrary recombination rates, thus allowing for recombination hotspots of different intensities at any position. fastsimcoal implements a new approximation to the ancestral recombination graph in the form of sequential Markov coalescent allowing it to very quickly generate genetic diversity for >100 Mb genomic segments. |