27
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to search→Model for Variant Calling
**P(non-reference heterozygote) = 0.01 * 2/3
**P(non-reference heterozygote) = 0.01 * 2/3
**P(non-reference homozygote) = 0.01 * 1/3
**P(non-reference homozygote) = 0.01 * 1/3
*Two options of mutation model:
*Mutation model: Transitions (C <-> T or A <-> G) accounts for most variants, while transversions account for minority of variants
**transition has 2/3 probability
**each transversion has 1/6 probability
**Uniform (uninformative) prior for transition to transversion ratio
*New implementation: Alternative mutation model with uniform (uninformative) prior for transition to transversion ratio
**each mutation has a 1/3 probability
**add --uniformTsTv in command line to enable this alternative mutation model
**download here
== Download ==
== Download ==
