annotate repex_tarean.xml @ 0:e2b8e71b85b9 draft

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