# HG changeset patch # User fcaramia # Date 1371185067 14400 # Node ID f8595c0382f66bbb30789fc907855395e386aef6 # Parent 84d0aae471b31ddc4db86313fdb15aecf7649478 Uploaded diff -r 84d0aae471b3 -r f8595c0382f6 varscan_somatic.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/varscan_somatic.xml Fri Jun 14 00:44:27 2013 -0400 @@ -0,0 +1,131 @@ + + + somatic mutation caller for cancer genomics + + + VarScan + + + varscan_somatic.pl + "COMMAND::java -jar \$JAVA_JAR_PATH/VarScan.v2.3.5.jar somatic" + "NORMAL::$normal" + "TUMOR::$tumor" + "OUTPUT::$output" + + "OPTION::--min-coverage $min_coverage" + "OPTION::--min-coverage-normal $min_coverage_normal" + "OPTION::--min-coverage-tumor $min_coverage_tumor" + + "OPTION::--min-var-freq $min_var_freq" + "OPTION::--min-freq-for-hom $min_freq_for_hom" + + "OPTION::--normal-purity $normal_purity" + "OPTION::--tumor-purity $tumor_purity" + + "OPTION::--p-value $p_value" + "OPTION::--somatic-p-value $somatic_p_value" + + "OPTION::--strand-filter $strand_filter" + "OPTION::--validation $validation" + "OPTION::--output-vcf 1" + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +.. class:: infomark + +**What it does** + +:: + + VarScan is a platform-independent mutation caller for targeted, exome, and whole-genome resequencing data generated on Illumina, SOLiD, Life/PGM, Roche/454, and similar instruments. The newest version, VarScan 2, is written in Java, so it runs on most operating systems. It can be used to detect different types of variation: + + Germline variants (SNPs an dindels) in individual samples or pools of samples. + Multi-sample variants (shared or private) in multi-sample datasets (with mpileup). + Somatic mutations, LOH events, and germline variants in tumor-normal pairs. + Somatic copy number alterations (CNAs) in tumor-normal exome data. + + +**Input** + +:: + + mpileup normal file - The SAMtools mpileup file for normal + mpileup tumor file - The SAMtools mpileup file for tumor + + +**Parameters** + +:: + + min-coverage + Minimum read depth at a position to make a call [8] + + min-coverage-normal + Minimum coverage in normal to call somatic [8] + + min-coverage-tumor + Minimum coverage in tumor to call somatic [6] + + min-var-freq + Minimum variant frequency to call a heterozygote [0.10] + + min-freq-for-hom + Minimum frequency to call homozygote [0.75] + + normal-purity + Estimated purity (non-tumor content) of normal sample [1.00] + + tumor-purity + Estimated purity (tumor content) of tumor sample [1.00] + + p-value + Default p-value threshold for calling variants [0.99] + + somatic-p-value + P-value threshold to call a somatic site [0.05] + + strand-filter + If set to 1, removes variants with >90% strand bias + + validation + If set to 1, outputs all compared positions even if non-variant + + output-vcf + If set to 1, outputs in VCF format [Default] + + + + + +