comparison somatic_sniper.xml @ 0:e81f5649bd4d draft default tip

Uploaded
author fcaramia
date Thu, 20 Jun 2013 00:06:01 -0400
parents
children
comparison
equal deleted inserted replaced
-1:000000000000 0:e81f5649bd4d
1 <tool id="somatic_sniper_tool" name="Somatic Sniper" version="1.0.2">
2 <description>: identify single nucleotide positions that are different between tumor and normal</description>
3 <requirements>
4 <requirement type="package" version="unstable">somatic-sniper</requirement>
5 </requirements>
6 <command interpreter="perl">
7 somatic_sniper_wrapper.pl
8
9 "NORMAL::$normal"
10 "TUMOR::$tumor"
11 "OUTPUT::$snp_output"
12 "OPTION::-F $output"
13
14 #if $option.option == "modify_parameters":
15 "OPTION::-q $option.readFilter"
16 "OPTION::-Q $option.somaticFilter"
17 "OPTION::-s $option.mutationPrior"
18 #if str($option.disablePriors) == "true"
19 "OPTION::-p"
20 #end if
21 #end if
22
23 "OPTION::-f $reference.fields.path"
24
25
26 </command>
27 <inputs>
28 <param name="reference" type="select" label="Select a reference genome">
29 <options from_data_table="all_fasta">
30 <filter type="sort_by" column="2" />
31 <validator type="no_options" message="No indexes are available" />
32 </options>
33 </param>
34
35 <param format="bam" name="normal" type="data" label="Normal sample" help=""/>
36 <param format="bam" name="tumor" type="data" label="Tumor Sample" help=""/>
37 <param name="output" type="select" label="Output Type" help="" optional="true">
38 <option value="classic" selected="true">Classic</option>
39 <option value="vcf">VCF</option>
40 <option value="bed">BED</option>
41 </param>
42
43
44 <conditional name="option">
45 <param name="option" type="select" label="Optional Parameters" help="" optional="true">
46 <option value="default_parameters" selected="true">Default Parameters</option>
47 <option value="modify_parameters">Modify Parameters</option>
48 </param>
49 <when value="modify_parameters">
50
51 <param name="readFilter" label="filtering reads with mapping quality less than" type="integer" value="0" optional="true" />
52 <param name="somaticFilter" label="filtering somatic snv output with somatic quality less than" type="integer" value="15" optional="true" />
53 <param name="disablePriors" type="select" label="disable priors in the somatic calculation. Increases sensitivity for solid tumors" help="" optional="true">
54 <option value="true" >true</option>
55 <option value="false" selected="true">false</option>
56 </param>
57 <param name="mutationPrior" label="prior probability of a somatic mutation" type="float" value="0.10000" optional="true" />
58
59 </when>
60
61 </conditional>
62
63 </inputs>
64 <outputs>
65 <data name="snp_output" format="text" label="${tool.name} result on ${on_string}" />
66 </outputs>
67 <help>
68 |
69
70
71 **Reference**
72
73 http://gmt.genome.wustl.edu/somatic-sniper/current/
74
75 -----
76
77 **What it does**
78
79 The purpose of this program is to identify single nucleotide positions that are different between tumor and normal
80 (or, in theory, any two bam files). It takes a tumor bam and a normal bam and compares the two to determine the
81 differences. It outputs a file in a format very similar to Samtools consensus format. It uses the genotype likelihood
82 model of MAQ (as implemented in Samtools) and then calculates the probability that the tumor and normal genotypes are
83 different. This probability is reported as a somatic score. The somatic score is the Phred-scaled probability (between 0 to 255)
84 that the Tumor and Normal genotypes are not different where 0 means there is no probability that the genotypes are different and
85 255 means there is a probability of 1 – 10(255/-10) that the genotypes are different between tumor and normal. This is consistent
86 with how the SAM format reports such probabilities.
87
88 bam-somaticsniper [options] -f ref.fasta tumor.bam normal.bam snp_output_file
89
90 Bam files must contain LB tag in @RG line.
91 Picard tools can be used to add lines to BAM headers.
92
93 -----
94
95 **Required Parameters**
96
97 ::
98
99 -f FILE REQUIRED reference sequence in the FASTA format
100
101 -----
102
103 **Options**
104
105 ::
106
107 -q INT filtering reads with mapping quality less than INT [0]
108
109 -Q INT filtering somatic snv output with somatic quality less than INT [15]
110
111 -p FLAG disable priors in the somatic calculation. Increases sensitivity for solid tumors
112
113 -J FLAG Use prior probabilities accounting for the somatic mutation rate
114
115 -s FLOAT prior probability of a somatic mutation (implies -J) [0.010000]
116
117 -T FLOAT theta in maq consensus calling model (for -c/-g) [0.850000]
118
119 -N INT number of haplotypes in the sample (for -c/-g) [2]
120
121 -r FLOAT prior of a difference between two haplotypes (for -c/-g) [0.001000]
122
123 -F STRING select output format [classic]
124 Available formats:
125 classic
126 vcf
127 bed
128
129 -----
130
131 **File Formats**
132
133 ::
134
135 Classic:
136
137 Each line contains the following tab-separated values:
138
139 1. Chromosome
140 2. Position
141 3. Reference base
142 4. IUB genotype of tumor
143 5. IUB genotype of normal
144 6. Somatic Score
145 7. Tumor Consensus quality
146 8. Tumor variant allele quality
147 9. Tumor mean mapping quality
148 10. Normal Consensus quality
149 11. Normal variant allele quality
150 12. Normal mean mapping quality
151 13. Depth in tumor (# of reads crossing the position)
152 14. Depth in normal (# of reads crossing the position)
153 15. Mean base quality of reads supporting reference in tumor
154 16. Mean mapping quality of reads supporting reference in tumor
155 17. Depth of reads supporting reference in tumor
156 18. Mean base quality of reads supporting variant(s) in tumor
157 19. Mean mapping quality of reads supporting variant(s) in tumor
158 20. Depth of reads supporting variant(s) in tumor
159 21. Mean base quality of reads supporting reference in normal
160 22. Mean mapping quality of reads supporting reference in normal
161 23. Depth of reads supporting reference in normal
162 24. Mean base quality of reads supporting variant(s) in normal
163 25. Mean mapping quality of reads supporting variant(s) in normal
164 26. Depth of reads supporting variant(s) in normal
165
166
167
168 </help>
169 </tool>
170
171