annotate readdepth2sequencingdepth.xml @ 0:07588b899c13 draft

Uploaded
author arkarachai-fungtammasan
date Wed, 01 Apr 2015 17:05:51 -0400
parents
children d5ed5c2e25c3
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
1 <tool id="readdepth2seqdepth" name="Convert informative read depth to sequencing depth" version="1.0.0">
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
2 <description>for flank-based mapping of microsatellites</description>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
3 <command interpreter="python2.7">sequencingdepthconversion_G.py $repeatlength $flanksize $readlength $infodepth $probprediction > $output </command>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
4
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
5 <inputs>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
6 <param name="repeatlength" type="integer" value="10" label="Repeat length (bp)" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
7 <param name="flanksize" type="integer" value="20" label="Required flank bases on each side in mapping" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
8 <param name="readlength" type="integer" value="100" label="Read length (treat all read as single end read)" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
9 <param name="infodepth" type="integer" value="5" label="Required read depth" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
10 <param name="probprediction" type="float" value="0.9" label="Proportion of genome that need certain level of read depth" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
11 </inputs>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
12 <outputs>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
13 <data format="input" name="output" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
14
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
15 </outputs>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
16 <tests>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
17 <!-- Test data with valid values -->
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
18 <test>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
19 <param name="repeatlength" value="10"/>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
20 <param name="flanksize" value="20" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
21 <param name="readlength" value="100" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
22 <param name="infodepth" value="5" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
23 <param name="probprediction" value="0.9" />
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
24 <output name="output" file="readdepth2seqdepth.out"/>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
25 </test>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
26
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
27 </tests>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
28 <help>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
29
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
30
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
31 .. class:: infomark
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
32
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
33 **What it does**
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
34
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
35 This tool is used to convert informative read depth (specified by user) to sequencing depth when the microsatellites is mapped using TRFM pipeline.
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
36 The locus specific sequencing depth is the sequencing depth that will make a certain loci have certain read depth based on uniform mapped of read. It is calculated as: ::
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
37
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
38 yrequired = ( X * L ) / (L - (2F+r-1))
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
39
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
40 Where X = read depth, L = read length, F = the number of flanked bases required on each flanking regions, r = the expected repeat length of microsatellite of interest.
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
41
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
42 The genome wide sequencing depth is the sequencing depth that will make certain percentage of genome (e.g. 90 percent or 95 percent) to have certain locus specific sequencing depth. It's calculated using numerical guessing to find smallest lambda that: ::
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
43
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
44 0.90 (or other proportion specified by user) &lt; = P(Y=0) + P(Y=1) + …+ P(Y=yrequired-1)
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
45
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
46 P(Y=y) = (lambda^(y) * e ^(-lambda)) /y!
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
47
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
48 y = specific level of sequencing depth. Lambda = genome wide sequencing depth
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
49
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
50
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
51 **Citation**
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
52
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
53 When you use this tool, please cite **Fungtammasan A, Ananda G, Hile SE, Su MS, Sun C, Harris R, Medvedev P, Eckert K, Makova KD. 2015. Accurate Typing of Short Tandem Repeats from Genome-wide Sequencing Data and its Applications, Genome Research**
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
54
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
55
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
56 </help>
07588b899c13 Uploaded
arkarachai-fungtammasan
parents:
diff changeset
57 </tool>