Correcting for DNA contamination in genotype calling
git clone https://github.com/hyunminkang/cleancall.git cd cleancall ./configure --prefix [/path/to/install] make make install [/path/to/install]/cctools pileup --loci [in.sites.vcf.gz] \\ --index [index.file.with.sample.id.and.bam.file.txt] \\ --out [output.pileup] --ref [ref.fa] --run [num.parallel.jobs] [/path/to/install]/cctools verify --index [output.pileup.ped] \\ --out [output.verify] --vcf [in.sites.vcf.gz] \\ --run [num.parallel.jobs] [/path/to/install]/cctools genotype --invcf [in.sites.vcf.gz] \\ --ped [output.verify.ped] --out [output.genotype.vcf.gz] \\ --ref [ref.fa] --run [num.parallel.jobs]
cleanCall is a software package for detecting, estimating, correcting for contamination in DNA sequence data. It contains one C++ program, cleanCall, as a core binary file that runs statistical algorithms. It contains three wrapper scripts, cctools-pileup for producing pileups from sequence data, cctools-verify for detecting and estimating contamination, and cctools-genotype for correcting for contaminations. These wrapper scripts allows to run cleanCall across a large number of sequenced genomes or exomes in a highly parallel manner.
The basic usage can be obtained by simply typing the command without any argument (e.g.
cctools-genotype). More detailed usage can be obtained with -man option (e.g.
cctools-genotype -man). If you have questions about cleanCall, you may send me an email at email@example.com
Flickinger M, Jun G, Abecasis GR, Boehnke M, Kang HM. Am J Hum Genet. (2015) Correcting for Sample Contamination in Genotype Calling of DNA Sequence Data. Am J Hum Genet 97(2):284-90. doi: 10.1016/j.ajhg.2015.07.002.
Jun G, Flickinger M, Hetrick KN, Romm JM, Doheny KF, Abecasis GR, Boehnke M, Kang HM (2012) Detecting and estimating contamination of human DNA samples in sequencing and array-based genotype data, Am J Hum Genet 91 (5), 839-848
Frequently Asked Questions (FAQs)
What is the purpose of the VCF file given to the -loci argument to cctools-pileup? What sort of data should this contain?
The variant site VCF file should contain VCF representing known common variant sites, such as HapMap, with allele frequency information available. For example, You can use HapMap3 variant sites available in our GotCloud Repository
What is the format of
The format of index file is simply tab-delimited files with the first column as sample ID and second column as full path of the BAM file corresponding to the sample.