annotate repex_full_clustering.xml @ 13:789a5e76c77c draft default tip

Deleted selected files
author petr-novak
date Wed, 29 Apr 2020 04:31:31 -0400
parents d08801a298ad
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
12
d08801a298ad Uploaded
petr-novak
parents: 11
diff changeset
1 <tool id="repeatexplorer2" name="RepeatExplorer2 clustering: " version="2.3.7" >
0
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
2 <stdio>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
3 <regex match="lastdb: can't open file: NEAR" source="stderr" level="fatal" description="Version of last is too old, use ver 956 or higher\n" />
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
4 <regex match="Traceback" source="stderr" level="fatal" description="Unknown error" />
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
5 <regex match="error" source="stderr" level="fatal" description="Unknown error" />
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
6 <regex match="Warning" source="stderr" level="warning" description="Unknown error" />
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
7 <exit_code range="1:" level="fatal" description="Error" />
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
8 </stdio>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
9 <description>Improved version or repeat discovery and characterization using graph based sequence clustering</description>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
10 <requirements>
2
394a697ffa49 Uploaded
petr-novak
parents: 0
diff changeset
11 <requirement type="package">last</requirement>
8
48704d3ce3a6 Uploaded
petr-novak
parents: 6
diff changeset
12 <requirement type="package">imagemagick</requirement>
0
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
13 <requirement type="package">mafft</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
14 <requirement type="package">blast</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
15 <requirement type="package">diamond</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
16 <requirement type="package">blast-legacy</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
17 <requirement type="package">r-igraph</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
18 <requirement type="package">r-data.tree</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
19 <requirement type="package">r-stringr</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
20 <requirement type="package">r-r2html</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
21 <requirement type="package">r-hwriter</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
22 <requirement type="package">r-dt</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
23 <requirement type="package">r-scales</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
24 <requirement type="package">r-plotrix</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
25 <requirement type="package">r-png</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
26 <requirement type="package">r-plyr</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
27 <requirement type="package">r-dplyr</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
28 <requirement type="package">r-optparse</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
29 <requirement type="package">r-dbi</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
30 <requirement type="package">r-rsqlite</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
31 <requirement type="package">r-rserve</requirement>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
32 <requirement type="package">bioconductor-biostrings</requirement>
12
d08801a298ad Uploaded
petr-novak
parents: 11
diff changeset
33 <requirement type="package" version="2.3.7">repex_tarean</requirement>
2
394a697ffa49 Uploaded
petr-novak
parents: 0
diff changeset
34 <requirement type="set_environment">REPEX</requirement>
11
3b1ebd937035 Uploaded
petr-novak
parents: 10
diff changeset
35 <requirement type="set_environment">REPEX_VERSION</requirement>
8
48704d3ce3a6 Uploaded
petr-novak
parents: 6
diff changeset
36 <requirement type="package" version="0.9.1" >pyrserve</requirement>
0
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
37 </requirements>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
38 <command >
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
39 export PYTHONHASHSEED=0;
2
394a697ffa49 Uploaded
petr-novak
parents: 0
diff changeset
40 \${REPEX}/seqclust --sample ${sample} --output_dir=tarean_output --logfile=${log} --cleanup $paired --taxon $taxon
0
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
41
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
42 #if $advanced_options.advanced:
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
43 --mincl $advanced_options.size_threshold $advanced_options.keep_names $advanced_options.automatic_filtering -D $advanced_options.blastx.options_blastx
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
44 --assembly_min $advanced_options.assembly_min_cluster_size
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
45
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
46 #if $advanced_options.comparative.options_comparative:
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
47 --prefix_length $advanced_options.comparative.prefix_length
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
48 #end if
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
49
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
50 #if $advanced_options.custom_library.options_custom_library:
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
51 -d $advanced_options.custom_library.library extra_database
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
52 #end if
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
53
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
54 #if $advanced_options.options.options:
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
55 -opt $advanced_options.options.options
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
56 #end if
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
57 #end if
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
58 ${FastaFile} >stdout.log 2> stderr.log ;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
59 echo "STDOUT CONTENT:" >> ${log} ;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
60 cat stdout.log >> ${log} ;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
61 echo "STDERR CONTENT:" >> ${log};
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
62 cat stderr.log >> ${log} &amp;&amp;
10
7d88d9fabe99 Uploaded
petr-novak
parents: 8
diff changeset
63 \${REPEX}/stderr_filter.py stderr.log &amp;&amp;
0
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
64 cd tarean_output &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
65 zip -r ${ReportArchive}.zip * &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
66 mv ${ReportArchive}.zip ${ReportArchive} &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
67 cp index.html ${ReportFile} &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
68 mkdir ${ReportFile.files_path} &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
69 cp -r --parents libdir ${ReportFile.files_path} &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
70 cp -r --parents seqclust/clustering/superclusters ${ReportFile.files_path} &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
71 cp -r --parents seqclust/clustering/clusters ${ReportFile.files_path} &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
72 cp seqclust/clustering/hitsort.cls ${ReportFile.files_path}/seqclust/clustering/hitsort.cls &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
73 cp *.png ${ReportFile.files_path}/ &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
74 cp *.csv ${ReportFile.files_path}/ &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
75 cp *.html ${ReportFile.files_path}/ &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
76 cp *.css ${ReportFile.files_path}/ &amp;&amp;
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
77 cp *.fasta ${ReportFile.files_path}/ 2>>$log &amp;&amp; rm -r ../tarean_output || :
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
78
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
79 </command>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
80 <inputs>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
81 <param name="FastaFile" label="NGS reads" type="data" format="fasta"
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
82 help="Input file must contain fasta-formatted NGS reads. If paired end reads are used, reads must be interlaced and all pairs must be complete. Example of input data format is provided in the help below. "/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
83 <param name="paired" type="boolean" truevalue="--paired" falsevalue="" checked="True" label="Paired-end reads" help="Check if you are using pair reads and input sequences contain both read mates and left mates alternate with their right mates" />
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
84
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
85 <param name="sample" label="Sample size" type="integer" value="500000" min="10000"/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
86 <param name="taxon" label="Select taxon and protein domain database version (REXdb)" type="select" help="Reference database of transposable element protein domains - REXdb - is used for annotation of repeats">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
87 <option value="VIRIDIPLANTAE3.0" selected="true">Viridiplantae version 3.0 </option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
88 <option value="VIRIDIPLANTAE2.2" selected="true">Viridiplantae version 2.2</option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
89 <option value="METAZOA3.0" >Metazoa version 3.0</option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
90 <option value="METAZOA2.0" >Metazoa version 2.0</option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
91 <!-- Modify setting in config.py accordingly -->
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
92 </param>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
93
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
94 <conditional name="advanced_options">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
95 <param name="advanced" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Advanced options" />
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
96 <when value="false">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
97 <!-- pass -->
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
98 </when>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
99 <when value="true">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
100 <conditional name="comparative">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
101 <param name="options_comparative" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Perform comparative analysis" help="Use this options when you want to compare sequences multiple groups"/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
102 <when value="false">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
103 <!-- do nothing here -->
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
104 </when>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
105 <when value="true">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
106 <param name="prefix_length" label="Group code length" type="integer" value="3" min="1" max="10" help="For comparative analysis, sequences are from individial groups distinguished by sample code which must be used as prefix for sequence name. See example below."/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
107 </when>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
108 </conditional>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
109
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
110 <conditional name="blastx">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
111 <param name="options_blastx" type="select" label="Select parameters for protein domain search">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
112 <option value="BLASTX_W2" selected="false">blastx with word size 2 (the most sensitive, slowest)</option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
113 <option value="BLASTX_W3" selected="true">blastx with word size 3 (default)</option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
114 <option value="DIAMOND" selected="false">diamond program (the least sensitive, fastest)</option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
115 </param>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
116 </conditional>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
117
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
118 <conditional name="options">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
119 <param name="options" type="select" label="Similarity search options" help="Different similarity search parameters are used depending on the used input data to adjust search to differences in length and error rate">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
120 <option value="ILLUMINA" selected="true">Illumina reads, read length 100nt or more </option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
121 <option value="ILLUMINA_SHORT" selected="false">Illumina reads, shorter than 100nt (Do not use reads shorter than 50nt!) </option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
122 <option value="ILLUMINA_DUST_OFF" selected="false">Illumina reads, no masking of low complexity repeats </option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
123 <option value="OXFORD_NANOPORE" selected="false">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
124 Pseudo short reads simulated from Oxford Nanopore data (experimental feature)
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
125 </option>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
126 </param>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
127 </conditional>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
128
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
129 <conditional name="custom_library">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
130 <param name="options_custom_library" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Use custom repeat database"/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
131 <when value="false">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
132 <!-- do nothing here -->
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
133 </when>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
134 <when value="true">
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
135 <param name="library" format="fasta" type="data" label="Custom library of repeats" help="Library of repeats as DNA sequences in fasta format. The required format for IDs in a custom library is : '>reapeatname#class/subclass'"/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
136 </when>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
137 </conditional>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
138 <param name="size_threshold" label="Cluster size threshold for detailed analysis" type="float" value="0.01" min="0.0001" max="100" help ="Minimal size (as percentage of input reads) of the smallest cluster which is analyzed, cluster with less than 20 reads are not considered at all."/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
139 <param name="automatic_filtering" label="Perform automatic filtering of abundant satellite repeats" help="Automatic filtering tries to identify the most abundant tandem repeats and remove such sequences partially from analysis. Removal of abundant tandem repeat can enable to analyze higher proportion of other less abundant repeats." type="boolean" truevalue="--automatic_filtering" falsevalue="" checked="false"/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
140 <param name="keep_names" label="Keep original sequences names" type="boolean" truevalue="--keep_names" falsevalue="" checked="false" help="By default sequence are relabeled using integers. If you want to keep original names, use this option."/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
141 <param name="assembly_min_cluster_size" type="integer" label="min cluster size for assembly" value="5" min="2" max="100"/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
142 </when>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
143 </conditional>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
144
12
d08801a298ad Uploaded
petr-novak
parents: 11
diff changeset
145
d08801a298ad Uploaded
petr-novak
parents: 11
diff changeset
146
0
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
147 </inputs>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
148 <outputs>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
149 <data name="log" format="txt" label="RepeatExplorer2 - log file"/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
150 <data name="ReportArchive" format="zip" label="RepeatExplorer2 - Archive with HTML report from data ${FastaFile.hid}"/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
151 <data name="ReportFile" format="html" label="RepeatExplorer2 - HTML report from data ${FastaFile.hid}"/>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
152 </outputs>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
153
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
154 <help>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
155 **HELP**
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
156
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
157 RepeatExplorer2 clustering is a computational pipeline for unsupervised
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
158 identification of repeats from unassembled sequence reads. The
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
159 pipeline uses low-pass whole genome sequence reads and performs graph-based
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
160 clustering. Resulting clusters, representing all types of repeats, are then
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
161 examined to identify and classify into repeats groups.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
162
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
163 **Input data**
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
164
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
165 The analysis requires either **single** or **paired-end reads** generated
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
166 by whole genome shotgun sequencing provided as a single fasta-formatted file.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
167 Generally, paired-end reads provide significantly better results than single
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
168 reads. Reads should be of uniform length (optimal size range is 100-200 nt) and
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
169 the number of analyzed reads should represent less than 1x genome equivalent
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
170 (genome coverage of 0.01 - 0.50 x is recommended). Reads should be
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
171 quality-filtered (recommended filtering : quality score >=10 over 95% of bases
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
172 and no Ns allowed) and only **complete read pairs** should be submitted for
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
173 analysis. When paired reads are used, input data must be **interlaced** format
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
174 as fasta file:
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
175
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
176 example of interlaced input format::
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
177
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
178 >0001_f
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
179 CGTAATATACATACTTGCTAGCTAGTTGGATGCATCCAACTTGCAAGCTAGTTTGATG
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
180 >0001_r
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
181 GATTTGACGGACACACTAACTAGCTAGTTGCATCTAAGCGGGCACACTAACTAACTAT
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
182 >0002_f
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
183 ACTCATTTGGACTTAACTTTGATAATAAAAACTTAAAAAGGTTTCTGCACATGAATCG
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
184 >0002_r
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
185 TATGTTGAAAAATTGAATTTCGGGACGAAACAGCGTCTATCGTCACGACATAGTGCTC
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
186 >0003_f
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
187 TGACATTTGTGAACGTTAATGTTCAACAAATCTTTCCAATGTCTTTTTATCTTATCAT
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
188 >0003_r
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
189 TATTGAAATACTGGACACAAATTGGAAATGAAACCTTGTGAGTTATTCAATTTATGTT
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
190 ...
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
191
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
192
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
193 **Comparative analysis**
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
194
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
195 For comparative analysis sequence names must contain code (prefix) for each group.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
196 Prefix in sequences names must be of fixed length.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
197
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
198 Example of labeling two groups with where **group code length** is 2 and is used to distinguish groups - AA and BB ::
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
199
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
200 >AA0001_f
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
201 CGTAATATACATACTTGCTAGCTAGTTGGATGCATCCAACTTGCAAGCTAGTTTGATG
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
202 >AA0001_r
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
203 GATTTGACGGACACACTAACTAGCTAGTTGCATCTAAGCGGGCACACTAACTAACTAT
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
204 >AA0002_f
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
205 ACTCATTTGGACTTAACTTTGATAATAAAAACTTAAAAAGGTTTCTGCACATGAATCG
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
206 >AA0002_r
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
207 TATGTTGAAAAATTGAATTTCGGGACGAAACAGCGTCTATCGTCACGACATAGTGCTC
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
208 >BB0001_f
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
209 TGACATTTGTGAACGTTAATGTTCAACAAATCTTTCCAATGTCTTTTTATCTTATCAT
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
210 >BB0001_r
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
211 TATTGAAATACTGGACACAAATTGGAAATGAAACCTTGTGAGTTATTCAATTTATGTT
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
212 >BB0002_f
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
213 TGACATTTGTGAACGTTAATGTTCAACAAATCTTTCCAATGTCTTTTTATCTTATCAT
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
214 >BB0002_r
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
215 TATTGAAATACTGGACACAAATTGGAAATGAAACCTTGTGAGTTATTCAATTTATGTT
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
216
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
217
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
218 To prepare quality filtered and interlaced input fasta file from fastq
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
219 files, use `Preprocessing of paired-reads`__ tool.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
220
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
221 .. __: tool_runner?tool_id=paired_fastq_filtering
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
222
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
223
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
224 **Additional parameters**
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
225
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
226 **Sample size** defines how many reads should be used in calculation.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
227 Default setting with 500,000 reads will enable detection of high copy
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
228 repeats within several hours of computation time. For higher
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
229 sensitivity the sample size can be set higher. Since sample size affects
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
230 the memory usage, this parameter may be automatically adjusted to lower
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
231 value during the run. Maximum sample size which can be processed depends on
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
232 the repetitiveness of analyzed genome.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
233
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
234
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
235 **Select taxon and protein domain database version (REXdb)**. Classification
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
236 of transposable elements is based on the similarity to our reference database
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
237 of transposable element protein domains (**REXdb**). Standalone database for Viridiplantae species
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
238 can be obtained on `repeatexplorer.org`__. Classification
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
239 system used in REXdb is described in article `Systematic survey of plant
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
240 LTR-retrotransposons elucidates phylogenetic relationships of their
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
241 polyprotein domains and provides a reference for element classification`__
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
242 Database for Metazoa species is still under development so use it with caution.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
243
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
244 .. __: http://repeatexplorer.org
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
245 .. __: https://doi.org/10.1186/s13100-018-0144-1
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
246
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
247 **Select parameters for protein domain search** REXdb is compared with s
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
248 equence clusters either using blastx or diamond aligner. Diamond program
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
249 is about three time faster than blastx with word size 3.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
250
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
251 **Similarity search options** By default sequence reads are compared using
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
252 mgblast program. Default threshold is explicitly set to 90% sequence
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
253 similarity spanning at least 55% of the read length (in the case of reads
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
254 differing in length it applies to the longer one). Additionally, sequence
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
255 overlap must be at least 55 nt. If you select option for shorter reads
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
256 than 100 nt, minimum overlap 55 nt is not required.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
257
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
258 By default,
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
259 mgblast search use DUST program to filter out
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
260 low-complexity sequences. If you want
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
261 to increase sensitivity of detection of satellites with shorter monomer
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
262 use option with '*no masking of low complexity repeats*'. Note that omitting
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
263 DUST filtering will significantly increase running times
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
264
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
265
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
266 **Automatic filtering of abundant satellite repeats** perform clustering on
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
267 smaller dataset of sequence reads to detect abundant high confidence
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
268 satellite repeats. If such satellites are detected, sequence reads derived
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
269 from these satellites are depleted from input dataset. This step enable more
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
270 sensitive detection of less abundant repeats as more reads can be used
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
271 in clustering step.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
272
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
273 **Use custom repeat database**. This option allows users to perform similarity
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
274 comparison of identified repeats to their custom databases. The repeat class must
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
275 be encoded in FASTA headers of database entries in order to allow correct
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
276 parsing of similarity hits. Required format for custom database sequence name is: ::
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
277
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
278 >reapeatname#class/subclass
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
279
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
280
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
281 **Output**
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
282
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
283 List of clusters identified as putative satellite repeats, their genomic
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
284 abundance and various cluster characteristics.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
285
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
286 Output includes a **HTML summary** with table listing of all analyzed
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
287 clusters. More detailed information about clusters is provided in
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
288 additional files and directories. All results are also provided as
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
289 downloadable **zip archive**. Additionally a **log file** reporting
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
290 the progress of the computational pipeline is provided.
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
291
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
292 </help>
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
293
1d1b9e1b2e2f Uploaded
petr-novak
parents:
diff changeset
294 </tool>