GotCloud: Variant Calling Pipeline

From Genome Analysis Wiki
Revision as of 17:49, 5 November 2012 by Mktrost (talk | contribs)
Jump to navigationJump to search

Back to the beginning [1]

The Variant Calling Pipeline (UMAKE) takes recalibrated BAM files and detects SNPs and calls their genotypes, producing VCF files.

Input Data:

  • Aligned/Processed/Recalibrated BAM files
  • Index file containing Sample IDs & BAM file names
  • Reference files
  • (Optional) Configuration file to override default options

BAM files

Index file

Reference Files

Configuration File

Running

Running umake is straightforward:

cd ~/myseq
/usr/local/biopipe/bin/umake --conf myconf ???
make -f [out-prefix].Makefile -j [# parallel jobs]