annotate oases_optimiser.xml @ 4:42fad5f9c78e draft

planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
author drosofff
date Mon, 19 Oct 2015 11:56:34 -0400
parents b4d1940ecc90
children 58bf34f98232
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
1 <tool id="oasesoptimiserv" name="Oases_optimiser" version="1.1.4">
1
529f5ef2bfa4 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 0
diff changeset
2 <description>Auto optimise de novo RNA-seq Oases/Velvet assembly</description>
0
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
3 <requirements>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
4 <requirement type="package" version="0.2.08_7a32460a60929b510037952ae815bb6e29b68123">oases</requirement>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
5 </requirements>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
6 <command interpreter="python">
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
7 oases_optimiser.py '$start_hash_length' '$end_hash_length'
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
8 '#for $i in $inputs
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
9 ${i.input}
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
10 #end for
4
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
11 '
0
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
12 '$transcripts'
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
13 </command>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
14
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
15 <inputs>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
16 <param label="Start Hash Length" name="start_hash_length" type="select" help="k-mer length in base pairs of the words being hashed. Shorter hash lengths (i.e. less than 31) may cause out-of-memory problems.">
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
17 <option value="11" selected="yes">11</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
18 <option value="13">13</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
19 <option value="15">15</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
20 <option value="17">17</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
21 <option value="19">19</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
22 <option value="21">21</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
23 <option value="23">23</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
24 <option value="25">25</option>
4
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
25 <option value="35">35</option>
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
26 <option value="45">45</option>
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
27 <option value="55">55</option>
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
28 <option value="65">65</option>
0
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
29 </param>
4
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
30 <param label="End Hash Length" name="end_hash_length" type="select" help="k-mer length in base pairs of the words being hashed. Value has to be higher than the Start Hash Length value">
0
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
31 <option value="25">25</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
32 <option value="27">27</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
33 <option value="29">29</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
34 <option value="31">31</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
35 <option value="33">33</option>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
36 <option value="35">35</option>
4
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
37 <option value="45">45</option>
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
38 <option value="55">55</option>
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
39 <option value="63">63</option>
42fad5f9c78e planemo upload for repository https://bitbucket.org/drosofff/gedtools/ commit 4ea89b80296509604d1b224ee1671a6497ddca14
drosofff
parents: 3
diff changeset
40 <option value="69">69</option>
0
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
41 </param>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
42
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
43
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
44
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
45 <repeat name="inputs" title="Input Files">
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
46 <param name="input" type="data" format="fasta" label="Dataset: short reads, fasta format"/>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
47 </repeat>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
48 </inputs>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
49
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
50 <outputs>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
51 <data format="fasta" name="transcripts" label="${tool.name} on ${on_string}: Denovo assembled transcripts"/>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
52 </outputs>
2
6d3b484e43b0 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 1
diff changeset
53 <tests>
6d3b484e43b0 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 1
diff changeset
54 <test>
6d3b484e43b0 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 1
diff changeset
55 <param name="input" value="input.fa" ftype="fasta"/>
6d3b484e43b0 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 1
diff changeset
56 <param name="start_hash_length" value="15" />
6d3b484e43b0 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 1
diff changeset
57 <param name="end_hash_length" value="35" />
6d3b484e43b0 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 1
diff changeset
58 <output name="transcripts" file="output.fa"/>
6d3b484e43b0 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 1
diff changeset
59 </test>
6d3b484e43b0 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents: 1
diff changeset
60 </tests>
0
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
61
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
62 <help>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
63 **Oases Optimiser Overview**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
64
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
65 Oases_ is a de novo transcriptome assembler specially designed to work with short reads. It is an extension of Velvet developed by Daniel Zerbino and Ewan Birney at the European Bioinformatics Institute (EMBL-EBI), near Cambridge, in the United Kingdom.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
66
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
67 Provide all the information about insert lengths and their standard deviation as
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
68 possible (identical to Velvet). In the advanced oases options you should set the -ins_length flag when using paired-end reads (e.g. -ins_length 200 -ins_length_sd 100). If you do not do this, then paired-end information will not be used by Oases!
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
69
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
70 ----------------------------------------------------------------------------------
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
71
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
72 Oases produces a number of output files, which correspond to the different algorithms
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
73 being run succesively on the data. In the above example, you would find:
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
74
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
75 transcripts.fa
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
76 A FASTA file containing the transcripts imputed directly from trivial
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
77 clusters of contigs (loci with less than two transcripts and Confidence Values = 1)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
78 and the highly expressed transcripts imputed by dynamic
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
79 programming (loci with more than 2 transcripts and Confidence Values less than 1).
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
80
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
81 splicing_events.txt
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
82 A hybrid file which describes the contigs contained in each locus in FASTA
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
83 format, interlaced with one line descriptions of splicing events using the
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
84 AStalavista nomenclature*.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
85
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
86 Read the Oases `documentation`__ for details on using the Oases Assembler.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
87
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
88 .. _Velvet: http://www.ebi.ac.uk/~zerbino/oases/
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
89
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
90 .. __: http://www.ebi.ac.uk/~zerbino/oases/Manual.txt
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
91
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
92 ------
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
93
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
94 **Other outputs (as per Velvet)**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
95
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
96
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
97 **Contigs**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
98
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
99 The *contigs.fa* file.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
100 This fasta file contains the sequences of the contigs longer than 2k, where k is the word-length used in velveth. If you have specified a min contig lgth threshold, then the contigs shorter than that value are omitted.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
101 Note that the length and coverage information provided in the header of each contig should therefore be understood in k-mers and in k-mer coverage (cf. 5.1) respectively.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
102 The N's in the sequence correspond to gaps between scaffolded contigs. The number of N's corresponds to the estimated length of the gap. For reasons of compatibility with the archives, any gap shorter than 10bp is represented by a sequence of 10 N's.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
103
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
104 **Stats**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
105
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
106 The *stats.txt* file.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
107 This file is a simple tabbed-delimited description of the nodes. The column names are pretty much self-explanatory. Note however that node lengths are given in k-mers. To obtain the length in nucleotides of each node you simply need to add k - 1, where k is the word-length used in velveth.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
108 The in and out columns correspond to the number of arcs on the 5' and 3' ends of the contig respectively.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
109 The coverages in columns short1 cov, short1 Ocov, short2 cov, and short2 Ocov are provided in k-mer coverage (5.1).
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
110 Also, the difference between # cov and # Ocov is the way these values are computed. In the first count, slightly divergent sequences are added to the coverage tally. However, in the second, stricter count, only the sequences which map perfectly onto the consensus sequence are taken into account.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
111
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
112 **LastGraph**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
113
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
114 The *LastGraph* file.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
115 This file describes in its entirety the graph produced by Velvet.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
116
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
117 **AMOS.afg**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
118 The *velvet_asm.afg* file.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
119 This file is mainly designed to be read by the open-source AMOS genome assembly package. Nonetheless, a number of programs are available to transform this kind of file into other assembly file formats (namely ACE, TIGR, Arachne and Celera). See http://amos.sourceforge.net/ for more information.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
120 The file describes all the contigs contained in the contigs.fa file (cf 4.2.1).
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
121
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
122 **Advanced options**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
123 -ins_length integer : expected distance between two paired-end reads in the first short-read dataset (default: no read pairing)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
124 -ins_length2 integer : expected distance between two paired-end reads in the second short-read dataset (default: no read pairing)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
125 -ins_length_long integer : expected distance between two long paired-end reads (default: no read pairing)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
126 -ins_length*_sd integer : est. standard deviation of respective dataset (default: 10% of corresponding length)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
127 -scaffolding yes|no : scaffolding of contigs used paired end information (default: on)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
128 -max_branch_length integer : maximum length in base pair of bubble (default: 100)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
129 -max_divergence floating-point : maximum divergence rate between two branches in a bubble (default: 0.2)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
130 -max_gap_count integer : maximum number of gaps allowed in the alignment of the two branches of a bubble (default: 3)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
131 -min_pair_count integer : minimum number of paired end connections to justify the scaffolding of two long contigs (default: 10)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
132 -max_coverage floating point : removal of high coverage nodes AFTER tour bus (default: no removal)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
133 -long_mult_cutoff int : minimum number of long reads required to merge contigs (default: 2)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
134 -min_trans_length
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
135 simple threshold on output transcript length
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
136 -cov_cutoff
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
137 minimum number of times a k-mer has to be observed to be used in the
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
138 assembly (just like in Velvet) [default=3]
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
139 -min_pair_cov
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
140 minimum number of times two contigs must be connected by reads or read pairs to be clustered together [default=4]
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
141 -paired_cutoff
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
142 minimum ratio between the numbers of observed and expected connecting
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
143 read pairs between two contigs [default=0.1]
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
144
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
145
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
146
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
147 **Hash Length**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
148
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
149 The hash length, also known as k-mer length, corresponds to the length, in base pairs, of the words being hashed.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
150
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
151 The hash length is the length of the k-mers being entered in the hash table. Firstly, you must observe three technical constraints::
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
152
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
153 # it must be an odd number, to avoid palindromes. If you put in an even number, Velvet will just decrement it and proceed.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
154 # it must be below or equal to MAXKMERHASH length (cf. 2.3.3, by default 31bp), because it is stored on 64 bits
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
155 # it must be strictly inferior to read length, otherwise you simply will not observe any overlaps between reads, for obvious reasons.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
156
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
157 Now you still have quite a lot of possibilities. As is often the case, it's a trade- off between specificity and sensitivity. Longer kmers bring you more specificity (i.e. less spurious overlaps) but lowers coverage (cf. below). . . so there's a sweet spot to be found with time and experience.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
158 We like to think in terms of "k-mer coverage", i.e. how many times has a k-mer been seen among the reads. The relation between k-mer coverage Ck and standard (nucleotide-wise) coverage C is Ck = C # (L - k + 1)/L where k is your hash length, and L you read length.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
159 Experience shows that this kmer coverage should be above 10 to start getting decent results. If Ck is above 20, you might be "wasting" coverage. Experience also shows that empirical tests with different values for k are not that costly to run! VelvetOptimiser automates these tests for you.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
160
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
161
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
162 **Velvetg options**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
163
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
164
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
165
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
166 **Input Files**
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
167
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
168 This Velvet wrapper works only with fasta format, non-paired-end reads.
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
169
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
170 Supported file formats are::
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
171
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
172 fasta
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
173
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
174 Read categories are::
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
175
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
176 short (default)
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
177
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
178 </help>
435711ba68e3 planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
diff changeset
179 </tool>