annotate bowtie_transcriptome_wrapper/bowtie_transcriptomic_wrapper.xml @ 16:bff7f8f0ba98 draft

Uploaded
author jackcurragh
date Thu, 09 Feb 2023 14:25:30 +0000
parents 9bef8486912e
children 816d9c5ef4a8
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
16
bff7f8f0ba98 Uploaded
jackcurragh
parents: 15
diff changeset
1 <tool id="bowtie_transcriptomic_wrapper" name="Bowtie Transcriptome Alignment" version="1.7.2">
9
15e6eb67d710 Uploaded
jackcurragh
parents: 8
diff changeset
2 <description>Align to Transcriptome using Bowtie</description>
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
3 <requirements>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
4 <requirement type="package" version="1.2.0">bowtie</requirement>
9
15e6eb67d710 Uploaded
jackcurragh
parents: 8
diff changeset
5 <requirement type="package" version="1.13">samtools</requirement>
15e6eb67d710 Uploaded
jackcurragh
parents: 8
diff changeset
6
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
7 </requirements>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
8 <version_command>bowtie --version</version_command>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
9 <command>
1
14fb4b4872e3 Uploaded
jackcurragh
parents: 0
diff changeset
10 python '$__tool_directory__/bowtie_transcriptome_wrapper.py'
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
11 ## Set number of threads
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
12 --threads="\${GALAXY_SLOTS:-4}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
13 ## Outputs
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
14 --output="${output}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
15 #if str( $singlePaired.sPaired ) == "single"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
16 #if $output_unmapped_reads_l
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
17 --output_unmapped_reads="${output_unmapped_reads_l}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
18 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
19 #if $output_suppressed_reads_l
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
20 --output_suppressed_reads="${output_suppressed_reads_l}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
21 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
22 --galaxy_input_format="${singlePaired.sInput1.ext}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
23 #else
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
24 #if $output_unmapped_reads_l and $output_unmapped_reads_r
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
25 --output_unmapped_reads_l="${output_unmapped_reads_l}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
26 --output_unmapped_reads_r="${output_unmapped_reads_r}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
27 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
28 #if $output_suppressed_reads_l and $output_suppressed_reads_l
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
29 --output_suppressed_reads_l="${output_suppressed_reads_l}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
30 --output_suppressed_reads_r="${output_suppressed_reads_r}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
31 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
32 --galaxy_input_format="${singlePaired.pInput1.ext}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
33 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
34 ## Inputs
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
35 --dataType="solexa" ##this indicates that nucleotide base space is used in the wrapper
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
36 --suppressHeader="${suppressHeader}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
37 --genomeSource="${refGenomeSource.genomeSource}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
38 #if $refGenomeSource.genomeSource == "history":
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
39 ##index already exists
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
40 #if $refGenomeSource.ownFile.extension.startswith( 'bowtie_' ):
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
41 ##user previously built
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
42 --ref="${refGenomeSource.ownFile.extra_files_path}/${refGenomeSource.ownFile.metadata.base_name}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
43 --do_not_build_index
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
44 #else:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
45 ##build index on the fly
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
46 --ref="${refGenomeSource.ownFile}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
47 --indexSettings="${refGenomeSource.indexParams.indexSettings}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
48 #if $refGenomeSource.indexParams.indexSettings == "indexFull":
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
49 --iautoB="${refGenomeSource.indexParams.autoBehavior.autoB}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
50 #if $refGenomeSource.indexParams.autoBehavior.autoB == "set":
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
51 --ipacked="${refGenomeSource.indexParams.autoBehavior.packed}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
52 --ibmax="${refGenomeSource.indexParams.autoBehavior.bmax}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
53 --ibmaxdivn="${refGenomeSource.indexParams.autoBehavior.bmaxdivn}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
54 --idcv="${refGenomeSource.indexParams.autoBehavior.dcv}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
55 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
56 --inodc="${refGenomeSource.indexParams.nodc}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
57 --inoref="${refGenomeSource.indexParams.noref}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
58 --ioffrate="${refGenomeSource.indexParams.offrate}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
59 --iftab="${refGenomeSource.indexParams.ftab}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
60 --intoa="${refGenomeSource.indexParams.ntoa}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
61 --iendian="${refGenomeSource.indexParams.endian}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
62 --iseed="${refGenomeSource.indexParams.seed}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
63 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
64 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
65 #else
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
66 ##use pre-built index
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
67 --ref="${refGenomeSource.index.fields.path}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
68 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
69 --paired="${singlePaired.sPaired}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
70 #if $singlePaired.sPaired == "single":
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
71 --input1="${singlePaired.sInput1}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
72 --params="${singlePaired.sParams.sSettingsType}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
73 #if $singlePaired.sParams.sSettingsType == "full":
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
74 --skip="${singlePaired.sParams.sSkip}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
75 --alignLimit="${singlePaired.sParams.sAlignLimit}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
76 --trimH="${singlePaired.sParams.sTrimH}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
77 --trimL="${singlePaired.sParams.sTrimL}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
78 #if $singlePaired.sParams.alignModeOption.alignMode == 'nMode'
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
79 --mismatchSeed="${singlePaired.sParams.alignModeOption.sMismatchSeed}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
80 --mismatchQual="${singlePaired.sParams.alignModeOption.sMismatchQual}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
81 --seedLen="${singlePaired.sParams.alignModeOption.sSeedLen}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
82 --rounding="${singlePaired.sParams.alignModeOption.sRounding}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
83 #else
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
84 --maxMismatches="${singlePaired.sParams.alignModeOption.maxMismatches}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
85 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
86 --forwardAlign="${singlePaired.sParams.sForwardAlign}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
87 --reverseAlign="${singlePaired.sParams.sReverseAlign}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
88 --tryHard="${singlePaired.sParams.sBestOption.sTryHardOption.sTryHard}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
89 --allValAligns="${singlePaired.sParams.sAllValAlignsOption.sAllValAligns}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
90 #if $singlePaired.sParams.sAllValAlignsOption.sAllValAligns == "noAllValAligns"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
91 --valAlign="${singlePaired.sParams.sAllValAlignsOption.sValAlign}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
92 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
93 --suppressAlign="${singlePaired.sParams.sSuppressAlign}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
94 --best="${singlePaired.sParams.sBestOption.sBest}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
95 #if $singlePaired.sParams.sBestOption.sBest == "doBest":
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
96 --strata="${singlePaired.sParams.sBestOption.sdStrata}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
97 #if $singlePaired.sParams.sBestOption.sTryHardOption.sTryHard == "noTryHard"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
98 --maxBacktracks="${singlePaired.sParams.sBestOption.sTryHardOption.sdMaxBacktracks}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
99 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
100 #else:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
101 #if $singlePaired.sParams.sBestOption.sTryHardOption.sTryHard == "noTryHard"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
102 --maxBacktracks="${singlePaired.sParams.sBestOption.sTryHardOption.snMaxBacktracks}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
103 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
104 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
105 --offrate="${singlePaired.sParams.sOffrate}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
106 --seed="${singlePaired.sParams.sSeed}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
107 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
108 #else:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
109 --input1="${singlePaired.pInput1}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
110 --input2="${singlePaired.pInput2}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
111 --maxInsert="${singlePaired.pMaxInsert}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
112 --mateOrient="${singlePaired.pMateOrient}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
113 --params="${singlePaired.pParams.pSettingsType}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
114 #if $singlePaired.pParams.pSettingsType == "full":
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
115 --skip="${singlePaired.pParams.pSkip}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
116 --alignLimit="${singlePaired.pParams.pAlignLimit}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
117 --trimH="${singlePaired.pParams.pTrimH}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
118 --trimL="${singlePaired.pParams.pTrimL}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
119 #if $singlePaired.pParams.alignModeOption.alignMode == 'nMode'
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
120 --mismatchSeed="${singlePaired.pParams.alignModeOption.pMismatchSeed}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
121 --mismatchQual="${singlePaired.pParams.alignModeOption.pMismatchQual}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
122 --seedLen="${singlePaired.pParams.alignModeOption.pSeedLen}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
123 --rounding="${singlePaired.pParams.alignModeOption.pRounding}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
124 #else
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
125 --maxMismatches="${singlePaired.pParams.alignModeOption.maxMismatches}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
126 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
127 --minInsert="${singlePaired.pParams.pMinInsert}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
128 --forwardAlign="${singlePaired.pParams.pForwardAlign}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
129 --reverseAlign="${singlePaired.pParams.pReverseAlign}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
130 --tryHard="${singlePaired.pParams.pBestOption.pTryHardOption.pTryHard}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
131 --allValAligns="${singlePaired.pParams.pAllValAlignsOption.pAllValAligns}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
132 #if $singlePaired.pParams.pAllValAlignsOption.pAllValAligns == "noAllValAligns"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
133 --valAlign="${singlePaired.pParams.pAllValAlignsOption.pValAlign}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
134 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
135 --suppressAlign="${singlePaired.pParams.pSuppressAlign}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
136 --best="${singlePaired.pParams.pBestOption.pBest}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
137 #if $singlePaired.pParams.pBestOption.pBest == "doBest":
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
138 --strata="${singlePaired.pParams.pBestOption.pdStrata}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
139 #if $singlePaired.pParams.pBestOption.pTryHardOption.pTryHard == "noTryHard"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
140 --maxAlignAttempt="${singlePaired.pParams.pBestOption.pTryHardOption.pMaxAlignAttempt}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
141 --maxBacktracks="${singlePaired.pParams.pBestOption.pTryHardOption.pdMaxBacktracks}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
142 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
143 #else:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
144 #if $singlePaired.pParams.pBestOption.pTryHardOption.pTryHard == "noTryHard"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
145 --maxAlignAttempt="${singlePaired.pParams.pBestOption.pTryHardOption.pMaxAlignAttempt}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
146 --maxBacktracks="${singlePaired.pParams.pBestOption.pTryHardOption.pnMaxBacktracks}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
147 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
148 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
149 --offrate="${singlePaired.pParams.pOffrate}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
150 --seed="${singlePaired.pParams.pSeed}"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
151 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
152 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
153 #if $save_mapping_stats
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
154 --output_mapping_stats="$mapping_stats"
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
155 #end if
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
156 </command>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
157 <inputs>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
158 <conditional name="refGenomeSource">
1
14fb4b4872e3 Uploaded
jackcurragh
parents: 0
diff changeset
159 <param name="genomeSource" type="select" label="Will you select a reference from your history or use a built-in index?" help="Built-ins were indexed using default options">
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
160 <option value="indexed">Use a built-in index</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
161 <option value="history">Use one from the history</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
162 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
163 <when value="indexed">
1
14fb4b4872e3 Uploaded
jackcurragh
parents: 0
diff changeset
164 <param name="index" type="select" label="Select a reference" help="if your reference of interest is not listed - contact RiboGalaxy team">
16
bff7f8f0ba98 Uploaded
jackcurragh
parents: 15
diff changeset
165 <options from_data_table="bowtie_transcriptome_indexes">
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
166 <filter type="sort_by" column="2" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
167 <validator type="no_options" message="No indexes are available" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
168 </options>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
169 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
170 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
171 <when value="history">
1
14fb4b4872e3 Uploaded
jackcurragh
parents: 0
diff changeset
172 <param name="ownFile" type="data" format="bowtie_base_index,fasta" label="Select the reference" />
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
173 <conditional name="indexParams">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
174 <param name="indexSettings" type="select" label="Choose whether to use Default options for building indices or to Set your own" help="These settings are ignored when using a prebuilt index">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
175 <option value="indexPreSet">Default</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
176 <option value="indexFull">Set your own</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
177 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
178 <when value="indexPreSet" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
179 <when value="indexFull">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
180 <conditional name="autoBehavior">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
181 <param name="autoB" type="select" label="Choose to use automatic or specified behavior for some parameters (-a)" help="Allows you to set --packed, --bmax, --bmaxdivn, and --dcv">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
182 <option value="auto">Automatic behavior</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
183 <option value="set">Set values (sets --noauto and allows others to be set)</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
184 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
185 <when value="auto" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
186 <when value="set">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
187 <param name="packed" type="select" label="Whether or not to use a packed representation for DNA strings (--packed)" help="Packed representation saves memory but makes indexing 2-3 times slower">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
188 <option value="unpacked">Use regular representation</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
189 <option value="packed">Use packed representation</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
190 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
191 <param name="bmax" type="integer" value="-1" label="Maximum number of suffixes allowed in a block (--bmax)" help="-1 for not specified. Must be at least 1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
192 <param name="bmaxdivn" type="integer" value="4" label="Maximum number of suffixes allowed in a block as a fraction of the length of the reference (--bmaxdivn)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
193 <param name="dcv" type="integer" value="1024" min="3" label="The period for the difference-cover sample (--dcv)" help="A larger period yields less memory overhead, but may make suffix sorting slower, especially if repeats are present" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
194 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
195 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
196 <param name="nodc" type="select" label="Whether or not to disable the use of the difference-cover sample (--nodc)" help="Suffix sorting becomes quadratic-time in the worst case (with a very repetitive reference)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
197 <option value="dc">Use difference-cover sample</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
198 <option value="nodc">Disable difference-cover sample</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
199 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
200 <param name="noref" type="select" label="Whether or not to build the part of the reference index used only in paired-end alignment (-r)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
201 <option value="ref">Build all index files</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
202 <option value="noref">Do not build paired-end alignment index files</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
203 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
204 <param name="offrate" type="integer" value="5" min="0" label="The indexer will mark every 2^n Burrows-Wheeler rows with their corresponding location on the genome (-o)" help="Marking more rows makes reference-position lookups faster, but requires more memory to hold the annotations at runtime" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
205 <param name="ftab" type="integer" value="10" min="1" label="The size of the ftab lookup table used to calculate an initial Burrows-Wheeler range with respect to the first n characters of the query (-t)" help="ftab size is 4^(n+1) bytes" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
206 <param name="ntoa" type="select" label="Whether or not to convert Ns in the reference sequence to As (--ntoa)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
207 <option value="no">Do not convert Ns</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
208 <option value="yes">Convert Ns to As</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
209 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
210 <param name="endian" type="select" label="Endianness to use when serializing integers to the index file (--big/--little)" help="Little is most appropriate for Intel- and AMD-based architecture">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
211 <option value="little">Little</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
212 <option value="big">Big</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
213 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
214 <param name="seed" type="integer" value="-1" label="Seed for the pseudorandom number generator (--seed)" help="Use -1 to use default" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
215 </when> <!-- indexFull -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
216 </conditional> <!-- indexParams -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
217 </when> <!-- history -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
218 </conditional> <!-- refGenomeSource -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
219 <conditional name="singlePaired">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
220 <param name="sPaired" type="select" label="Is this library mate-paired?">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
221 <option value="single">Single-end</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
222 <option value="paired">Paired-end</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
223 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
224 <when value="single">
10
5f3a3790cc79 Uploaded
jackcurragh
parents: 9
diff changeset
225 <param name="sInput1" type="data" format="fastq,fastqsanger,fastqillumina,fastqsolexa" label="FASTQ file" help="Must have ASCII encoded quality scores"/>
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
226 <conditional name="sParams">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
227 <param name="sSettingsType" type="select" label="Bowtie settings to use" help="For most mapping needs use Commonly used settings. If you want full control use Full parameter list">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
228 <option value="preSet">Commonly used</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
229 <option value="full" selected="true">Full parameter list</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
230 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
231 <when value="preSet" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
232 <when value="full">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
233 <param name="sSkip" type="integer" value="0" label="Skip the first n reads (-s)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
234 <param name="sAlignLimit" type="integer" value="-1" label="Only align the first n reads (-u)" help="-1 for off" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
235 <param name="sTrimH" type="integer" value="0" label="Trim n bases from high-quality (left) end of each read before alignment (-5)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
236 <param name="sTrimL" type="integer" value="0" label="Trim n bases from low-quality (right) end of each read before alignment (-3)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
237 <conditional name="alignModeOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
238 <param name="alignMode" type="select" label="Alignment mode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
239 <option value="nMode">Maq-like: quality-aware, limit mismatches in seed (-n)</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
240 <option value="vMode" selected="true">ignore qualities, limit end-to-end mismatches (-v)</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
241 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
242 <when value="nMode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
243 <param name="sMismatchSeed" type="integer" value="2" min="0" max="3" label="Maximum number of mismatches permitted in the seed (-n)" help="May be 0, 1, 2, or 3" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
244 <param name="sMismatchQual" type="integer" value="70" min="1" label="Maximum permitted total of quality values at all mismatched read positions (-e)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
245 <param name="sSeedLen" type="integer" value="28" min="5" label="Seed length (-l)" help="Minimum value is 5" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
246 <param name="sRounding" type="select" label="Whether or not to round to the nearest 10 and saturating at 30 (--nomaqround)" help="Maq accepts quality values in the Phred quality scale, but internally rounds values to the nearest 10, with a maximum of 30. By default, bowtie also rounds this way">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
247 <option value="round">Round to nearest 10</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
248 <option value="noRound">Do not round to nearest 10</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
249 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
250 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
251 <when value="vMode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
252 <param name="maxMismatches" type="integer" value="3" min="0" max="3" label="Maximum number of mismatches (-v)" help="May be 0, 1, 2, or 3" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
253 <param name="sSeedLen" type="integer" value="25" min="5" label="Seed length (-l)" help="Minimum value is 5" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
254 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
255 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
256 <param name="sForwardAlign" type="select" label="Choose whether or not to attempt to align against the forward reference strand (--nofw)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
257 <option value="forward">Align against the forward reference strand</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
258 <option value="noForward">Do not align against the forward reference strand</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
259 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
260 <param name="sReverseAlign" type="select" label="Choose whether or not to attempt to align against the reverse-complement reference strand (--norc)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
261 <option value="reverse">Align against the reverse-complement reference strand</option>
1
14fb4b4872e3 Uploaded
jackcurragh
parents: 0
diff changeset
262 <option value="noReverse" selected="true">Do not align against the reverse-complement reference strand</option>
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
263 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
264 <conditional name="sBestOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
265 <param name="sBest" type="select" label="Whether or not to make Bowtie guarantee that reported singleton alignments are 'best' in terms of stratum and in terms of the quality values at the mismatched positions (--best)" help="Removes all strand bias. Only affects which alignments are reported by Bowtie. Runs slower with best option">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
266 <option value="noBest">Do not use best</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
267 <option value="doBest">Use best</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
268 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
269 <when value="noBest">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
270 <conditional name="sTryHardOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
271 <param name="sTryHard" type="select" label="Whether or not to try as hard as possible to find valid alignments when they exist (-y)" help="Tryhard mode is much slower than regular mode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
272 <option value="noTryHard">Do not try hard</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
273 <option value="doTryHard">Try hard</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
274 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
275 <when value="noTryHard">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
276 <param name="snMaxBacktracks" type="integer" value="125" min="0" label="Maximum number of backtracks permitted when aligning a read (--maxbts)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
277 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
278 <when value="doTryHard" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
279 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
280 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
281 <when value="doBest">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
282 <param name="sdStrata" type="select" label="Whether or not to report only those alignments that fall in the best stratum if many valid alignments exist and are reportable (--strata)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
283 <option value="noStrata">Do not use strata option</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
284 <option value="doStrata">Use strata option</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
285 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
286 <conditional name="sTryHardOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
287 <param name="sTryHard" type="select" label="Whether or not to try as hard as possible to find valid alignments when they exist (-y)" help="Tryhard mode is much slower than regular mode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
288 <option value="noTryHard">Do not try hard</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
289 <option value="doTryHard">Try hard</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
290 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
291 <when value="noTryHard">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
292 <param name="sdMaxBacktracks" type="integer" value="800" min="0" label="Maximum number of backtracks permitted when aligning a read (--maxbts)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
293 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
294 <when value="doTryHard" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
295 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
296 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
297 </conditional> <!-- bestOption -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
298 <conditional name="sAllValAlignsOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
299 <param name="sAllValAligns" type="select" label="Whether or not to report all valid alignments per read (-a)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
300 <option value="noAllValAligns">Do not report all valid alignments</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
301 <option value="doAllValAligns" selected="true">Report all valid alignments</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
302 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
303 <when value="noAllValAligns">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
304 <param name="sValAlign" type="integer" value="1" min="1" label="Report up to n valid alignments per read (-k)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
305 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
306 <when value="doAllValAligns" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
307 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
308 <param name="sSuppressAlign" type="integer" value="-1" label="Suppress all alignments for a read if more than n reportable alignments exist (-m)" help="-1 for no limit" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
309 <param name="sMaxFile" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Write all reads with a number of valid alignments exceeding the limit set with the -m option to a file (--max)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
310 <param name="sUnmappedFile" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Write all reads that could not be aligned to a file (--un)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
311 <param name="sOffrate" type="integer" value="-1" label="Override the offrate of the index to n (-o)" help="-1 for default" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
312 <param name="sSeed" type="integer" value="-1" label="Seed for pseudo-random number generator (--seed)" help="-1 for default" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
313 </when> <!-- full -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
314 </conditional> <!-- sParams -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
315 </when> <!-- single -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
316 <when value="paired">
10
5f3a3790cc79 Uploaded
jackcurragh
parents: 9
diff changeset
317 <param name="pInput1" type="data" format="fastq,fastqsanger,fastqillumina,fastqsolexa" label="Forward FASTQ file" help="Must have ASCII encoded quality scores"/>
5f3a3790cc79 Uploaded
jackcurragh
parents: 9
diff changeset
318 <param name="pInput2" type="data" format="fastq,fastqsanger,fastqillumina,fastqsolexa" label="Reverse FASTQ file" help="File format must match the Forward FASTQ file">
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
319 <options options_filter_attribute="ext" from_parameter="tool.app.datatypes_registry.datatypes_by_extension" transform_lines="obj.keys()">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
320 <column name="name" index="0"/>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
321 <column name="value" index="0"/>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
322 <filter type="param_value" ref="pInput1" ref_attribute="ext" column="0"/>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
323 </options>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
324 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
325 <param name="pMaxInsert" type="integer" value="1000" label="Maximum insert size for valid paired-end alignments (-X)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
326 <param name="pMateOrient" type="select" label="The upstream/downstream mate orientation for valid paired-end alignment against the forward reference strand (--fr/--rf/--ff)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
327 <option value="fr">FR (for Illumina)</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
328 <option value="rf">RF</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
329 <option value="ff">FF (for SOLiD)</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
330 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
331 <conditional name="pParams">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
332 <param name="pSettingsType" type="select" label="Bowtie settings to use" help="For most mapping needs use Commonly used settings. If you want full control use Full parameter list">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
333 <option value="preSet">Commonly used</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
334 <option value="full">Full parameter list</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
335 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
336 <when value="preSet" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
337 <when value="full">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
338 <param name="pSkip" type="integer" value="0" label="Skip the first n pairs (-s)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
339 <param name="pAlignLimit" type="integer" value="-1" label="Only align the first n pairs (-u)" help="-1 for off" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
340 <param name="pTrimH" type="integer" value="0" label="Trim n bases from high-quality (left) end of each read before alignment (-5)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
341 <param name="pTrimL" type="integer" value="0" label="Trim n bases from low-quality (right) end of each read before alignment (-3)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
342 <conditional name="alignModeOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
343 <param name="alignMode" type="select" label="Alignment mode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
344 <option value="nMode">Maq-like: quality-aware, limit mismatches in seed (-n)</option>
1
14fb4b4872e3 Uploaded
jackcurragh
parents: 0
diff changeset
345 <option value="vMode" selected="true">ignore qualities, limit end-to-end mismatches (-v)</option>
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
346 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
347 <when value="nMode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
348 <param name="pMismatchSeed" type="integer" value="2" min="0" max="3" label="Maximum number of mismatches permitted in the seed (-n)" help="May be 0, 1, 2, or 3" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
349 <param name="pMismatchQual" type="integer" value="70" min="1" label="Maximum permitted total of quality values at all mismatched read positions (-e)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
350 <param name="pSeedLen" type="integer" value="28" min="5" label="Seed length (-l)" help="Minimum value is 5" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
351 <param name="pRounding" type="select" label="Whether or not to round to the nearest 10 and saturating at 30 (--nomaqround)" help="Maq accepts quality values in the Phred quality scale, but internally rounds values to the nearest 10, with a maximum of 30. By default, bowtie also rounds this way">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
352 <option value="round">Round to nearest 10</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
353 <option value="noRound">Do not round to nearest 10</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
354 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
355 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
356 <when value="vMode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
357 <param name="maxMismatches" type="integer" value="" min="0" max="3" label="Maximum number of mismatches (-v)" help="May be 0, 1, 2, or 3" />
1
14fb4b4872e3 Uploaded
jackcurragh
parents: 0
diff changeset
358 <param name="pSeedLen" type="integer" value="25" min="5" label="Seed length (-l)" help="Minimum value is 5" />
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
359 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
360 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
361 <param name="pMinInsert" type="integer" value="0" label="Minimum insert size for valid paired-end alignments (-I)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
362 <param name="pForwardAlign" type="select" label="Choose whether or not to attempt to align against the forward reference strand (--nofw)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
363 <option value="forward">Align against the forward reference strand</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
364 <option value="noForward">Do not align against the forward reference strand</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
365 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
366 <param name="pReverseAlign" type="select" label="Choose whether or not to attempt to align against the reverse-complement reference strand (--norc)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
367 <option value="reverse">Align against the reverse-complement reference strand</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
368 <option value="noReverse" selected="true">Do not align against the reverse-complement reference strand</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
369 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
370 <conditional name="pBestOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
371 <param name="pBest" type="select" label="Whether or not to make Bowtie guarantee that reported singleton alignments are 'best' in terms of stratum and in terms of the quality values at the mismatched positions (--best)" help="Removes all strand bias. Only affects which alignments are reported by Bowtie. Runs slower with best option">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
372 <option value="noBest">Do not use best</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
373 <option value="doBest">Use best</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
374 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
375 <when value="noBest">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
376 <conditional name="pTryHardOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
377 <param name="pTryHard" type="select" label="Whether or not to try as hard as possible to find valid alignments when they exist (-y)" help="Tryhard mode is much slower than regular mode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
378 <option value="noTryHard">Do not try hard</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
379 <option value="doTryHard">Try hard</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
380 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
381 <when value="noTryHard">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
382 <param name="pMaxAlignAttempt" type="integer" value="100" min="1" label="Maximum number of attempts Bowtie will make to match an alignment for one mate with an alignment for the opposite mate (--pairtries)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
383 <param name="pnMaxBacktracks" type="integer" value="125" min="0" label="Maximum number of backtracks permitted when aligning a read (--maxbts)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
384 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
385 <when value="doTryHard" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
386 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
387 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
388 <when value="doBest">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
389 <param name="pdStrata" type="select" label="Whether or not to report only those alignments that fall in the best stratum if many valid alignments exist and are reportable (--strata)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
390 <option value="noStrata">Do not use strata option</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
391 <option value="doStrata">Use strata option</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
392 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
393 <conditional name="pTryHardOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
394 <param name="pTryHard" type="select" label="Whether or not to try as hard as possible to find valid alignments when they exist (-y)" help="Tryhard mode is much slower than regular mode">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
395 <option value="noTryHard">Do not try hard</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
396 <option value="doTryHard">Try hard</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
397 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
398 <when value="noTryHard">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
399 <param name="pMaxAlignAttempt" type="integer" value="100" min="1" label="Maximum number of attempts Bowtie will make to match an alignment for one mate with an alignment for the opposite mate (--pairtries)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
400 <param name="pdMaxBacktracks" type="integer" value="800" min="0" label="Maximum number of backtracks permitted when aligning a read (--maxbts)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
401 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
402 <when value="doTryHard" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
403 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
404 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
405 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
406 <conditional name="pAllValAlignsOption">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
407 <param name="pAllValAligns" type="select" label="Whether or not to report all valid alignments per pair (-a)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
408 <option value="noAllValAligns">Do not report all valid alignments</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
409 <option value="doAllValAligns">Report all valid alignments</option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
410 </param>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
411 <when value="noAllValAligns">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
412 <param name="pValAlign" type="integer" value="1" min="1" label="Report up to n valid alignments per pair (-k)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
413 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
414 <when value="doAllValAligns" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
415 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
416 <param name="pSuppressAlign" type="integer" value="-1" label="Suppress all alignments for a pair if more than n reportable alignments exist (-m)" help="-1 for no limit" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
417 <param name="pMaxFile" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Write all reads with a number of valid alignments exceeding the limit set with the -m option to a file (--max)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
418 <param name="pUnmappedFile" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Write all reads that could not be aligned to a file (--un)" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
419 <param name="pOffrate" type="integer" value="-1" label="Override the offrate of the index to n (-o)" help="-1 for default" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
420 <param name="pSeed" type="integer" value="-1" label="Seed for pseudo-random number generator (--seed)" help="-1 for default" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
421 </when> <!-- full -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
422 </conditional> <!-- pParams -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
423 </when> <!-- paired -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
424 </conditional> <!-- singlePaired -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
425 <param name="save_mapping_stats" type="boolean" checked="False" label="Save the bowtie mapping statistics to the history" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
426 <param name="suppressHeader" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Suppress the header in the output SAM file (--sam-nohead)" help="Bowtie produces SAM with several lines of header information by default" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
427 </inputs>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
428 <outputs>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
429 <data format="sam" name="output" label="${tool.name} on ${on_string}: mapped reads">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
430 <actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
431 <conditional name="refGenomeSource.genomeSource">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
432 <when value="indexed">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
433 <action type="metadata" name="dbkey">
7
d57209761277 Uploaded
jackcurragh
parents: 3
diff changeset
434 <option type="from_data_table" name="bowtie_transcriptome_indexes" column="1" offset="0">
0
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
435 <filter type="param_value" column="0" value="#" compare="startswith" keep="False"/>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
436 <filter type="param_value" ref="refGenomeSource.index" column="0"/>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
437 </option>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
438 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
439 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
440 <when value="history">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
441 <action type="metadata" name="dbkey">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
442 <option type="from_param" name="refGenomeSource.ownFile" param_attribute="dbkey" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
443 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
444 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
445 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
446 </actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
447 </data>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
448 <data format="txt" name="mapping_stats" label="${tool.name} on ${on_string}: mapping stats">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
449 <filter>save_mapping_stats is True</filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
450 </data>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
451 <data format="fastq" name="output_suppressed_reads_l" label="${tool.name} on ${on_string}: suppressed reads (L)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
452 <filter>((
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
453 singlePaired['sPaired'] == "single" and
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
454 singlePaired['sParams']['sSettingsType'] == "full" and
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
455 singlePaired['sParams']['sMaxFile'] is True
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
456 ) or (
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
457 singlePaired['sPaired'] == "paired" and
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
458 singlePaired['pParams']['pSettingsType'] == "full" and
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
459 singlePaired['pParams']['pMaxFile'] is True
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
460 ))
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
461 </filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
462 <actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
463 <conditional name="singlePaired.sPaired">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
464 <when value="single">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
465 <action type="format">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
466 <option type="from_param" name="singlePaired.sInput1" param_attribute="ext" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
467 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
468 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
469 <when value="paired">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
470 <action type="format">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
471 <option type="from_param" name="singlePaired.pInput1" param_attribute="ext" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
472 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
473 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
474 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
475 </actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
476 </data>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
477 <data format="fastq" name="output_suppressed_reads_r" label="${tool.name} on ${on_string}: suppressed reads (R)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
478 <filter>singlePaired['sPaired'] == "paired"</filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
479 <filter>singlePaired['pParams']['pSettingsType'] == "full"</filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
480 <filter>singlePaired['pParams']['pMaxFile'] is True</filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
481 <actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
482 <conditional name="singlePaired.sPaired">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
483 <when value="single">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
484 <action type="format">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
485 <option type="from_param" name="singlePaired.sInput1" param_attribute="ext" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
486 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
487 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
488 <when value="paired">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
489 <action type="format">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
490 <option type="from_param" name="singlePaired.pInput1" param_attribute="ext" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
491 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
492 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
493 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
494 </actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
495 </data>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
496 <data format="fastq" name="output_unmapped_reads_l" label="${tool.name} on ${on_string}: unmapped reads (L)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
497 <filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
498 ((
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
499 singlePaired['sPaired'] == "single" and
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
500 singlePaired['sParams']['sSettingsType'] == "full" and
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
501 singlePaired['sParams']['sUnmappedFile'] is True
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
502 ) or (
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
503 singlePaired['sPaired'] == "paired" and
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
504 singlePaired['pParams']['pSettingsType'] == "full" and
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
505 singlePaired['pParams']['pUnmappedFile'] is True
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
506 ))
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
507 </filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
508 <actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
509 <conditional name="singlePaired.sPaired">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
510 <when value="single">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
511 <action type="format">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
512 <option type="from_param" name="singlePaired.sInput1" param_attribute="ext" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
513 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
514 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
515 <when value="paired">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
516 <action type="format">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
517 <option type="from_param" name="singlePaired.pInput1" param_attribute="ext" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
518 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
519 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
520 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
521 </actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
522 </data>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
523 <data format="fastq" name="output_unmapped_reads_r" label="${tool.name} on ${on_string}: unmapped reads (R)">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
524 <filter>singlePaired['sPaired'] == "paired"</filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
525 <filter>singlePaired['pParams']['pSettingsType'] == "full"</filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
526 <filter>singlePaired['pParams']['pUnmappedFile'] is True</filter>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
527 <actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
528 <conditional name="singlePaired.sPaired">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
529 <when value="single">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
530 <action type="format">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
531 <option type="from_param" name="singlePaired.sInput1" param_attribute="ext" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
532 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
533 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
534 <when value="paired">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
535 <action type="format">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
536 <option type="from_param" name="singlePaired.pInput1" param_attribute="ext" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
537 </action>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
538 </when>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
539 </conditional>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
540 </actions>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
541 </data>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
542 </outputs>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
543 <tests>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
544 <test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
545 <!--
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
546 Bowtie command:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
547 bowtie -q -p 4 -S +sam-nohead chrM_base test-data/bowtie_in2.fastqsanger > bowtie_out6_u.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
548 sort bowtie_out6_u.sam > bowtie_out6.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
549 -p is the number of threads. You need to replace the + with 2 dashes.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
550 chrM_base needs to be the base location/name of the index files.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
551 -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
552 <param name="genomeSource" value="indexed" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
553 <!-- this is the backwards-compatible "unique value" for this index, not an actual path -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
554 <param name="index" value="equCab2chrM" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
555 <param name="sPaired" value="single" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
556 <param name="sInput1" ftype="fastqsanger" value="bowtie_in2.fastqsanger" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
557 <param name="sSettingsType" value="preSet" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
558 <param name="suppressHeader" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
559 <output name="output" ftype="sam" file="bowtie_out6.sam" sort="True">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
560 <metadata name="dbkey" value="equCab2" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
561 </output>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
562 </test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
563 <test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
564 <!--
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
565 Bowtie command:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
566 bowtie-build -f test-data/phiX.fasta phiX_base
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
567 bowtie -q -X 1000 +ff -p 4 -S +sam-nohead -n 2 -e 70 -l 28 +pairtries 100 +maxbts 800 +best +un bowtie_out8_u.fastq phiX_base -1 test-data/bowtie_in5.fastqsanger -2 test-data/bowtie_in6.fastqsanger > bowtie_out7_u.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
568 sort bowtie_out7_u.sam > bowtie_out7.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
569 sort bowtie_out8_u_1.sam > bowtie_out8_1.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
570 sort bowtie_out8_u_2.sam > bowtie_out8_2.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
571 Then also need to modify bowtie_out8_1.sam and bowtie_out8_2.sam so that all @ lines come before sequence lines.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
572 -p is the number of threads. You need to replace the + with 2 dashes.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
573 The two unmapped output files will be named bowtie_out8_1.fastq and bowtie_out8_2.fastq.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
574 chrM_base is the index files' location/base name.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
575 -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
576 <param name="genomeSource" value="history" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
577 <param name="ownFile" value="phiX.fasta" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
578 <param name="indexSettings" value="indexPreSet" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
579 <param name="sPaired" value="paired" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
580 <param name="pInput1" ftype="fastqsanger" value="bowtie_in5.fastqsanger" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
581 <param name="pInput2" ftype="fastqsanger" value="bowtie_in6.fastqsanger" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
582 <param name="pMaxInsert" value="1000" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
583 <param name="pMateOrient" value="ff" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
584 <param name="pSettingsType" value="full" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
585 <param name="pSkip" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
586 <param name="pAlignLimit" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
587 <param name="pTrimH" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
588 <param name="pTrimL" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
589 <param name="alignMode" value="nMode" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
590 <param name="pMismatchSeed" value="2" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
591 <param name="pMismatchQual" value="70" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
592 <param name="pSeedLen" value="28" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
593 <param name="pRounding" value="round" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
594 <param name="pMinInsert" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
595 <param name="pMaxAlignAttempt" value="100" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
596 <param name="pForwardAlign" value="forward" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
597 <param name="pReverseAlign" value="reverse" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
598 <param name="pTryHard" value="noTryHard" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
599 <param name="pValAlign" value="1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
600 <param name="pAllValAligns" value="noAllValAligns" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
601 <param name="pSuppressAlign" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
602 <param name="pUnmappedFile" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
603 <param name="pMaxFile" value="false" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
604 <param name="pBest" value="doBest" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
605 <param name="pdMaxBacktracks" value="800" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
606 <param name="pdStrata" value="noStrata" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
607 <param name="pOffrate" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
608 <param name="pSeed" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
609 <param name="suppressHeader" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
610 <output name="output" ftype="sam" file="bowtie_out7.sam" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
611 <output name="output_unmapped_reads_l" ftype="fastqsanger" file="bowtie_out8_1.fastq" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
612 <output name="output_unmapped_reads_r" ftype="fastqsanger" file="bowtie_out8_2.fastq" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
613 </test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
614 <!-- start testing of non-sanger variant fastq reads -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
615 <test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
616 <param name="genomeSource" value="history" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
617 <param name="ownFile" value="phiX.fasta" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
618 <param name="indexSettings" value="indexPreSet" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
619 <param name="sPaired" value="paired" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
620 <param name="pInput1" ftype="fastqillumina" value="bowtie_in5.fastqillumina" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
621 <param name="pInput2" ftype="fastqillumina" value="bowtie_in6.fastqillumina" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
622 <param name="pMaxInsert" value="1000" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
623 <param name="pMateOrient" value="ff" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
624 <param name="pSettingsType" value="full" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
625 <param name="pSkip" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
626 <param name="pAlignLimit" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
627 <param name="pTrimH" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
628 <param name="pTrimL" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
629 <param name="alignMode" value="nMode" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
630 <param name="pMismatchSeed" value="2" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
631 <param name="pMismatchQual" value="70" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
632 <param name="pSeedLen" value="28" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
633 <param name="pRounding" value="round" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
634 <param name="pMinInsert" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
635 <param name="pMaxAlignAttempt" value="100" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
636 <param name="pForwardAlign" value="forward" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
637 <param name="pReverseAlign" value="reverse" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
638 <param name="pTryHard" value="noTryHard" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
639 <param name="pValAlign" value="1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
640 <param name="pAllValAligns" value="noAllValAligns" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
641 <param name="pSuppressAlign" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
642 <param name="pUnmappedFile" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
643 <param name="pMaxFile" value="false" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
644 <param name="pBest" value="doBest" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
645 <param name="pdMaxBacktracks" value="800" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
646 <param name="pdStrata" value="noStrata" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
647 <param name="pOffrate" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
648 <param name="pSeed" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
649 <param name="suppressHeader" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
650 <output name="output" ftype="sam" file="bowtie_out7.sam" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
651 <output name="output_unmapped_reads_l" ftype="fastqillumina" file="bowtie_out8_1.fastqillumina.sorted" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
652 <output name="output_unmapped_reads_r" ftype="fastqillumina" file="bowtie_out8_2.fastqillumina.sorted" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
653 </test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
654 <test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
655 <param name="genomeSource" value="history" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
656 <param name="ownFile" value="phiX.fasta" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
657 <param name="indexSettings" value="indexPreSet" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
658 <param name="sPaired" value="paired" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
659 <param name="pInput1" ftype="fastqsolexa" value="bowtie_in5.fastqsolexa" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
660 <param name="pInput2" ftype="fastqsolexa" value="bowtie_in6.fastqsolexa" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
661 <param name="pMaxInsert" value="1000" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
662 <param name="pMateOrient" value="ff" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
663 <param name="pSettingsType" value="full" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
664 <param name="pSkip" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
665 <param name="pAlignLimit" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
666 <param name="pTrimH" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
667 <param name="pTrimL" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
668 <param name="alignMode" value="nMode" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
669 <param name="pMismatchSeed" value="2" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
670 <param name="pMismatchQual" value="70" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
671 <param name="pSeedLen" value="28" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
672 <param name="pRounding" value="round" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
673 <param name="pMinInsert" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
674 <param name="pMaxAlignAttempt" value="100" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
675 <param name="pForwardAlign" value="forward" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
676 <param name="pReverseAlign" value="reverse" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
677 <param name="pTryHard" value="noTryHard" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
678 <param name="pValAlign" value="1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
679 <param name="pAllValAligns" value="noAllValAligns" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
680 <param name="pSuppressAlign" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
681 <param name="pUnmappedFile" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
682 <param name="pMaxFile" value="false" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
683 <param name="pBest" value="doBest" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
684 <param name="pdMaxBacktracks" value="800" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
685 <param name="pdStrata" value="noStrata" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
686 <param name="pOffrate" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
687 <param name="pSeed" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
688 <param name="suppressHeader" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
689 <output name="output" ftype="sam" file="bowtie_out7.sam" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
690 <output name="output_unmapped_reads_l" ftype="fastqsolexa" file="bowtie_out8_1.fastqsolexa.sorted" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
691 <output name="output_unmapped_reads_r" ftype="fastqsolexa" file="bowtie_out8_2.fastqsolexa.sorted" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
692 </test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
693 <!-- end testing of non-sanger variant fastq reads -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
694 <test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
695 <!--
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
696 Bowtie command:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
697 bowtie -q -p 4 -S +sam-nohead -n 2 -e 70 -l 28 -y -k 1 chrM_base test-data/bowtie_in2.fastqsanger > bowtie_out9_u.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
698 sort bowtie_out9_u.sam > bowtie_out9.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
699 -p is the number of threads. You need to replace the + with 2 dashes.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
700 chrM_base is the index files' location/base name.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
701 -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
702 <param name="genomeSource" value="indexed" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
703 <!-- this is the backwards-compatible "unique value" for this index, not an actual path -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
704 <param name="index" value="equCab2chrM" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
705 <param name="sPaired" value="single" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
706 <param name="sInput1" ftype="fastqsanger" value="bowtie_in2.fastqsanger" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
707 <param name="sSettingsType" value="full" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
708 <param name="sSkip" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
709 <param name="sAlignLimit" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
710 <param name="sTrimH" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
711 <param name="sTrimL" value="0" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
712 <param name="alignMode" value="nMode" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
713 <param name="sMismatchSeed" value="2" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
714 <param name="sMismatchQual" value="70" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
715 <param name="sSeedLen" value="28" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
716 <param name="sRounding" value="round" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
717 <param name="sForwardAlign" value="forward" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
718 <param name="sReverseAlign" value="reverse" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
719 <param name="sTryHard" value="doTryHard" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
720 <param name="sValAlign" value="1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
721 <param name="sAllValAligns" value="noAllValAligns" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
722 <param name="sSuppressAlign" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
723 <param name="sUnmappedFile" value="false" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
724 <param name="sMaxFile" value="false" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
725 <param name="sBest" value="noBest" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
726 <param name="sOffrate" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
727 <param name="sSeed" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
728 <param name="suppressHeader" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
729 <output name="output" ftype="sam" file="bowtie_out9.sam" sort="True">
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
730 <metadata name="dbkey" value="equCab2" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
731 </output>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
732 </test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
733 <test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
734 <!--
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
735 Bowtie command:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
736 bowtie-build +offrate 5 +ftabchars 10 +little -f test-data/phiX.fasta phiX_base
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
737 bowtie -q -X 1000 +ff -p 4 -S +sam-nohead phiX_base -1 test-data/bowtie_in5.fastqsanger -2 test-data/bowtie_in6.fastqsanger > bowtie_out10_u.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
738 sort bowtie_out10_u.sam > bowtie_out10.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
739 -p is the number of threads. You need to replace the + with 2 dashes.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
740 chrM_base is the index files' location/base name.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
741 -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
742 <param name="genomeSource" value="history" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
743 <param name="ownFile" value="phiX.fasta" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
744 <param name="indexSettings" value="indexFull" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
745 <param name="autoB" value="auto" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
746 <param name="nodc" value="dc" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
747 <param name="noref" value="ref" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
748 <param name="offrate" value="5" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
749 <param name="ftab" value="10" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
750 <param name="ntoa" value="no" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
751 <param name="endian" value="little" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
752 <param name="seed" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
753 <param name="sPaired" value="paired" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
754 <param name="pInput1" ftype="fastqsanger" value="bowtie_in5.fastqsanger" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
755 <param name="pInput2" ftype="fastqsanger" value="bowtie_in6.fastqsanger" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
756 <param name="pMaxInsert" value="1000" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
757 <param name="pMateOrient" value="ff" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
758 <param name="pSettingsType" value="preSet" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
759 <param name="suppressHeader" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
760 <output name="output" ftype="sam" file="bowtie_out10.sam" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
761 </test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
762 <test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
763 <!--
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
764 Bowtie command:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
765 bowtie-build +offrate 5 +ftabchars 10 +little -f test-data/phiX.fasta phiX_base
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
766 bowtie -q -X 1000 +ff -p 4 -S +sam-nohead phiX_base -1 test-data/bowtie_in5.fastqsanger -2 test-data/bowtie_in6.fastqsanger > bowtie_out10_u.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
767 sort bowtie_out10_u.sam > bowtie_out10.sam
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
768 -p is the number of threads. You need to replace the + with 2 dashes.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
769 chrM_base is the index files' location/base name.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
770 -->
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
771 <param name="genomeSource" value="history" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
772 <param name="ownFile" value="phiX.fasta" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
773 <param name="indexSettings" value="indexFull" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
774 <param name="autoB" value="auto" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
775 <param name="nodc" value="dc" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
776 <param name="noref" value="ref" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
777 <param name="offrate" value="5" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
778 <param name="ftab" value="10" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
779 <param name="ntoa" value="no" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
780 <param name="endian" value="little" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
781 <param name="seed" value="-1" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
782 <param name="sPaired" value="paired" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
783 <param name="pInput1" ftype="fastqsanger" value="bowtie_in5.fastqsanger" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
784 <param name="pInput2" ftype="fastqsanger" value="bowtie_in6.fastqsanger" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
785 <param name="pMaxInsert" value="1000" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
786 <param name="pMateOrient" value="ff" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
787 <param name="pSettingsType" value="preSet" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
788 <param name="suppressHeader" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
789 <param name="save_mapping_stats" value="true" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
790 <output name="output" ftype="sam" file="bowtie_out10.sam" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
791 <output name="mapping_stats" ftype="txt" file="bowtie_out11.txt" sort="True" />
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
792 </test>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
793 </tests>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
794
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
795 <help>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
796
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
797 **What it does**
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
798
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
799 Bowtie_ is a short read aligner designed to be ultrafast and memory-efficient. It is developed by Ben Langmead and Cole Trapnell. Please cite: Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biology 10:R25.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
800
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
801 .. _Bowtie: http://bowtie-bio.sourceforge.net/index.shtml
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
802
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
803 ------
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
804
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
805 **Know what you are doing**
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
806
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
807 .. class:: warningmark
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
808
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
809 There is no such thing (yet) as an automated gearshift in short read mapping. 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.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
810
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
811 .. __: http://bowtie-bio.sourceforge.net/index.shtml
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
812
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
813 ------
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
814
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
815 **Input formats**
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
816
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
817 Bowtie accepts files in Sanger FASTQ format. Use the FASTQ Groomer to prepare your files.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
818
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
819 ------
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
820
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
821 **A Note on Built-in Reference Genomes**
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
822
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
823 The default variant for all genomes is "Full", defined as all primary chromosomes (or scaffolds/contigs) including mitochondrial plus associated unmapped, plasmid, and other segments. When only one version of a genome is available in this tool, it represents the default "Full" variant. Some genomes will have more than one variant available. The "Canonical Male" or sometimes simply "Canonical" variant contains the primary chromosomes for a genome. For example a human "Canonical" variant contains chr1-chr22, chrX, chrY, and chrM. The "Canonical Female" variant contains the primary chromosomes excluding chrY.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
824
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
825 ------
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
826
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
827 **Outputs**
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
828
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
829 The output is in SAM format, and has the following columns::
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
830
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
831 Column Description
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
832 -------- --------------------------------------------------------
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
833 1 QNAME Query (pair) NAME
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
834 2 FLAG bitwise FLAG
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
835 3 RNAME Reference sequence NAME
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
836 4 POS 1-based leftmost POSition/coordinate of clipped sequence
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
837 5 MAPQ MAPping Quality (Phred-scaled)
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
838 6 CIGAR extended CIGAR string
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
839 7 MRNM Mate Reference sequence NaMe ('=' if same as RNAME)
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
840 8 MPOS 1-based Mate POSition
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
841 9 ISIZE Inferred insert SIZE
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
842 10 SEQ query SEQuence on the same strand as the reference
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
843 11 QUAL query QUALity (ASCII-33 gives the Phred base quality)
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
844 12 OPT variable OPTional fields in the format TAG:VTYPE:VALUE
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
845
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
846 The flags are as follows::
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
847
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
848 Flag Description
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
849 ------ -------------------------------------
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
850 0x0001 the read is paired in sequencing
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
851 0x0002 the read is mapped in a proper pair
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
852 0x0004 the query sequence itself is unmapped
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
853 0x0008 the mate is unmapped
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
854 0x0010 strand of the query (1 for reverse)
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
855 0x0020 strand of the mate
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
856 0x0040 the read is the first read in a pair
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
857 0x0080 the read is the second read in a pair
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
858 0x0100 the alignment is not primary
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
859
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
860 It looks like this (scroll sideways to see the entire example)::
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
861
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
862 QNAME FLAG RNAME POS MAPQ CIAGR MRNM MPOS ISIZE SEQ QUAL OPT
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
863 HWI-EAS91_1_30788AAXX:1:1:1761:343 4 * 0 0 * * 0 0 AAAAAAANNAAAAAAAAAAAAAAAAAAAAAAAAAAACNNANNGAGTNGNNNNNNNGCTTCCCACAGNNCTGG hhhhhhh;;hhhhhhhhhhh^hOhhhhghhhfhhhgh;;h;;hhhh;h;;;;;;;hhhhhhghhhh;;Phhh
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
864 HWI-EAS91_1_30788AAXX:1:1:1578:331 4 * 0 0 * * 0 0 GTATAGANNAATAAGAAAAAAAAAAATGAAGACTTTCNNANNTCTGNANNNNNNNTCTTTTTTCAGNNGTAG hhhhhhh;;hhhhhhhhhhhhhhhhhhhhhhhhhhhh;;h;;hhhh;h;;;;;;;hhhhhhhhhhh;;hhVh
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
865
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
866 -------
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
867
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
868 **Bowtie settings**
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
869
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
870 All of the options have a default value. You can change any of them. Most of the options in Bowtie have been implemented here.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
871
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
872 ------
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
873
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
874 **Bowtie parameter list**
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
875
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
876 This is an exhaustive list of Bowtie options:
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
877
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
878 For indexing (bowtie-build)::
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
879
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
880 -a No auto behavior. Disable the default behavior where bowtie automatically
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
881 selects values for --bmax/--bmaxdivn/--dcv/--packed parameters according
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
882 to the memory available. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
883 --packed Packing. Use a packed representation for DNA strings. [auto]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
884 --bmax INT Suffix maximum. The maximum number of suffixes allowed in a block. [auto]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
885 --bmaxdivn INT Suffix maximum fraction. The maximum number of suffixes allowed in a block
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
886 expressed as a fraction of the length of the reference. [4]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
887 --dcv INT Difference-cover sample. Use INT as the period for the difference-cover
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
888 sample. [1024]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
889 --nodc INT No difference-cover sample. Disable the difference-cover sample. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
890 -r No reference indexes. Do not build the NAME.3.ebwt and NAME.4.ebwt portions
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
891 of the index. Used only for paired-end alignment. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
892 -o Offrate. How many Burrows-Wheeler rows get marked by the indexer. The
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
893 indexer will mark every 2^INT rows. The marked rows correspond to rows on
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
894 the genome. [5]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
895 -t INT The ftab lookup table used to calculate an initial Burrows-Wheeler range
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
896 with respect to the first INT characters of the query. Ftab size is 4^(INT+1)
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
897 bytes. [10]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
898 --ntoa N conversion. Convert Ns to As before building the index. Otherwise, Ns are
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
899 simply excluded from the index and Bowtie will not find alignments that
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
900 overlap them. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
901 --big Endianness. Endianness to use when serializing integers to the index file. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
902 --little Endianness. [--little]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
903 --seed INT Random seed. Use INT as the seed for the pseudo-random number generator. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
904
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
905 For aligning (bowtie)::
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
906
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
907 -s INT Skip. Do not align the first INT reads or pairs in the input. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
908 -u INT Align limit. Only align the first INT reads/pairs from the input. [no limit]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
909 -5 INT High-quality trim. Trim INT bases from the high-quality (left) end of each
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
910 read before alignment. [0]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
911 -3 INT Low-quality trim. Trim INT bases from the low-quality (right) end of each
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
912 read before alignment. [0]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
913 -n INT Mismatch seed. Maximum number of mismatches permitted in the seed (defined
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
914 with seed length option). Can be 0, 1, 2, or 3. [2]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
915 -e INT Mismatch quality. Maximum permitted total of quality values at mismatched
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
916 read positions. Bowtie rounds quality values to the nearest 10 and saturates
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
917 at 30. [70]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
918 -l INT Seed length. The number of bases on the high-quality end of the read to
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
919 which the -n ceiling applies. Must be at least 5. [28]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
920 --nomaqround Suppress Maq rounding. Values are internally rounded to the nearest 10 and
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
921 saturate at 30. This options turns off that rounding. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
922 -v INT Maq- or SOAP-like alignment policy. This option turns off the default
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
923 Maq-like alignment policy in favor of a SOAP-like one. End-to-end alignments
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
924 with at most INT mismatches. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
925 -I INT Minimum insert. The minimum insert size for valid paired-end alignments.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
926 Does checking on untrimmed reads if -5 or -3 is used. [0]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
927 -X INT Maximum insert. The maximum insert size for valid paired-end alignments.
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
928 Does checking on untrimmed reads if -5 or -3 is used. [250]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
929 --fr Mate orientation. The upstream/downstream mate orientations for a valid
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
930 paired-end alignment against the forward reference strand. [--fr]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
931 --rf Mate orientation. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
932 --ff Mate orientation. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
933 --pairtries INT Maximum alignment attempts for paired-end data. [100]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
934 --nofw No forward aligning. Choosing this option means that Bowtie will not attempt
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
935 to align against the forward reference strand. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
936 --norc No reverse-complement aligning. Setting this will mean that Bowtie will not
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
937 attempt to align against the reverse-complement reference strand. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
938 --un FILENAME Write all reads that could not be aligned to file [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
939 --max FILENAME Write all reads with a number of valid alignments exceeding the limit
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
940 set with the -m option to file [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
941 --maxbts INT Maximum backtracks. The maximum number of backtracks permitted when aligning
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
942 a read in -n 2 or -n 3 mode. [125 without --best] [800 with --best]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
943 -y Try hard. Try as hard as possible to find valid alignments when they exist,
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
944 including paired-end alignments. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
945 --chunkmbs INT Thread memory. The number of megabytes of memory a given thread is given to
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
946 store path descriptors in --best mode. [32]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
947 -k INT Valid alignments. The number of valid alignments per read or pair. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
948 -a All valid alignments. Choosing this means that all valid alignments per read
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
949 or pair will be reported. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
950 -m INT Suppress alignments. Suppress all alignments for a particular read or pair
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
951 if more than INT reportable alignments exist for it. [no limit]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
952 --best Best mode. Make Bowtie guarantee that reported singleton alignments are
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
953 "best" in terms of stratum (the number of mismatches) and quality values at
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
954 mismatched position. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
955 --strata Best strata. When running in best mode, report alignments that fall into the
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
956 best stratum if there are ones falling into more than one. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
957 -o INT Offrate override. Override the offrate of the index with INT. Some row
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
958 markings are discarded when index read into memory. INT must be greater than
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
959 the value used to build the index (default: 5). [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
960 --seed INT Random seed. Use INT as the seed for the pseudo-random number generator. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
961 --snpphred INT Use INT as the SNP penalty for decoding colorspace alignments. True ratio of
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
962 SNPs per base in the subject genome. [see --snpfrac]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
963 --snpfrac DEC Use DEC as the estimated ratio of SNPs per base when decoding colorspace
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
964 alignments. [0.001]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
965 --col-keepends Keep the extreme-end nucleotides and qualities when decoding colorspace
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
966 alignments. [off]
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
967
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
968 </help>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
969 <citations>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
970 <citation type="doi">10.1186/gb-2009-10-3-r25</citation>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
971 </citations>
1b109d213a42 Uploaded
jackcurragh
parents:
diff changeset
972 </tool>