Mercurial > repos > geert-vandeweyer > varscan_wrapper
diff varscan/varscan_processSomatic.xml @ 0:848f3dc54593 draft
Uploaded
author | geert-vandeweyer |
---|---|
date | Fri, 07 Mar 2014 06:17:32 -0500 |
parents | |
children |
line wrap: on
line diff
--- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/varscan/varscan_processSomatic.xml Fri Mar 07 06:17:32 2014 -0500 @@ -0,0 +1,154 @@ +<tool id="varscan_processSomatic" name="ProcessSomatic" version="2.3.5"> + <description> + Extract HC calls from Somatic Caller + </description> + <requirements> + <requirement type="package" version="2.3.5">VarScan</requirement> + </requirements> + <command interpreter="perl"> + + varscan_processSomatic.pl + "COMMAND::java -jar \$JAVA_JAR_PATH/VarScan.v2.3.5.jar processSomatic" + "INPUT::$input" + "LOG::$log" + + "OPTION::--min-tumor-freq $min_tumor_freq" + "OPTION::--max-normal-freq $max_normal_freq" + "OPTION::--p-value $p_value" + "OUTPUT::--loh $loh" + "OUTPUT::--loh_hc $loh_hc" + "OUTPUT::--germ $germ" + "OUTPUT::--germ_hc $germ_hc" + "OUTPUT::--som $som" + "OUTPUT::--som_hc $som_hc" + "OUTPUT::--som_hc_vcf $som_hc_vcf" + "OUTPUT::--germ_hc_vcf $germ_hc_vcf" + "OUTPUT::--loh_hc_vcf $loh_hc_vcf" + + </command> + + <inputs> + + <param name="input" type="data" format="txt" label="Input File" help="VarScan Somatics output file in native format" /> + + <param name="min_tumor_freq" type="float" label="min-tumor-freq" help="Minimum variant allele frequency in tumor" optional="true" value="0.10"/> + <param name="max_normal_freq" type="float" label="max-normal-freq" help="Maximum variant allele frequency in normal" optional="true" value="0.05"/> + <param name="p_value" type="text" label="p-value" help="P-value for high-confidence calling" optional="true" value="0.07"/> + <param name="outtype" type="select" label="Output Type:" default="1"> + <option value="0">Native VarScan Tables</option> + <option value="1">VCF format (only High Confidence)</option> + <option value="2">Both</option> + </param> + + </inputs> + <outputs> + <data type="data" format="txt" name="log" label="${tool.name} result on ${on_string} (log) " /> + <data type="data" format="txt" name="loh" label="${tool.name} result on ${on_string} (loh) " > + <filter>outtype != "1"</filter> + </data> + <data type="data" format="txt" name="loh_hc" label="${tool.name} result on ${on_string} (loh_hc)"> + <filter>outtype != "1"</filter> + </data> + + <data type="data" format="txt" name="germ" label="${tool.name} result on ${on_string} (germline)" > + <filter>outtype != "1"</filter> + </data> + + <data type="data" format="txt" name="germ_hc" label="${tool.name} result on ${on_string} (germline_hc)"> + <filter>outtype != "1"</filter> + </data> + + <data type="data" format="txt" name="som" label="${tool.name} result on ${on_string} (somatic)" > + <filter>outtype != "1"</filter> + </data> + + <data type="data" format="txt" name="som_hc" label="${tool.name} result on ${on_string} (somatic_hc)" > + <filter>outtype != "1"</filter> + </data> + + <data type="data" format="vcf" name="som_hc_vcf" label="${tool.name} result on ${on_string} (Somatic_HC.vcf)" > + <filter>outtype != "0"</filter> + </data> + + <data type="data" format="vcf" name="loh_hc_vcf" label="${tool.name} result on ${on_string} (LOH_HC.vcf)" > + <filter>outtype != "0"</filter> + </data> + + <data type="data" format="vcf" name="germ_hc_vcf" label="${tool.name} result on ${on_string} (Germline_HC.vcf)" > + <filter>outtype != "0"</filter> + </data> + + + + + </outputs> + + <help> + +.. 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] + + + + </help> +</tool> +