This software takes paired allele counts from somatic and germline samples to estimate the distortion effect size and the map position of segregation distortion elements in the germline sample.
It takes a single input data file with read counts supporting either parental chromosome by site. The file is tab delimited, with no header line and has the format:
- Chromosome ID.
- Position in basepairs.
- Map position in centimorgans.
- Read counts for the first parental type in the somatic library.
- Read counts for the second parental type in the somatic library.
- Read counts for the first parental type in the gamete library.
- Read counts for the second parental type in the gamete library.
Basic Usage:
map_sd -d <data_file>
This will produce an output file with the following format
- Chromosome ID
- Distorter Position
- Ancestry ratio in the somatic library.
- Ancestry ratio in the gamete library.
- lnL of somatic ratio in gamete library.
- lnL of the gamete ratio in the gamete library.
Additional Options/Advanced Usage
-w [int] will produce windowed estimates of SD across each chromosome with window size -w
-b [int] will produce bootstrap estimates of SD position via resampling
-e [float] uniform error rate per read
-t [float] terminate optimation if successive points are within likelihood units.