This summary statistic returns the actual DNA sequences from finite sites simulations. It can not be calculated together with other summary statistics or when assuming an infinite sites mutation model. No outgroup is needed for it, and the outgroup sequences will also be returned if present.

1 | ```
sumstat_dna(name = "dna", transformation = identity)
``` |

`name` |
The name of the summary statistic. When simulating a model, the value of the statistics are written to an entry of the returned list with this name. Summary statistic names must be unique in a model. |

`transformation` |
An optional function for transforming the results of the statistic. If specified, the results of the transformation are returned instead of the original values. |

A list of sequences for each locus. Each entries is a character matrix decoding the sequences. Each row is an individual, and each column is a genetic position.

To create a demographic model: `coal_model`

To calculate this statistic from data: `calc_sumstats_from_data`

Other summary statistics: `sumstat_file`

,
`sumstat_four_gamete`

,
`sumstat_ihh`

, `sumstat_jsfs`

,
`sumstat_mcmf`

,
`sumstat_nucleotide_div`

,
`sumstat_omega`

,
`sumstat_seg_sites`

,
`sumstat_sfs`

,
`sumstat_tajimas_d`

,
`sumstat_trees`

1 2 3 4 | ```
model <- coal_model(5, 1, 10) +
feat_mutation(5, model = "GTR", gtr_rates = rep(1, 6)) +
sumstat_dna()
## Not run: simulate(model)$dna
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.