From Genome Analysis Wiki
Jump to navigationJump to search
12 bytes added
, 16:04, 15 October 2010
Line 15: |
Line 15: |
| This variant calling pipeline has two steps. (step 1) promotion of a set of potential polymorphisms; and (step 2) genotype/haplotype calling using LD information. | | This variant calling pipeline has two steps. (step 1) promotion of a set of potential polymorphisms; and (step 2) genotype/haplotype calling using LD information. |
| | | |
− | (step 1) Site promotion using software glfMultiples [https://www.sph.umich.edu/csg/yli/GPT_Freq.011.source.tgz GPT_Freq]. | + | === (step 1) Site promotion using software glfMultiples [https://www.sph.umich.edu/csg/yli/GPT_Freq.011.source.tgz GPT_Freq] === |
| | | |
| GPT_Freq -b my.out -p 0.9 --minDepth 10 --maxDepth 1000 *.glf | | GPT_Freq -b my.out -p 0.9 --minDepth 10 --maxDepth 1000 *.glf |
Line 21: |
Line 21: |
| 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. | | 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://www.sph.umich.edu/csg/yli/thunder/thunder.V009.source.tgz thunder_glf_freq]. | + | === (step 2) Genotype/haplotype calling using thunder [https://www.sph.umich.edu/csg/yli/thunder/thunder.V009.source.tgz thunder_glf_freq] === |
| | | |
| thunder_glf_freq --shotgun my.out.$chr -r 100 --states 200 --dosage --phase --interim 25 -o my.final.out | | thunder_glf_freq --shotgun my.out.$chr -r 100 --states 200 --dosage --phase --interim 25 -o my.final.out |