--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/varscan_mpileup.xml	Sat Nov 16 18:48:12 2013 -0500
@@ -0,0 +1,155 @@
+<tool id="varscan" name="Varscan" version="0.1">
+    <description>for variant detection</description>
+
+    <requirements>
+        <requirement type="package" version="2.3.6">varscan</requirement>
+    </requirements>
+
+    <!-- 
+    The version command string is not yet a template that can be filled in, so version command is not yet possible.
+    <version_command>java -jar ${GALAXY_DATA_INDEX_DIR}/shared/jars/varscan/VarScan.jar 2>&amp;1 | head -n 1</version_command>
+    -->
+
+    <command>
+        ## Set up samples list file.
+        #if $sample_names.strip() != '':
+	       echo $sample_names | awk -F ',' '{ for (i = 1; i &lt;= NF; i++) { print \$i; } }' > samples_list.txt;
+        #end if
+
+        ## Set up command + input.
+        java -jar \$JAVA_JAR_PATH/VarScan.v2.3.6.jar ${cmd} ${input}
+        --min-coverage ${min_coverage} 
+        --min-reads2 ${min_supporting_reads} 
+        --min-avg-qual ${min_avg_qual}
+        --min-var-freq ${min_var_freq}
+        --min-freq-for-hom ${min_freq_for_hom}
+        --p-value ${p_value}
+        #if str($strand_filter) == 'yes':
+          --strand-filter 1
+        #end if
+
+        ## Report only variants in consensus.
+        #if str($cmd) == 'mpileup2cns':
+          --variants
+        #end if
+        
+        ## Set up outputs.
+        --output-vcf 1 > $output
+
+        #if $sample_names.strip() != '':
+            --vcf-sample-list samples_list.txt
+        #end if
+    </command>
+
+    <inputs>
+        <param format="pileup" name="input" type="data" label="Pileup dataset" help=""/>
+
+        <param name="cmd" type="select" label="Analysis type">
+          <option value="mpileup2snp" selected="True">single nucleotide variation</option>
+          <option value="mpileup2indel">insertions and deletions</option>
+          <option value="mpileup2cns">consensus genotype</option>
+        </param>
+
+        <param name="min_coverage" type="integer" value="8" min="1" max="200" label="Minimum read depth" help="Minimum depth at a position to make a call"/>
+        <param name="min_supporting_reads" type="integer" value="2" min="1" max="200" label="Minimum supporting reads" help="Minimum supporting reads at a position to make a call"/>
+        <param name="min_avg_qual" type="integer" value="15" min="1" max="50" label="Minimum base quality at a position to count a read"/>
+        <param name="min_var_freq" type="float" value="0.01" min="0" max="1" label="Minimum variant allele frequency threshold"/>
+        <param name="min_freq_for_hom" type="float" value="0.75" min="0" max="1" label="Minimum frequency to call homozygote"/>
+        <param name="p_value" type="float" value="0.99" min="0" max="1" label="p-value threshold for calling variants"/>
+        <param name="strand_filter" type="select" label="Ignore variants with >90% support on one strand">
+            <option value="no" selected="True">no</option>
+            <option value="yes">yes</option>
+        </param>
+        <param name="sample_names" type="text" value="" help="Separate sample names by comma; leave blank to use default sample names."/>
+    </inputs>
+
+    <stdio>
+        <regex match="Exception" source="both" level="fatal" description="Tool exception"/>
+        <regex match=".*" source="both" level="log" description="tool progress"/>
+    </stdio>
+
+    <outputs>
+        <data name="output" format="vcf"/>
+    </outputs>
+
+    <trackster_conf>
+    </trackster_conf>
+
+    <tests>
+        <test>
+            <param name="input" value="test_in1.pileup" />
+            <param name="cmd" value="mpileup2cns" />
+            <param name="min_coverage" value="8" />
+            <param name="min_supporting_reads" value="2" />
+            <param name="min_avg_qual" value="15" />
+            <param name="min_var_freq" value="0.01" />
+            <param name="min_freq_for_hom" value="0.75" />
+            <param name="p_value" value="0.99" />
+            <param name="strand_filter" value="no" />
+            <param name="sample_names" value="" />
+            <output name="output" file="test_out1.vcf" lines_diff="0" />
+        </test>
+    </tests>
+
+    <help>
+**VarScan Overview**
+
+VarScan_ performs variant detection for massively parallel sequencing data, such as exome, WGS, and transcriptome data. It calls variants from a mpileup dataset and produces a VCF 4.1 Full documentation is available online_.
+
+Please cite: Koboldt, D., Zhang, Q., Larson, D., Shen, D., McLellan, M., Lin, L., Miller, C., Mardis, E., Ding, L., and Wilson, R. (2012). VarScan 2: Somatic mutation and copy number alteration discovery in cancer by exome sequencing Genome Research DOI: 10.1101/gr.129684.111 
+
+.. _VarScan: http://varscan.sourceforge.net/
+.. _online: http://varscan.sourceforge.net/using-varscan.html
+
+**Input**
+
+::
+
+  mpileup file - The SAMtools mpileup file
+ 
+
+**Output**
+
+VarScan produces a VCF 4.1 dataset as output.
+
+**Parameters**
+
+::
+
+  analysis type
+    single nucleotide detection     Identify SNPs from an mpileup file
+    insertions and deletion       Identify indels an mpileup file
+    consensus genotype     Call consensus and variants from an mpileup file
+
+  min-coverage  
+    Minimum read depth at a position to make a call [8]
+
+  min-reads2    
+    Minimum supporting reads at a position to call variants [2]
+
+  min-avg-qual  
+    Minimum base quality at a position to count a read [15]
+
+  min-var-freq  
+        Minimum variant allele frequency threshold [0.01]
+
+  min-freq-for-hom
+    Minimum frequency to call homozygote [0.75]
+  
+  p-value
+    Default p-value threshold for calling variants [99e-02]
+  
+  strand-filter
+    Ignore variants with >90% support on one strand [1]
+  
+  output-vcf
+    If set to 1, outputs in VCF format
+  
+  vcf-sample-list
+    For VCF output, a list of sample names in order, one per line
+  
+  variants
+    Report only variant (SNP/indel) positions [0]
+
+    </help>
+</tool>