From Genome Analysis Wiki
Jump to navigationJump to search
282 bytes added
, 11:34, 26 September 2013
Line 19: |
Line 19: |
| --maxDepth ''threshold'' Positions where the read depth exceeds this threshold will be excluded. | | --maxDepth ''threshold'' Positions where the read depth exceeds this threshold will be excluded. |
| --reference Positions called as homozygous reference will be included in the output. | | --reference Positions called as homozygous reference will be included in the output. |
− | --uTsTv Use a uniform prior for transition-to-transversion ratio.
| |
| | | |
| To learn about default values for these options, simply run the program with no arguments. | | To learn about default values for these options, simply run the program with no arguments. |
Line 32: |
Line 31: |
| **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 |
− | **Transitions (C <-> T or A <-> G) accounts for most variants, while transversions account for minority of variants
| + | **transition has 2/3 probability |
− | **Uniform prior for transition to transversion ratio (1:1) | + | **each transversion has 1/6 probability |
| + | |
| + | *'''New implementation''': Alternative mutation model with uniform (uninformative) prior for transition to transversion ratio |
| + | **updated by Yancy Lo, 9/24/2012 |
| + | **each mutation has a 1/3 probability |
| + | **add --uniformTsTv in command line to enable this alternative mutation model |
| + | **download glfSingle with this new implementation here: [[File:Generic-glfSingle-2013-09-25.tar.gz]] |
| | | |
| == Download == | | == Download == |