Difference between revisions of "Thunder"
| Line 1: | Line 1: | ||
| − | This page documents how to perform variant calling from low-coverage sequencing data using glfmultiples and thunder. | + | This page documents how to perform variant calling from low-coverage sequencing data using glfmultiples and thunder. The pipeline was originally developed by [mailto:yunli@med.unc.edu Yun Li] for the 1000 Genomes Low Coverage Pilot Project. |
== Input Data == | == Input Data == | ||
| Line 12: | Line 12: | ||
== How to Run == | == 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://www.sph.umich.edu/csg/yli/GPT_Freq.011.source.tgz GPT_Freq]. | ||
| + | |||
== Important Filters == | == Important Filters == | ||
Revision as of 12:43, 15 October 2010
This page documents how to perform variant calling from low-coverage sequencing data using glfmultiples and thunder. The pipeline was originally developed by Yun Li for the 1000 Genomes Low Coverage Pilot Project.
Input Data
To get started, you will need glf files in the standard format glf format.
If you do not have glf files, you can generate them from bam files (bam format also specified in glf format bam format) using the following command line:
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 GPT_Freq.
