Changes

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 ==
27

edits

Navigation menu