Changes

From Genome Analysis Wiki
Jump to: navigation, search

Thunder

283 bytes added, 12:08, 21 February 2017
(step 2) Genotype/haplotype calling using thunder thunder_glf_freq
samtools pileup -g -T 1 -f ref.fa my.bam > my.glf
Note: you will need the reference fasta file ref.fa to create glf file from bam file.
== How to Run ==
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://wwwcsg.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
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.V010V011.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://www.sph.umich.edu/csg/yliabecasis/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://wwwsph.umich.edu/csg/abecasis/thunder/thunder.V011.source.tgz thunder_glf_freq]. == Example Showing the Whole Pipeline ==In the thunder [https://csg.sph.umich.edu/csg/yli/thunder/thunder.V010V011.source.tgz thunder_glf_freq]tarball, you can find under example/thunder/ folder, input files extracted from real data and a C-shell script that executes the whole analysis pipeline.
== Ligate Haplotypes ==
Please use [http://wwwcsg.sph.umich.edu/csg/yli/ligateHap.V004.tgz ligateHaplotypes].
== Important Filters ==
96
edits

Navigation menu