comparison readdepth2sequencingdepth.xml @ 5:b27006b0a953

update to latest version
author devteam@galaxyproject.org
date Wed, 22 Apr 2015 12:19:28 -0400
parents 20ab85af9505
children
comparison
equal deleted inserted replaced
4:ecfc9041bcc5 5:b27006b0a953
30 30
31 .. class:: infomark 31 .. class:: infomark
32 32
33 **What it does** 33 **What it does**
34 34
35 This tool is used to convert informative read depth (specified by user) to sequencing depth when the microsatellites is mapped using TRFM pipeline. 35 This tool is used to convert informative read depth (specified by user) to sequencing depth when the STRs is mapped using STR-FM pipeline.
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: :: 36 The locus specific sequencing depth (yrequired) is the sequencing depth that will make an STR locus to have a certain informative read depth based on uniform mapping of reads. It is calculated as follows: ::
37 37
38 yrequired = ( X * L ) / (L - (2F+r-1)) 38 yrequired = ( X * L ) / (L - (2F+r-1))
39 39
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. 40 where X = informative read depth, L = read length, F = the number of flanking bases required on either side, r = the expected repeat length of the STR of interest.
41 41
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: :: 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: ::
43 43
44 0.90 (or other proportion specified by user) < = P(Y=0) + P(Y=1) + …+ P(Y=yrequired-1) 44 0.90 (or other proportion specified by user) < = P(Y=0) + P(Y=1) + …+ P(Y=yrequired-1)
45 45
46 P(Y=y) = (lambda^(y) * e ^(-lambda)) /y! 46 where P(Y=y) = (lambda^(y) * e ^(-lambda)) /y!
47 47
48 y = specific level of sequencing depth. Lambda = genome wide sequencing depth 48 y = specific level of sequencing depth. Lambda = genome wide sequencing depth
49 49
50 Please refer the Methods section of the paper cited below for further details.
50 51
51 **Citation** 52 **Citation**
52 53
53 When you use this tool, please cite **Arkarachai Fungtammasan and Guruprasad Ananda (2014).** 54 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**
54 55
55 56
56 </help> 57 </help>
57 </tool> 58 </tool>