annotate repex_tarean.xml @ 25:a27646c27090 draft default tip

Uploaded
author petr-novak
date Wed, 05 Jun 2024 08:47:44 +0000
parents 5d6ff2684c5b
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
25
a27646c27090 Uploaded
petr-novak
parents: 23
diff changeset
1 <tool id="tarean" name="Tandem Repeat Analyzer (TAREAN)" version="2.3.11.1" >
2
349b197133dc Uploaded
petr-novak
parents:
diff changeset
2 <stdio>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
3 <regex match="Traceback" source="stderr" level="fatal" description="Unknown error" />
349b197133dc Uploaded
petr-novak
parents:
diff changeset
4 <regex match="error" source="stderr" level="fatal" description="Unknown error" />
349b197133dc Uploaded
petr-novak
parents:
diff changeset
5 <regex match="warning" source="stderr" level="warning" description="Unknown warning" />
349b197133dc Uploaded
petr-novak
parents:
diff changeset
6 <exit_code range="1:" level="fatal" description="Error" />
349b197133dc Uploaded
petr-novak
parents:
diff changeset
7 </stdio>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
8 <description>Identification of genomic tandem repeats from NGS data</description>
4
72c028e3e70e Uploaded
petr-novak
parents: 3
diff changeset
9 <requirements>
25
a27646c27090 Uploaded
petr-novak
parents: 23
diff changeset
10 <container type="singularity">library://repeatexplorer/default/repex_tarean:0.3.11-579a65d</container>
4
72c028e3e70e Uploaded
petr-novak
parents: 3
diff changeset
11 </requirements>
2
349b197133dc Uploaded
petr-novak
parents:
diff changeset
12 <command detect_errors="exit_code">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
13 export PYTHONHASHSEED=0;
16
20ef42e81eeb Uploaded
petr-novak
parents: 12
diff changeset
14 export TAREAN_CPU=\$GALAXY_SLOTS;
19
3d0ba5734776 Uploaded
petr-novak
parents: 16
diff changeset
15 export TAREAN_MAX_MEM=\$((GALAXY_MEMORY_MB*1024));
4
72c028e3e70e Uploaded
petr-novak
parents: 3
diff changeset
16 seqclust --paired --sample ${read_sampling.sample} --output_dir=tarean_output --logfile=${log} --cleanup --tarean_mode
2
349b197133dc Uploaded
petr-novak
parents:
diff changeset
17 #if $advanced_options.advanced:
349b197133dc Uploaded
petr-novak
parents:
diff changeset
18 --mincl $advanced_options.size_threshold $advanced_options.keep_names $advanced_options.automatic_filtering -M $advanced_options.merging
349b197133dc Uploaded
petr-novak
parents:
diff changeset
19 #if $advanced_options.custom_library.options_custom_library :
349b197133dc Uploaded
petr-novak
parents:
diff changeset
20 -d $advanced_options.custom_library.library extra_database
349b197133dc Uploaded
petr-novak
parents:
diff changeset
21 #end if
349b197133dc Uploaded
petr-novak
parents:
diff changeset
22 #if $advanced_options.options.options:
349b197133dc Uploaded
petr-novak
parents:
diff changeset
23 -opt $advanced_options.options.options
349b197133dc Uploaded
petr-novak
parents:
diff changeset
24 #end if
349b197133dc Uploaded
petr-novak
parents:
diff changeset
25 #else:
349b197133dc Uploaded
petr-novak
parents:
diff changeset
26 -M 0.2
349b197133dc Uploaded
petr-novak
parents:
diff changeset
27
349b197133dc Uploaded
petr-novak
parents:
diff changeset
28 #end if
349b197133dc Uploaded
petr-novak
parents:
diff changeset
29 ${FastaFile} >stdout.log 2> stderr.log ;
349b197133dc Uploaded
petr-novak
parents:
diff changeset
30 echo "STDOUT CONTENT:" >> ${log} ;
349b197133dc Uploaded
petr-novak
parents:
diff changeset
31 cat stdout.log >> ${log} ;
349b197133dc Uploaded
petr-novak
parents:
diff changeset
32 echo "STDERR CONTENT:" >> ${log} ;
349b197133dc Uploaded
petr-novak
parents:
diff changeset
33 cat stderr.log >> ${log} &amp;&amp;
4
72c028e3e70e Uploaded
petr-novak
parents: 3
diff changeset
34 /opt/repex_tarean/stderr_filter.py stderr.log &amp;&amp;
2
349b197133dc Uploaded
petr-novak
parents:
diff changeset
35 cd tarean_output &amp;&amp;
23
5d6ff2684c5b Uploaded
petr-novak
parents: 21
diff changeset
36 zip -fz -r ${ReportArchive}.zip * &amp;&amp;
2
349b197133dc Uploaded
petr-novak
parents:
diff changeset
37 mv ${ReportArchive}.zip ${ReportArchive} &amp;&amp;
349b197133dc Uploaded
petr-novak
parents:
diff changeset
38 cp index.html ${ReportFile} &amp;&amp;
11
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
39 mkdir -p ${ReportFile.extra_files_path} &amp;&amp;
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
40 cp -r --parents libdir ${ReportFile.extra_files_path} &amp;&amp;
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
41 cp -r --parents seqclust/clustering/superclusters ${ReportFile.extra_files_path} &amp;&amp;
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
42 cp -r --parents seqclust/clustering/clusters ${ReportFile.extra_files_path} &amp;&amp;
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
43 cp seqclust/clustering/hitsort.cls ${ReportFile.extra_files_path}/seqclust/clustering/hitsort.cls &amp;&amp;
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
44 cp *.png ${ReportFile.extra_files_path}/ &amp;&amp;
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
45 cp *.csv ${ReportFile.extra_files_path}/ &amp;&amp;
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
46 cp *.html ${ReportFile.extra_files_path}/ &amp;&amp;
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
47 cp *.css ${ReportFile.extra_files_path}/ &amp;&amp;
1d1dc045695e Uploaded
petr-novak
parents: 7
diff changeset
48 cp *.fasta ${ReportFile.extra_files_path}/ 2>>$log &amp;&amp; rm -r ../tarean_output || :
2
349b197133dc Uploaded
petr-novak
parents:
diff changeset
49
349b197133dc Uploaded
petr-novak
parents:
diff changeset
50
349b197133dc Uploaded
petr-novak
parents:
diff changeset
51 </command>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
52
349b197133dc Uploaded
petr-novak
parents:
diff changeset
53 <inputs>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
54 <param name="FastaFile" label="Paired-end Illumina reads" type="data" format="fasta"
349b197133dc Uploaded
petr-novak
parents:
diff changeset
55 help="Input file must contain FASTA-formatted interlaced read pairs from paired-end sequencing. All pairs must be complete. Example of the input data format is provided in the help below."/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
56
349b197133dc Uploaded
petr-novak
parents:
diff changeset
57 <conditional name="read_sampling">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
58 <param name="do_sampling" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Read sampling" help="Use this option if you want to analyze only a part of the reads" />
349b197133dc Uploaded
petr-novak
parents:
diff changeset
59 <when value="false">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
60 <!-- pass -->
16
20ef42e81eeb Uploaded
petr-novak
parents: 12
diff changeset
61 <param name="sample" label="Sample size" type="hidden" value="0" help="Number of analyzed reads"/>
2
349b197133dc Uploaded
petr-novak
parents:
diff changeset
62 </when>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
63 <when value="true">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
64 <param name="sample" label="Sample size" type="integer" value="500000" min="10000" help="Number of analyzed reads"/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
65 </when>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
66 </conditional>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
67
349b197133dc Uploaded
petr-novak
parents:
diff changeset
68 <conditional name="advanced_options">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
69 <param name="advanced" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Advanced options" />
349b197133dc Uploaded
petr-novak
parents:
diff changeset
70 <when value="false">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
71 <!-- pass -->
349b197133dc Uploaded
petr-novak
parents:
diff changeset
72 </when>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
73 <when value="true">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
74 <param name="merging" type="boolean" truevalue="0.2" falsevalue="0" checked="True" label="Perform cluster merging" help="By default, clusters connected through paired-end reads are merged"/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
75 <conditional name="custom_library">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
76 <param name="options_custom_library" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Use custom repeat database"/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
77 <when value="false">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
78 <!-- do nothing here -->
349b197133dc Uploaded
petr-novak
parents:
diff changeset
79 </when>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
80 <when value="true">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
81 <param name="library" format="fasta" type="data" label="Use custom repeat database" help="Perform additional similarity search to user-provided repeat database. The database should contain FASTA-formatted DNA sequences with headers (sequence names) in the format: '>reapeatname#class/subclass'"/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
82 </when>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
83 </conditional>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
84 <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."/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
85 <param name="automatic_filtering" label="Perform automatic filtering of abundant satellite repeats" type="boolean" truevalue="--automatic_filtering" falsevalue="" checked="false"/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
86 <param name="keep_names" label="Keep original read names" type="boolean" truevalue="--keep_names" falsevalue="" checked="false" help="By default, reads are renamed using integers. Use this option if you want to keep original names."/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
87 <conditional name="options">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
88 <param name="options" type="select" label="Similarity search options">
349b197133dc Uploaded
petr-novak
parents:
diff changeset
89 <option value="ILLUMINA" selected="true">Default </option>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
90 <option value="ILLUMINA_DUST_OFF" selected="false">Masking of low complexity repeats disabled </option>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
91
349b197133dc Uploaded
petr-novak
parents:
diff changeset
92 <!-- <option value="ILLUMINA_SENSITIVE_MGBLAST" selected="false">Illumina reads, sensitive search (search parameters: mgblast, min PID 80, -W8) slow, experimental feature!</option> -->
349b197133dc Uploaded
petr-novak
parents:
diff changeset
93 <!-- <option value="ILLUMINA_SENSITIVE_BLASTPLUS" selected="false">Illumina reads, more sensitive search (search parameters: blastn, min PID 80, -W6) extremely slow, experimental feature!</option> -->
349b197133dc Uploaded
petr-novak
parents:
diff changeset
94 <!-- <option value="OXFORD_NANOPORE" selected="false"> -->
349b197133dc Uploaded
petr-novak
parents:
diff changeset
95 <!-- Pseudo short reads simulated from Oxford Nanopore data, experimental feature! -->
349b197133dc Uploaded
petr-novak
parents:
diff changeset
96 <!-- </option> -->
349b197133dc Uploaded
petr-novak
parents:
diff changeset
97 </param>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
98 </conditional>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
99 </when>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
100 </conditional>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
101
16
20ef42e81eeb Uploaded
petr-novak
parents: 12
diff changeset
102 <param name="queue_select" type="select" label="Select queue">
20ef42e81eeb Uploaded
petr-novak
parents: 12
diff changeset
103 <option value="basic_fast_queue">basic (max runtime 2 days, 4 GB RAM)</option>
20ef42e81eeb Uploaded
petr-novak
parents: 12
diff changeset
104 <option value="long_slow_queue">long (max runtime 2 weeks, 64 GB RAM)</option>
20ef42e81eeb Uploaded
petr-novak
parents: 12
diff changeset
105 <option value="extra_long_slow_queue">extra long (max runtime 4 weeks, 64 GB RAM)</option>
20ef42e81eeb Uploaded
petr-novak
parents: 12
diff changeset
106 </param>
2
349b197133dc Uploaded
petr-novak
parents:
diff changeset
107
349b197133dc Uploaded
petr-novak
parents:
diff changeset
108 </inputs>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
109 <outputs>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
110 <data name="log" format="txt" label="TAREAN log file"/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
111 <data name="ReportArchive" format="zip" label="TAREAN Archive with HTML report from data ${FastaFile.hid}"/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
112 <data name="ReportFile" format="html" label="TAREAN HTML report from data ${FastaFile.hid}"/>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
113 </outputs>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
114
349b197133dc Uploaded
petr-novak
parents:
diff changeset
115 <help>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
116 **HELP**
349b197133dc Uploaded
petr-novak
parents:
diff changeset
117
349b197133dc Uploaded
petr-novak
parents:
diff changeset
118 TAREAN - TAndem REpeat ANalyzer is a computational pipeline for
349b197133dc Uploaded
petr-novak
parents:
diff changeset
119 **unsupervised identification of satellite repeats** from unassembled
349b197133dc Uploaded
petr-novak
parents:
diff changeset
120 sequence reads. The pipeline uses low-pass paired-end whole genome
349b197133dc Uploaded
petr-novak
parents:
diff changeset
121 sequence reads and performs graph-based clustering. The resulting
349b197133dc Uploaded
petr-novak
parents:
diff changeset
122 clusters, representing all types of repeats present in the genome, are
349b197133dc Uploaded
petr-novak
parents:
diff changeset
123 then examined to identify those containing circular structures indicative
349b197133dc Uploaded
petr-novak
parents:
diff changeset
124 of tandem repeats. A poster summarizing TAREAN principles and
349b197133dc Uploaded
petr-novak
parents:
diff changeset
125 implementation can be found `here.`__
349b197133dc Uploaded
petr-novak
parents:
diff changeset
126
349b197133dc Uploaded
petr-novak
parents:
diff changeset
127
349b197133dc Uploaded
petr-novak
parents:
diff changeset
128 .. __: http://w3lamc.umbr.cas.cz/lamc/?page_id=312
349b197133dc Uploaded
petr-novak
parents:
diff changeset
129
349b197133dc Uploaded
petr-novak
parents:
diff changeset
130 **Input data**
349b197133dc Uploaded
petr-novak
parents:
diff changeset
131
349b197133dc Uploaded
petr-novak
parents:
diff changeset
132
349b197133dc Uploaded
petr-novak
parents:
diff changeset
133 The analysis requires **paired-end reads** generated by whole genome
349b197133dc Uploaded
petr-novak
parents:
diff changeset
134 shotgun sequencing. The data should be provided as a single input file in
349b197133dc Uploaded
petr-novak
parents:
diff changeset
135 fasta format with the reads interlaced (see example below). All the pairs
349b197133dc Uploaded
petr-novak
parents:
diff changeset
136 must be complete, i.e. both "forward" and "reverse" sequence reads must be
349b197133dc Uploaded
petr-novak
parents:
diff changeset
137 present. The reads should all be trimmed to the same length. The optimal
349b197133dc Uploaded
petr-novak
parents:
diff changeset
138 size range is between 100 and 200 nucleotides. The number of reads to be
349b197133dc Uploaded
petr-novak
parents:
diff changeset
139 analyzed should not exceed 1x coverage of the genome. Genome coverage
349b197133dc Uploaded
petr-novak
parents:
diff changeset
140 between 0.01 and 0.5x is recommended. The reads should be filtered for
349b197133dc Uploaded
petr-novak
parents:
diff changeset
141 quality. The recommended quality filtering is as follows: each read should
349b197133dc Uploaded
petr-novak
parents:
diff changeset
142 have a quality score >=10 for 95% of the bases, i.e. if your reads are 100
349b197133dc Uploaded
petr-novak
parents:
diff changeset
143 base pairs long, then a read only passes this quality threshold if 95
349b197133dc Uploaded
petr-novak
parents:
diff changeset
144 bases have a quality of 10 or higher. Additionally, any reads containing
349b197133dc Uploaded
petr-novak
parents:
diff changeset
145 indeterminate base pairs (indicated as N in the reads) should be removed.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
146 Finally, if either one of the reads in a pair fails to meet the
349b197133dc Uploaded
petr-novak
parents:
diff changeset
147 aforementioned thresholds, **both** sequences should be removed.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
148 example of interlaced input format::
349b197133dc Uploaded
petr-novak
parents:
diff changeset
149
349b197133dc Uploaded
petr-novak
parents:
diff changeset
150 >0001_f
349b197133dc Uploaded
petr-novak
parents:
diff changeset
151 CGTAATATACATACTTGCTAGCTAGTTGGATGCATCCAACTTGCAAGCTAGTTTGATG
349b197133dc Uploaded
petr-novak
parents:
diff changeset
152 >0001_r
349b197133dc Uploaded
petr-novak
parents:
diff changeset
153 GATTTGACGGACACACTAACTAGCTAGTTGCATCTAAGCGGGCACACTAACTAACTAT
349b197133dc Uploaded
petr-novak
parents:
diff changeset
154 >0002_f
349b197133dc Uploaded
petr-novak
parents:
diff changeset
155 ACTCATTTGGACTTAACTTTGATAATAAAAACTTAAAAAGGTTTCTGCACATGAATCG
349b197133dc Uploaded
petr-novak
parents:
diff changeset
156 >0002_r
349b197133dc Uploaded
petr-novak
parents:
diff changeset
157 TATGTTGAAAAATTGAATTTCGGGACGAAACAGCGTCTATCGTCACGACATAGTGCTC
349b197133dc Uploaded
petr-novak
parents:
diff changeset
158 >0003_f
349b197133dc Uploaded
petr-novak
parents:
diff changeset
159 TGACATTTGTGAACGTTAATGTTCAACAAATCTTTCCAATGTCTTTTTATCTTATCAT
349b197133dc Uploaded
petr-novak
parents:
diff changeset
160 >0003_r
349b197133dc Uploaded
petr-novak
parents:
diff changeset
161 TATTGAAATACTGGACACAAATTGGAAATGAAACCTTGTGAGTTATTCAATTTATGTT
349b197133dc Uploaded
petr-novak
parents:
diff changeset
162 ...
349b197133dc Uploaded
petr-novak
parents:
diff changeset
163
349b197133dc Uploaded
petr-novak
parents:
diff changeset
164
349b197133dc Uploaded
petr-novak
parents:
diff changeset
165 To perform the quality filtering on your fastQ formatted data as described
349b197133dc Uploaded
petr-novak
parents:
diff changeset
166 above, and to interlace your paired-end sequence reads,
349b197133dc Uploaded
petr-novak
parents:
diff changeset
167 please use the `Preprocessing of paired-reads`__ tool.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
168
349b197133dc Uploaded
petr-novak
parents:
diff changeset
169 .. __: tool_runner?tool_id=paired_fastq_filtering
349b197133dc Uploaded
petr-novak
parents:
diff changeset
170
349b197133dc Uploaded
petr-novak
parents:
diff changeset
171
349b197133dc Uploaded
petr-novak
parents:
diff changeset
172 **Additional parameters**
349b197133dc Uploaded
petr-novak
parents:
diff changeset
173
349b197133dc Uploaded
petr-novak
parents:
diff changeset
174 **Sample size** defines how many reads will be used during the computation.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
175 The default setting of 500,000 reads will enable detection of high copy
349b197133dc Uploaded
petr-novak
parents:
diff changeset
176 number satellites within several hours. For higher
349b197133dc Uploaded
petr-novak
parents:
diff changeset
177 sensitivity the sample size can be increased. Since the sample size affects
349b197133dc Uploaded
petr-novak
parents:
diff changeset
178 memory usage, this parameter may be automatically adjusted to a lower value
349b197133dc Uploaded
petr-novak
parents:
diff changeset
179 during the run. The maximum sample size which can be processed depends on the
349b197133dc Uploaded
petr-novak
parents:
diff changeset
180 repetitiveness of the analyzed genome. This significantly limits the number of reads
349b197133dc Uploaded
petr-novak
parents:
diff changeset
181 that can be analyzed with the TAREAN pipeline.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
182
349b197133dc Uploaded
petr-novak
parents:
diff changeset
183 **Perform cluster merging**. Families of repetitive elements are
349b197133dc Uploaded
petr-novak
parents:
diff changeset
184 frequently split into multiple clusters rather than being represented as a
349b197133dc Uploaded
petr-novak
parents:
diff changeset
185 single one. If you do not want to merge clusters based on the presence
349b197133dc Uploaded
petr-novak
parents:
diff changeset
186 of broken read pairs, disable this option.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
187
349b197133dc Uploaded
petr-novak
parents:
diff changeset
188 **Use custom repeat database**. This option allows users to perform similarity
349b197133dc Uploaded
petr-novak
parents:
diff changeset
189 comparison of identified repeats to their custom databases. The repeat class should
349b197133dc Uploaded
petr-novak
parents:
diff changeset
190 be encoded in FASTA headers of database entries in order to allow correct
349b197133dc Uploaded
petr-novak
parents:
diff changeset
191 parsing of similarity hits.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
192
349b197133dc Uploaded
petr-novak
parents:
diff changeset
193 **Similarity search options** By default sequence reads are compared using
349b197133dc Uploaded
petr-novak
parents:
diff changeset
194 mgblast program. Default threshold is explicitly set to 90% sequence
349b197133dc Uploaded
petr-novak
parents:
diff changeset
195 similarity spanning at least 55% of the read length (in the case of reads
349b197133dc Uploaded
petr-novak
parents:
diff changeset
196 differing in length it applies to the longer one). Additionally, sequence
349b197133dc Uploaded
petr-novak
parents:
diff changeset
197 overlap must be at least 55 nt. If you select option for shorter reads
349b197133dc Uploaded
petr-novak
parents:
diff changeset
198 than 100 nt, minimum overlap 55 nt is not required.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
199
349b197133dc Uploaded
petr-novak
parents:
diff changeset
200 By default,
349b197133dc Uploaded
petr-novak
parents:
diff changeset
201 mgblast search use DUST program to filter out
349b197133dc Uploaded
petr-novak
parents:
diff changeset
202 low-complexity sequences. If you want
349b197133dc Uploaded
petr-novak
parents:
diff changeset
203 to increase sensitivity of detection of satellites with shorter monomer
349b197133dc Uploaded
petr-novak
parents:
diff changeset
204 use option with '*no masking of low complexity repeats*'. Note that omitting
349b197133dc Uploaded
petr-novak
parents:
diff changeset
205 DUST filtering will significantly increase running times
349b197133dc Uploaded
petr-novak
parents:
diff changeset
206
349b197133dc Uploaded
petr-novak
parents:
diff changeset
207 **Output**
349b197133dc Uploaded
petr-novak
parents:
diff changeset
208
349b197133dc Uploaded
petr-novak
parents:
diff changeset
209 A list of clusters identified as putative satellite repeats, their genomic
349b197133dc Uploaded
petr-novak
parents:
diff changeset
210 abundance and various cluster characteristics are provided. Length and
349b197133dc Uploaded
petr-novak
parents:
diff changeset
211 consensus sequences of reconstructed monomers are also shown and
349b197133dc Uploaded
petr-novak
parents:
diff changeset
212 accompanied by a detailed output from kmer-based reconstruction including
349b197133dc Uploaded
petr-novak
parents:
diff changeset
213 sequences and sequence logos of alternative variants of monomer sequences.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
214
349b197133dc Uploaded
petr-novak
parents:
diff changeset
215 The output includes an **HTML summary** with a table listing all analyzed
349b197133dc Uploaded
petr-novak
parents:
diff changeset
216 clusters. More detailed information about clusters is provided in
349b197133dc Uploaded
petr-novak
parents:
diff changeset
217 additional files and directories. All results are also provided as a
349b197133dc Uploaded
petr-novak
parents:
diff changeset
218 downloadable **zip archive**. Since read clustering results in
349b197133dc Uploaded
petr-novak
parents:
diff changeset
219 thousands of clusters, the search for satellite repeats is limited to
349b197133dc Uploaded
petr-novak
parents:
diff changeset
220 a subset of the largest ones corresponding to the most abundant genomic
349b197133dc Uploaded
petr-novak
parents:
diff changeset
221 repeats. The default setting of the pipeline is to analyze all clusters containing at least
349b197133dc Uploaded
petr-novak
parents:
diff changeset
222 0.01% of the input reads. Besides the satellite repeats, three other
349b197133dc Uploaded
petr-novak
parents:
diff changeset
223 groups of clusters are reported in the output (1) LTR-retrotransposons,
349b197133dc Uploaded
petr-novak
parents:
diff changeset
224 (2) 45S and 5S rDNA and (3) all remaining clusters passing the size
349b197133dc Uploaded
petr-novak
parents:
diff changeset
225 threshold. As (1) and (2) contain sequences with circular
349b197133dc Uploaded
petr-novak
parents:
diff changeset
226 graphs, their consensus is calculated in the same way as for satellite
349b197133dc Uploaded
petr-novak
parents:
diff changeset
227 repeats. Additionally a **log file** reporting the progress of the
349b197133dc Uploaded
petr-novak
parents:
diff changeset
228 computational pipeline is provided.
349b197133dc Uploaded
petr-novak
parents:
diff changeset
229
349b197133dc Uploaded
petr-novak
parents:
diff changeset
230
349b197133dc Uploaded
petr-novak
parents:
diff changeset
231 </help>
349b197133dc Uploaded
petr-novak
parents:
diff changeset
232
349b197133dc Uploaded
petr-novak
parents:
diff changeset
233 </tool>