annotate tophat2_aurora.xml @ 0:589de4bdaca4 default tip

for Apache Aurora jobrunner
author Makoto Nasuno <nasuno@ascade.co.jp>
date Fri, 26 Feb 2016 14:23:22 +0900
parents
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
1 <tool id="tophat2_parallel" name="Tophat parallel for aurora" version="0.7">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
2 <!-- Wrapper compatible with Tophat version 2.0.0+ -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
3 <description>Gapped-read mapper for RNA-seq data</description>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
4 <!-- <version_command>tophat2 --version</version_command> -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
5 <requirements>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
6 <container type="docker">nasuno/tophat2:2.0.9</container>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
7 </requirements>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
8
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
9 <command>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
10 ##
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
11 ## Set path to index, building the reference if necessary.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
12 ##
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
13
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
14 #set index_path = ''
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
15 #if $refGenomeSource.genomeSource == "history":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
16 bowtie2-build "$refGenomeSource.ownFile" genome ; ln -s "$refGenomeSource.ownFile" genome.fa ;
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
17 #set index_path = 'genome'
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
18 #else:
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
19 #set index_path = $refGenomeSource.index.fields.path
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
20 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
21
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
22 ##
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
23 ## Run tophat.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
24 ##
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
25
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
26 tophat2
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
27
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
28 ## Change this to accommodate the number of threads you have available.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
29 --num-threads \${GALAXY_SLOTS:-4}
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
30
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
31 --output-dir ./tophat_out
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
32
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
33 ## Set params.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
34 #if $params.settingsType == "full":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
35 --read-mismatches $params.read_mismatches
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
36 #if str($params.bowtie_n) == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
37 --bowtie-n
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
38 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
39
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
40 --read-edit-dist $params.read_edit_dist
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
41 --read-realign-edit-dist $params.read_realign_edit_dist
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
42 -a $params.anchor_length
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
43 -m $params.splice_mismatches
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
44 -i $params.min_intron_length
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
45 -I $params.max_intron_length
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
46 -g $params.max_multihits
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
47 --min-segment-intron $params.min_segment_intron
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
48 --max-segment-intron $params.max_segment_intron
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
49 --segment-mismatches $params.seg_mismatches
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
50 --segment-length $params.seg_length
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
51 --library-type $params.library_type
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
52
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
53 ## Indel search.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
54 #if $params.indel_search.allow_indel_search == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
55 ## --allow-indels
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
56 --max-insertion-length $params.indel_search.max_insertion_length
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
57 --max-deletion-length $params.indel_search.max_deletion_length
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
58 #else:
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
59 --no-novel-indels
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
60 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
61
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
62 ## Supplying junctions parameters.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
63 #if $params.own_junctions.use_junctions == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
64 #if $params.own_junctions.gene_model_ann.use_annotations == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
65 -G $params.own_junctions.gene_model_ann.gene_annotation_model
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
66 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
67 #if $params.own_junctions.raw_juncs.use_juncs == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
68 -j $params.own_junctions.raw_juncs.raw_juncs
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
69 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
70 #if str($params.own_junctions.no_novel_juncs) == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
71 --no-novel-juncs
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
72 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
73 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
74
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
75 #if $params.coverage_search.use_search == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
76 --coverage-search
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
77 --min-coverage-intron $params.coverage_search.min_coverage_intron
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
78 --max-coverage-intron $params.coverage_search.max_coverage_intron
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
79 #else:
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
80 --no-coverage-search
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
81 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
82
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
83 #if str($params.microexon_search) == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
84 --microexon-search
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
85 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
86
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
87 #if $params.fusion_search.do_search == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
88 --fusion-search
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
89 --fusion-anchor-length $params.fusion_search.anchor_len
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
90 --fusion-min-dist $params.fusion_search.min_dist
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
91 --fusion-read-mismatches $params.fusion_search.read_mismatches
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
92 --fusion-multireads $params.fusion_search.multireads
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
93 --fusion-multipairs $params.fusion_search.multipairs
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
94 --fusion-ignore-chromosomes "$params.fusion_search.ignore_chromosomes"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
95 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
96
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
97 #if $params.bowtie2_settings.b2_settings == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
98 #if $params.bowtie2_settings.preset.b2_preset == "Yes":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
99 --b2-$params.bowtie2_settings.preset.b2_preset_select
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
100 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
101 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
102
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
103 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
104
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
105 ## Read group information.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
106 #if $readGroup.specReadGroup == "yes"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
107 --rg-id "$readGroup.rgid"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
108 --rg-library "$readGroup.rglb"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
109 --rg-platform "$readGroup.rgpl"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
110 --rg-sample "$readGroup.rgsm"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
111 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
112
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
113 ## Set index path, inputs and parameters specific to paired data.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
114 #if $singlePaired.sPaired != "single"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
115 -r $singlePaired.mate_inner_distance
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
116 --mate-std-dev=$singlePaired.mate_std_dev
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
117
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
118 #if str($singlePaired.report_discordant_pairs) == "No":
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
119 --no-discordant
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
120 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
121
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
122 #if $singlePaired.sPaired == "paired"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
123 ${index_path} "$singlePaired.input1" "$singlePaired.input2"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
124 #else
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
125 ${index_path} "$singlePaired.input.forward" "$singlePaired.input.reverse"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
126 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
127 #else
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
128 ${index_path} "$singlePaired.input1"
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
129 #end if
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
130 </command>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
131
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
132 <inputs>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
133 <param name="job_annotation" type="text" value="tophat.parallel" label="Job Annotation" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
134 <conditional name="singlePaired">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
135 <param name="sPaired" type="select" label="Is this library mate-paired?">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
136 <option value="single">Single-end</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
137 <option value="paired">Paired-end (as individual datasets)</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
138 <option value="paired_collection">Paired-end (as collection)</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
139 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
140 <when value="single">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
141 <param format="fastqsanger" name="input1" type="data" label="RNA-Seq FASTQ file" help="Nucleotide-space: Must have Sanger-scaled quality values with ASCII offset 33"/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
142 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
143 <when value="paired">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
144 <param format="fastqsanger" name="input1" type="data" label="RNA-Seq FASTQ file, forward reads" help="Nucleotide-space: Must have Sanger-scaled quality values with ASCII offset 33" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
145 <param format="fastqsanger" name="input2" type="data" label="RNA-Seq FASTQ file, reverse reads" help="Nucleotide-space: Must have Sanger-scaled quality values with ASCII offset 33" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
146 <expand macro="paired_parameters" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
147 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
148 <when value="paired_collection">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
149 <param format="fastqsanger" name="input" type="data_collection" collection_type="paired" label="RNA-Seq FASTQ paired reads" help="Nucleotide-space: Must have Sanger-scaled quality values with ASCII offset 33" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
150 <expand macro="paired_parameters" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
151 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
152 </conditional>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
153 <expand macro="refGenomeSourceConditional">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
154 <options from_data_table="tophat2_indexes">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
155 <filter type="sort_by" column="2"/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
156 <validator type="no_options" message="No genomes are available for the selected input dataset"/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
157 </options>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
158 </expand>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
159 <conditional name="params">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
160 <param name="settingsType" type="select" label="TopHat settings to use" help="You can use the default settings or set custom values for any of Tophat's parameters.">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
161 <option value="preSet">Use Defaults</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
162 <option value="full">Full parameter list</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
163 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
164 <when value="preSet" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
165 <!-- Full/advanced params. -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
166 <when value="full">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
167 <param name="read_realign_edit_dist" type="integer" value="1000" label="Max realign edit distance" help="Some of the reads spanning multiple exons may be mapped incorrectly as a contiguous alignment to the genome even though the correct alignment should be a spliced one - this can happen in the presence of processed pseudogenes that are rarely (if at all) transcribed or expressed. This option can direct TopHat to re-align reads for which the edit distance of an alignment obtained in a previous mapping step is above or equal to this option value. If you set this option to 0, TopHat will map every read in all the mapping steps (transcriptome if you provided gene annotations, genome, and finally splice variants detected by TopHat), reporting the best possible alignment found in any of these mapping steps. This may greatly increase the mapping accuracy at the expense of an increase in running time. The default value for this option is set such that TopHat will not try to realign reads already mapped in earlier steps." />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
168
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
169 <param name="read_edit_dist" type="integer" value="2" label="Max edit distance" help="Final read alignments having more than these many edit distance are discarded." />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
170
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
171 <param name="library_type" type="select" label="Library Type" help="TopHat will treat the reads as strand specific. Every read alignment will have an XS attribute tag. Consider supplying library type options below to select the correct RNA-seq protocol.">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
172 <option value="fr-unstranded">FR Unstranded</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
173 <option value="fr-firststrand">FR First Strand</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
174 <option value="fr-secondstrand">FR Second Strand</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
175 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
176 <param name="read_mismatches" type="integer" value="2" label="Final read mismatches" help="Final read alignments having more than these many mismatches are discarded." />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
177 <param name="bowtie_n" type="select" label="Use bowtie -n mode">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
178 <option selected="true" value="No">No</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
179 <option value="Yes">Yes</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
180 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
181 <param name="anchor_length" type="integer" value="8" label="Anchor length (at least 3)" help="Report junctions spanned by reads with at least this many bases on each side of the junction." />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
182 <param name="splice_mismatches" type="integer" value="0" label="Maximum number of mismatches that can appear in the anchor region of spliced alignment" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
183 <param name="min_intron_length" type="integer" value="70" label="The minimum intron length" help="TopHat will ignore donor/acceptor pairs closer than this many bases apart." />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
184 <param name="max_intron_length" type="integer" value="500000" label="The maximum intron length" help="When searching for junctions ab initio, TopHat will ignore donor/acceptor pairs farther than this many bases apart, except when such a pair is supported by a split segment alignment of a long read." />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
185 <expand macro="indel_searchConditional" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
186 alignments (number of reads divided by average depth of coverage)" help="0.0 to 1.0 (0 to turn off)" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
187 <param name="max_multihits" type="integer" value="20" label="Maximum number of alignments to be allowed" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
188 <param name="min_segment_intron" type="integer" value="50" label="Minimum intron length that may be found during split-segment (default) search" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
189 <param name="max_segment_intron" type="integer" value="500000" label="Maximum intron length that may be found during split-segment (default) search" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
190 <param name="seg_mismatches" type="integer" min="0" max="3" value="2" label="Number of mismatches allowed in each segment alignment for reads mapped independently" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
191 <param name="seg_length" type="integer" value="25" label="Minimum length of read segments" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
192
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
193 <!-- Options for supplying own junctions. -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
194 <expand macro="own_junctionsConditional" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
195 <!-- Coverage search. -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
196 <conditional name="coverage_search">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
197 <param name="use_search" type="select" label="Use Coverage Search" help="Enables the coverage based search for junctions. Use when coverage search is disabled by default (such as for reads 75bp or longer), for maximum sensitivity.">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
198 <option selected="true" value="No">No</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
199 <option value="Yes">Yes</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
200 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
201 <when value="Yes">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
202 <param name="min_coverage_intron" type="integer" value="50" label="Minimum intron length that may be found during coverage search" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
203 <param name="max_coverage_intron" type="integer" value="20000" label="Maximum intron length that may be found during coverage search" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
204 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
205 <when value="No" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
206 </conditional>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
207
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
208 <!-- Microexon search params -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
209 <param name="microexon_search" type="select" label="Use Microexon Search" help="With this option, the pipeline will attempt to find alignments incident to microexons. Works only for reads 50bp or longer.">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
210 <option value="No">No</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
211 <option value="Yes">Yes</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
212 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
213
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
214 <!-- Fusion mapping. -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
215 <conditional name="fusion_search">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
216 <param name="do_search" type="select" label="Do Fusion Search">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
217 <option selected="true" value="No">No</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
218 <option value="Yes">Yes</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
219 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
220 <when value="No" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
221 <when value="Yes">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
222 <param name="anchor_len" type="integer" value="20" label="Anchor Length" help="A 'supporting' read must map to both sides of a fusion by at least this many bases."/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
223 <param name="min_dist" type="integer" value="10000000" label="Minimum Distance" help="For intra-chromosomal fusions, TopHat-Fusion tries to find fusions separated by at least this distance."/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
224 <param name="read_mismatches" type="integer" value="2" label="Read Mismatches" help="Reads support fusions if they map across fusion with at most this many mismatches."/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
225 <param name="multireads" type="integer" value="2" label="Multireads" help="Reads that map to more than this many places will be ignored. It may be possible that a fusion is supported by reads (or pairs) that map to multiple places."/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
226 <param name="multipairs" type="integer" value="2" label="Multipairs" help="Pairs that map to more than this many places will be ignored."/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
227 <param name="ignore_chromosomes" type="text" value='' label="Ignore some chromosomes such as chrM when detecting fusion break points"/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
228 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
229 </conditional>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
230
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
231 <!-- Bowtie2 settings. -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
232 <conditional name="bowtie2_settings">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
233 <param name="b2_settings" type="select" label="Set Bowtie2 settings">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
234 <option selected="true" value="No">No</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
235 <option value="Yes">Yes</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
236 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
237 <when value="No" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
238 <when value="Yes">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
239 <conditional name="preset">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
240 <param name="b2_preset" type="select" label="Use Preset options">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
241 <option selected="true" value="Yes">Yes</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
242 <option value="No">No</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
243 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
244 <when value="Yes">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
245 <param name="b2_preset_select" type="select" label="Preset option">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
246 <option value="very-fast">Very fast</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
247 <option value="fast">Fast</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
248 <option selected="true" value="sensitive">Sensitive</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
249 <option value="very-sensitive">Very sensitive</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
250 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
251 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
252 <!-- TODO: -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
253 <when value="No" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
254 </conditional>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
255 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
256 </conditional>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
257 </when> <!-- full -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
258 </conditional> <!-- params -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
259 <conditional name="readGroup">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
260 <param name="specReadGroup" type="select" label="Specify read group?">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
261 <option value="yes">Yes</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
262 <option value="no" selected="True">No</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
263 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
264 <when value="yes">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
265 <param name="rgid" type="text" size="25" label="Read group identifier (ID). Each @RG line must have a unique ID. The value of ID is used in the RG tags of alignment records. Must be unique among all read groups in header section." help="Required if RG specified. Read group IDs may be modified when merging SAM files in order to handle collisions." />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
266 <param name="rglb" type="text" size="25" label="Library name (LB)" help="Required if RG specified" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
267 <param name="rgpl" type="text" size="25" label="Platform/technology used to produce the reads (PL)" help="Required if RG specified. Valid values : CAPILLARY, LS454, ILLUMINA, SOLID, HELICOS, IONTORRENT and PACBIO" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
268 <param name="rgsm" type="text" size="25" label="Sample (SM)" help="Required if RG specified. Use pool name where a pool is being sequenced" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
269 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
270 <when value="no" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
271 </conditional> <!-- readGroup -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
272 </inputs>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
273
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
274 <stdio>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
275 <regex match="Exception|Error" source="both" level="fatal" description="Tool execution failed"/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
276 <regex match=".*" source="both" level="log" description="tool progress"/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
277 </stdio>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
278
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
279 <outputs>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
280 <!--
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
281 <data format="txt" name="align_summary" label="${tool.name} on ${on_string}: align_summary" from_work_dir="tophat_out/align_summary.txt"/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
282 <data format="tabular" name="fusions" label="${tool.name} on ${on_string}: fusions" from_work_dir="tophat_out/fusions.out">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
283 <filter>(params['settingsType'] == 'full' and params['fusion_search']['do_search'] == 'Yes')</filter>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
284 </data>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
285 <data format="bed" name="insertions" label="${tool.name} on ${on_string}: insertions" from_work_dir="tophat_out/insertions.bed">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
286 <expand macro="dbKeyActions" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
287 </data>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
288 <data format="bed" name="deletions" label="${tool.name} on ${on_string}: deletions" from_work_dir="tophat_out/deletions.bed">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
289 <expand macro="dbKeyActions" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
290 </data>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
291 <data format="bed" name="junctions" label="${tool.name} on ${on_string}: splice junctions" from_work_dir="tophat_out/junctions.bed">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
292 <expand macro="dbKeyActions" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
293 </data>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
294 -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
295 <data format="bam" name="accepted_hits" label="${tool.name} on ${on_string}: accepted_hits" from_work_dir="tophat_out/accepted_hits.bam">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
296 <expand macro="dbKeyActions" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
297 </data>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
298 </outputs>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
299
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
300 <macros>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
301 <import>tophat_macros.xml</import>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
302 <xml name="paired_parameters">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
303 <param name="mate_inner_distance" type="integer" value="300" label="Mean Inner Distance between Mate Pairs" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
304 <param name="mate_std_dev" type="integer" value="20" label="Std. Dev for Distance between Mate Pairs" help="The standard deviation for the distribution on inner distances between mate pairs."/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
305 <!-- Discordant pairs. -->
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
306 <param name="report_discordant_pairs" type="select" label="Report discordant pair alignments?">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
307 <option value="No">No</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
308 <option selected="True" value="Yes">Yes</option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
309 </param>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
310 </xml>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
311 <macro name="dbKeyActions">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
312 <actions>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
313 <conditional name="refGenomeSource.genomeSource">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
314 <when value="indexed">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
315 <action type="metadata" name="dbkey">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
316 <option type="from_data_table" name="tophat2_indexes" column="1" offset="0">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
317 <filter type="param_value" column="0" value="#" compare="startswith" keep="False"/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
318 <filter type="param_value" ref="refGenomeSource.index" column="0"/>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
319 </option>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
320 </action>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
321 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
322 <when value="history">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
323 <action type="metadata" name="dbkey">
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
324 <option type="from_param" name="refGenomeSource.ownFile" param_attribute="dbkey" />
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
325 </action>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
326 </when>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
327 </conditional>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
328 </actions>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
329 </macro>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
330 </macros>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
331
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
332 <help>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
333 **Tophat Overview**
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
334
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
335 TopHat_ is a fast splice junction mapper for RNA-Seq reads. It aligns RNA-Seq reads to mammalian-sized genomes using the ultra high-throughput short read aligner Bowtie(2), and then analyzes the mapping results to identify splice junctions between exons. Please cite: Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, and Salzberg SL. TopHat2: accurate alignment
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
336 of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol 14:R36, 2013.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
337
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
338 .. _Tophat: http://ccb.jhu.edu/software/tophat/
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
339
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
340 ------
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
341
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
342 **Know what you are doing**
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
343
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
344 .. class:: warningmark
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
345
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
346 There is no such thing (yet) as an automated gearshift in splice junction identification. It is all like stick-shift driving in San Francisco. In other words, running this tool with default parameters will probably not give you meaningful results. A way to deal with this is to **understand** the parameters by carefully reading the `documentation`__ and experimenting. Fortunately, Galaxy makes experimenting easy.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
347
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
348 .. __: http://ccb.jhu.edu/software/tophat/manual.shtml
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
349
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
350 ------
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
351
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
352 **Input formats**
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
353
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
354 Tophat accepts files in Sanger FASTQ format. Use the FASTQ Groomer to prepare your files.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
355
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
356 ------
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
357
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
358 **Outputs**
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
359
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
360 Tophat produces two output files:
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
361
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
362 - junctions -- A UCSC BED_ track of junctions reported by TopHat. Each junction consists of two connected BED blocks, where each block is as long as the maximal overhang of any read spanning the junction. The score is the number of alignments spanning the junction.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
363 - accepted_hits -- A list of read alignments in BAM_ format.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
364
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
365 .. _BED: http://genome.ucsc.edu/FAQ/FAQformat.html#format1
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
366 .. _BAM: http://samtools.sourceforge.net/
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
367
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
368 Two other possible outputs, depending on the options you choose, are insertions and deletions, both of which are in BED format.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
369
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
370 -------
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
371
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
372 **Tophat settings**
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
373
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
374 All of the options have a default value. You can change any of them. Some of the options in Tophat have been implemented here.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
375
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
376 ------
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
377
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
378 **Tophat parameter list**
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
379
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
380 This is a list of implemented Tophat options::
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
381
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
382 -r This is the expected (mean) inner distance between mate pairs. For, example, for paired end runs with fragments
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
383 selected at 300bp, where each end is 50bp, you should set -r to be 200. There is no default, and this parameter
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
384 is required for paired end runs.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
385 --mate-std-dev INT The standard deviation for the distribution on inner distances between mate pairs. The default is 20bp.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
386 -a/--min-anchor-length INT The "anchor length". TopHat will report junctions spanned by reads with at least this many bases on each side of the junction. Note that individual spliced
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
387 alignments may span a junction with fewer than this many bases on one side. However, every junction involved in spliced alignments is supported by at least one
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
388 read with this many bases on each side. This must be at least 3 and the default is 8.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
389 -m/--splice-mismatches INT The maximum number of mismatches that may appear in the "anchor" region of a spliced alignment. The default is 0.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
390 -i/--min-intron-length INT The minimum intron length. TopHat will ignore donor/acceptor pairs closer than this many bases apart. The default is 70.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
391 -I/--max-intron-length INT The maximum intron length. When searching for junctions ab initio, TopHat will ignore donor/acceptor pairs farther than this many bases apart, except when such a pair is supported by a split segment alignment of a long read. The default is 500000.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
392 -g/--max-multihits INT Instructs TopHat to allow up to this many alignments to the reference for a given read, and suppresses all alignments for reads with more than this many
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
393 alignments. The default is 40.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
394 -G/--GTF [GTF 2.2 file] Supply TopHat with a list of gene model annotations. TopHat will use the exon records in this file to build a set of known splice junctions for each gene, and will attempt to align reads to these junctions even if they would not normally be covered by the initial mapping.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
395 -j/--raw-juncs [juncs file] Supply TopHat with a list of raw junctions. Junctions are specified one per line, in a tab-delimited format. Records look like: [chrom] [left] [right] [+/-], left and right are zero-based coordinates, and specify the last character of the left sequenced to be spliced to the first character of the right sequence, inclusive.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
396 -no-novel-juncs Only look for junctions indicated in the supplied GFF file. (ignored without -G)
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
397 --no-coverage-search Disables the coverage based search for junctions.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
398 --coverage-search Enables the coverage based search for junctions. Use when coverage search is disabled by default (such as for reads 75bp or longer), for maximum sensitivity.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
399 --microexon-search With this option, the pipeline will attempt to find alignments incident to microexons. Works only for reads 50bp or longer.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
400 --segment-mismatches Read segments are mapped independently, allowing up to this many mismatches in each segment alignment. The default is 2.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
401 --segment-length Each read is cut up into segments, each at least this long. These segments are mapped independently. The default is 25.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
402 --min-coverage-intron The minimum intron length that may be found during coverage search. The default is 50.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
403 --max-coverage-intron The maximum intron length that may be found during coverage search. The default is 20000.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
404 --min-segment-intron The minimum intron length that may be found during split-segment search. The default is 50.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
405 --max-segment-intron The maximum intron length that may be found during split-segment search. The default is 500000.
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
406 </help>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
407 <citations>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
408 <citation type="doi">10.1186/gb-2013-14-4-r36</citation>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
409 </citations>
589de4bdaca4 for Apache Aurora jobrunner
Makoto Nasuno <nasuno@ascade.co.jp>
parents:
diff changeset
410 </tool>