96
edits
Changes
From Genome Analysis Wiki
→(step 2) Genotype/haplotype calling using thunder thunder_glf_freq
minDepth and maxDepth are the cutoffs on total depth (across all individuals). We have found it useful to exclude sites with extremely low and high total depth. Please see Important Filters below.
=== (step 2) Genotype/haplotype calling using thunder [https://wwwcsg.sph.umich.edu/csg/yli/thunder/thunder.V011.source.tgz thunder_glf_freq] ===
thunder_glf_freq --shotgun my.out.$chr --detailedInput -r 100 --states 200 --dosage --phase --interim 25 -o my.final.out
Notes:
(1) The program thunder used in step 2 is an extension of MaCH, the genotype imputation software we have previously developed. For details regarding the shared options, please check out [http://wwwcsg.sph.umich.edu/csg/yli/mach/index.html MaCH website] and [http://genome.sph.umich.edu/wiki/Mach MaCH wiki].
(2) Check out example files and command lines under examples/thunder/ in the thunder package [https://wwwcsg.sph.umich.edu/csg/yli/thunder/thunder.V011.source.tgz thunder_glf_freq].
== Example Showing the Whole Pipeline ==
