annotate bsmap.xml @ 3:91e88de226a3 draft

Uploaded
author eiriche
date Thu, 29 Nov 2012 10:10:16 -0500
parents b35020882aad
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
0
b35020882aad Uploaded
eiriche
parents:
diff changeset
1 <tool id="bsmap" name="BSMAP Mapper">
b35020882aad Uploaded
eiriche
parents:
diff changeset
2 <requirements>
b35020882aad Uploaded
eiriche
parents:
diff changeset
3 <requirement type='package'>
b35020882aad Uploaded
eiriche
parents:
diff changeset
4 bsmap
b35020882aad Uploaded
eiriche
parents:
diff changeset
5 </requirement>
b35020882aad Uploaded
eiriche
parents:
diff changeset
6 </requirements>
b35020882aad Uploaded
eiriche
parents:
diff changeset
7 <command interpreter="bash">
b35020882aad Uploaded
eiriche
parents:
diff changeset
8 bsmap_wrapper.sh
b35020882aad Uploaded
eiriche
parents:
diff changeset
9 ##Reference genome
b35020882aad Uploaded
eiriche
parents:
diff changeset
10 ref="${reference.fields.path}"
b35020882aad Uploaded
eiriche
parents:
diff changeset
11 ##Output files (SAM output, BSMAP summary)
b35020882aad Uploaded
eiriche
parents:
diff changeset
12 mapped=$mapped
b35020882aad Uploaded
eiriche
parents:
diff changeset
13 ##Temp directory
b35020882aad Uploaded
eiriche
parents:
diff changeset
14 tempdir=$mapped.files_path
b35020882aad Uploaded
eiriche
parents:
diff changeset
15 summary=$summary
b35020882aad Uploaded
eiriche
parents:
diff changeset
16 #if str($singlePaired.sPaired) == "single":
b35020882aad Uploaded
eiriche
parents:
diff changeset
17 library="single"
b35020882aad Uploaded
eiriche
parents:
diff changeset
18 mate1=$singlePaired.sInput1
b35020882aad Uploaded
eiriche
parents:
diff changeset
19 #if str($singlePaired.sParams.sSettingsType) == "full":
b35020882aad Uploaded
eiriche
parents:
diff changeset
20 fullparam=true
b35020882aad Uploaded
eiriche
parents:
diff changeset
21 qual=$singlePaired.sParams.qual
b35020882aad Uploaded
eiriche
parents:
diff changeset
22 threshold=$singlePaired.sParams.threshold
b35020882aad Uploaded
eiriche
parents:
diff changeset
23 lowqual=$singlePaired.sParams.lowqual
b35020882aad Uploaded
eiriche
parents:
diff changeset
24 adapter=$singlePaired.sParams.adapter
b35020882aad Uploaded
eiriche
parents:
diff changeset
25 firstn=$singlePaired.sParams.firstn
b35020882aad Uploaded
eiriche
parents:
diff changeset
26 repeat_reads=$singlePaired.sParams.repeat_reads
b35020882aad Uploaded
eiriche
parents:
diff changeset
27 seed_size=$singlePaired.sParams.seed_size
b35020882aad Uploaded
eiriche
parents:
diff changeset
28 mismatch=$singlePaired.sParams.mismatch
b35020882aad Uploaded
eiriche
parents:
diff changeset
29 equal_best=$singlePaired.sParams.equal_best
b35020882aad Uploaded
eiriche
parents:
diff changeset
30 start=$singlePaired.sParams.start
b35020882aad Uploaded
eiriche
parents:
diff changeset
31 end=$singlePaired.sParams.end
b35020882aad Uploaded
eiriche
parents:
diff changeset
32 index_interval=$singlePaired.sParams.index_interval
b35020882aad Uploaded
eiriche
parents:
diff changeset
33 seed_random=$singlePaired.sParams.seed_random
b35020882aad Uploaded
eiriche
parents:
diff changeset
34 rrbs=$singlePaired.sParams.rrbs
b35020882aad Uploaded
eiriche
parents:
diff changeset
35 mode=$singlePaired.sParams.mode
b35020882aad Uploaded
eiriche
parents:
diff changeset
36 align_info=$singlePaired.sParams.align_info
b35020882aad Uploaded
eiriche
parents:
diff changeset
37 #end if
b35020882aad Uploaded
eiriche
parents:
diff changeset
38 #else:
b35020882aad Uploaded
eiriche
parents:
diff changeset
39 library="paired"
b35020882aad Uploaded
eiriche
parents:
diff changeset
40 mate1=$singlePaired.pInput1
b35020882aad Uploaded
eiriche
parents:
diff changeset
41 mate2=$singlePaired.pInput2
b35020882aad Uploaded
eiriche
parents:
diff changeset
42 unpaired=$unpaired
b35020882aad Uploaded
eiriche
parents:
diff changeset
43 #if str($singlePaired.pParams.pSettingsType) == "full":
b35020882aad Uploaded
eiriche
parents:
diff changeset
44 fullparam=true
b35020882aad Uploaded
eiriche
parents:
diff changeset
45 qual=$singlePaired.pParams.qual
b35020882aad Uploaded
eiriche
parents:
diff changeset
46 threshold=$singlePaired.pParams.threshold
b35020882aad Uploaded
eiriche
parents:
diff changeset
47 lowqual=$singlePaired.pParams.lowqual
b35020882aad Uploaded
eiriche
parents:
diff changeset
48 adapter=$singlePaired.pParams.adapter
b35020882aad Uploaded
eiriche
parents:
diff changeset
49 firstn=$singlePaired.pParams.firstn
b35020882aad Uploaded
eiriche
parents:
diff changeset
50 repeat_reads=$singlePaired.pParams.repeat_reads
b35020882aad Uploaded
eiriche
parents:
diff changeset
51 seed_size=$singlePaired.pParams.seed_size
b35020882aad Uploaded
eiriche
parents:
diff changeset
52 mismatch=$singlePaired.pParams.mismatch
b35020882aad Uploaded
eiriche
parents:
diff changeset
53 equal_best=$singlePaired.pParams.equal_best
b35020882aad Uploaded
eiriche
parents:
diff changeset
54 start=$singlePaired.pParams.start
b35020882aad Uploaded
eiriche
parents:
diff changeset
55 end=$singlePaired.pParams.end
b35020882aad Uploaded
eiriche
parents:
diff changeset
56 index_interval=$singlePaired.pParams.index_interval
b35020882aad Uploaded
eiriche
parents:
diff changeset
57 seed_random=$singlePaired.pParams.seed_random
b35020882aad Uploaded
eiriche
parents:
diff changeset
58 rrbs=$singlePaired.pParams.rrbs
b35020882aad Uploaded
eiriche
parents:
diff changeset
59 mode=$singlePaired.pParams.mode
b35020882aad Uploaded
eiriche
parents:
diff changeset
60 align_info=$singlePaired.pParams.align_info
b35020882aad Uploaded
eiriche
parents:
diff changeset
61 maxinsert=$singlePaired.pParams.maxinsert
b35020882aad Uploaded
eiriche
parents:
diff changeset
62 mininsert=$singlePaired.pParams.mininsert
b35020882aad Uploaded
eiriche
parents:
diff changeset
63 #end if
b35020882aad Uploaded
eiriche
parents:
diff changeset
64 #end if
b35020882aad Uploaded
eiriche
parents:
diff changeset
65 </command>
b35020882aad Uploaded
eiriche
parents:
diff changeset
66 <inputs>
b35020882aad Uploaded
eiriche
parents:
diff changeset
67 <param name="reference" type="select" label="Select a reference genome">
b35020882aad Uploaded
eiriche
parents:
diff changeset
68 <options from_data_table="all_fasta">
b35020882aad Uploaded
eiriche
parents:
diff changeset
69 <filter type="sort_by" column="2" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
70 <validator type="no_options" message="No reference genomes are available" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
71 </options>
b35020882aad Uploaded
eiriche
parents:
diff changeset
72 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
73
b35020882aad Uploaded
eiriche
parents:
diff changeset
74 <conditional name="singlePaired">
b35020882aad Uploaded
eiriche
parents:
diff changeset
75 <param name="sPaired" type="select" label="Is this library mate-paired?">
b35020882aad Uploaded
eiriche
parents:
diff changeset
76 <option value="single">Single-end</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
77 <option value="paired">Paired-end</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
78 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
79 <when value="single">
b35020882aad Uploaded
eiriche
parents:
diff changeset
80 <param name="sInput1" type="data" format="fastq,fasta" label="FASTQ file" help="Must have ASCII encoded quality scores"/>
b35020882aad Uploaded
eiriche
parents:
diff changeset
81 <conditional name="sParams">
b35020882aad Uploaded
eiriche
parents:
diff changeset
82 <param name="sSettingsType" type="select" label="BSMAP settings to use" help="For most mapping needs use Commonly used settings. If you want full control use Full parameter list">
b35020882aad Uploaded
eiriche
parents:
diff changeset
83 <option value="preSet">Commonly used</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
84 <option value="full">Full parameter list</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
85 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
86 <when value="preSet" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
87 <when value="full">
b35020882aad Uploaded
eiriche
parents:
diff changeset
88 <param name="qual" type="select" label="Select the type of FastQ qualities">
b35020882aad Uploaded
eiriche
parents:
diff changeset
89 <option value="33">phred33-quals</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
90 <option value="64">phred64-quals</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
91 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
92 <param name="threshold" type="integer" value="0" label="Quality threshold in trimming" help="0-40, default=0 (no trim)" min="0" max="40" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
93 <param name="lowqual" type="integer" value="5" label="Filter low-quality reads containing >n Ns" help="default=5" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
94 <param name="adapter" type="text" value="none" label="3-end adapter sequence" help="default: none (no trim)" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
95 <param name="firstn" type="integer" value="144" label="Map the first N nucleotides of the read" help="default:144 (map the whole read)" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
96
b35020882aad Uploaded
eiriche
parents:
diff changeset
97 <param name="repeat_reads" type="select" label="How to report repeat hits">
b35020882aad Uploaded
eiriche
parents:
diff changeset
98 <option value="0">none(unique hit only)</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
99 <option value="1">random one</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
100 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
101
b35020882aad Uploaded
eiriche
parents:
diff changeset
102 <param name="seed_size" type="integer" value="16" label="Seed size" min="8" max="16" help="default=16(WGBS mode), 12(RRBS mode)" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
103 <param name="mismatch" type="integer" value="2" label="Maximum number of mismatches allowed on a read" max="15" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
104 <param name="equal_best" type="integer" value="20" label="Maximum number of equal best hits to count" max="1000" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
105 <param name="start" type="integer" value="1" label="Start from the Nth read or read pair" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
106 <param name="end" type="integer" value="4294967295" label="End at the Nth read or read pair" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
107 <param name="index_interval" type="integer" value="4" label="Index interval" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
108 <param name="seed_random" type="integer" value="-1" label="Seed for random number generation used in selecting multiple hits" help="other seed values generate pseudo random number based on read index number, to allow reproducible mapping results" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
109 <param name="rrbs" type="text" value="none" label="Activating RRBS mapping mode and set restriction enzyme digestion sites" help="digestion position marked by '-', example: -D C-CGG for MspI digestion. default: none (whole genome shotgun bisulfite mapping mode)" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
110 <param name="mode" type="select" label="Set mapping strand information">
b35020882aad Uploaded
eiriche
parents:
diff changeset
111 <option value="0">only map to 2 forward strands</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
112 <option value="1">map SE or PE reads to all 4 strands</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
113 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
114 <param name="align_info" type="text" value="none" label="Set alignment information for the additional nucleotide transition" help="is in the form of two different nucleotides N1N2,indicating N1 in the reads could be mapped to N2 in the reference sequences. default: -M TC, corresponds to C=>U(T) transition in bisulfite conversion. example: -M GA could be used to detect A=>I(G) transition in RNA editing." />
b35020882aad Uploaded
eiriche
parents:
diff changeset
115 </when> <!-- full -->
b35020882aad Uploaded
eiriche
parents:
diff changeset
116 </conditional> <!-- sParams -->
b35020882aad Uploaded
eiriche
parents:
diff changeset
117 </when> <!-- single -->
b35020882aad Uploaded
eiriche
parents:
diff changeset
118
b35020882aad Uploaded
eiriche
parents:
diff changeset
119 <when value="paired">
b35020882aad Uploaded
eiriche
parents:
diff changeset
120 <param name="pInput1" type="data" format="fastq,fasta" label="Forward FASTQ file" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
121 <param name="pInput2" type="data" format="fastq,fasta" label="Reverse FASTQ file" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
122
b35020882aad Uploaded
eiriche
parents:
diff changeset
123 <conditional name="pParams">
b35020882aad Uploaded
eiriche
parents:
diff changeset
124 <param name="pSettingsType" type="select" label="BSMAP settings to use" help="For most mapping needs use Commonly used settings. If you want full control use Full parameter list">
b35020882aad Uploaded
eiriche
parents:
diff changeset
125 <option value="preSet">Commonly used</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
126 <option value="full">Full parameter list</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
127 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
128 <when value="preSet" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
129 <when value="full">
b35020882aad Uploaded
eiriche
parents:
diff changeset
130 <param name="qual" type="select" label="Select the type of FastQ qualities">
b35020882aad Uploaded
eiriche
parents:
diff changeset
131 <option value="33">phred33-quals</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
132 <option value="64">phred64-quals</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
133 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
134
b35020882aad Uploaded
eiriche
parents:
diff changeset
135 <param name="mininsert" type="integer" value="28" label="Minimal insert size allowed" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
136 <param name="maxinsert" type="integer" value="500" label="Maximal insert size allowed" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
137
b35020882aad Uploaded
eiriche
parents:
diff changeset
138 <param name="threshold" type="integer" value="0" label="Quality threshold in trimming" help="0-40, default=0 (no trim)" min="0" max="40" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
139 <param name="lowqual" type="integer" value="5" label="Filter low-quality reads containing >n Ns" help="default=5" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
140 <param name="adapter" type="text" value="none" label="3-end adapter sequence" help="default: none (no trim)" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
141 <param name="firstn" type="integer" value="144" label="Map the first N nucleotides of the read" help="default:144 (map the whole read)" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
142
b35020882aad Uploaded
eiriche
parents:
diff changeset
143 <param name="repeat_reads" type="select" label="How to report repeat hits">
b35020882aad Uploaded
eiriche
parents:
diff changeset
144 <option value="0">none(unique hit only)</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
145 <option value="1">random one</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
146 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
147
b35020882aad Uploaded
eiriche
parents:
diff changeset
148 <param name="seed_size" type="integer" value="16" label="Seed size" min="8" max="16" help="default=16(WGBS mode), 12(RRBS mode)" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
149 <param name="mismatch" type="integer" value="2" label="Maximum number of mismatches allowed on a read" max="15" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
150 <param name="equal_best" type="integer" value="20" label="Maximum number of equal best hits to count" max="1000" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
151 <param name="start" type="integer" value="1" label="Start from the Nth read or read pair" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
152 <param name="end" type="integer" value="4294967295" label="End at the Nth read or read pair" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
153 <param name="index_interval" type="integer" value="4" label="Index interval" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
154 <param name="seed_random" type="integer" value="-1" label="Seed for random number generation used in selecting multiple hits" help="other seed values generate pseudo random number based on read index number, to allow reproducible mapping results" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
155 <param name="rrbs" type="text" value="none" label="Activating RRBS mapping mode and set restriction enzyme digestion sites" help="digestion position marked by '-', example: -D C-CGG for MspI digestion. default: none (whole genome shotgun bisulfite mapping mode)" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
156 <param name="mode" type="select" label="Set mapping strand information">
b35020882aad Uploaded
eiriche
parents:
diff changeset
157 <option value="0">only map to 2 forward strands</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
158 <option value="1">map SE or PE reads to all 4 strands</option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
159 </param>
b35020882aad Uploaded
eiriche
parents:
diff changeset
160 <param name="align_info" type="text" value="none" label="Set alignment information for the additional nucleotide transition" help="is in the form of two different nucleotides N1N2,indicating N1 in the reads could be mapped to N2 in the reference sequences. default: -M TC, corresponds to C=>U(T) transition in bisulfite conversion. example: -M GA could be used to detect A=>I(G) transition in RNA editing." />
b35020882aad Uploaded
eiriche
parents:
diff changeset
161
b35020882aad Uploaded
eiriche
parents:
diff changeset
162
b35020882aad Uploaded
eiriche
parents:
diff changeset
163 </when> <!-- full -->
b35020882aad Uploaded
eiriche
parents:
diff changeset
164 </conditional> <!-- pParams -->
b35020882aad Uploaded
eiriche
parents:
diff changeset
165 </when> <!-- paired -->
b35020882aad Uploaded
eiriche
parents:
diff changeset
166 </conditional> <!-- singlePaired -->
b35020882aad Uploaded
eiriche
parents:
diff changeset
167
b35020882aad Uploaded
eiriche
parents:
diff changeset
168
b35020882aad Uploaded
eiriche
parents:
diff changeset
169 </inputs>
b35020882aad Uploaded
eiriche
parents:
diff changeset
170 <outputs>
b35020882aad Uploaded
eiriche
parents:
diff changeset
171 <data name="mapped" format="sam" label="BSMAP Mapped Reads">
b35020882aad Uploaded
eiriche
parents:
diff changeset
172 <actions>
b35020882aad Uploaded
eiriche
parents:
diff changeset
173 <action type="metadata" name="dbkey">
b35020882aad Uploaded
eiriche
parents:
diff changeset
174 <option type="from_data_table" name="bsmap_fasta" column="1" offset="0">
b35020882aad Uploaded
eiriche
parents:
diff changeset
175 <filter type="param_value" column="0" value="#" compare="startswith" keep="False"/>
b35020882aad Uploaded
eiriche
parents:
diff changeset
176 <filter type="param_value" ref="reference" column="0"/>
b35020882aad Uploaded
eiriche
parents:
diff changeset
177 </option>
b35020882aad Uploaded
eiriche
parents:
diff changeset
178 </action>
b35020882aad Uploaded
eiriche
parents:
diff changeset
179 </actions>
b35020882aad Uploaded
eiriche
parents:
diff changeset
180 </data>
b35020882aad Uploaded
eiriche
parents:
diff changeset
181 <data name="summary" format="txt" label="BSMAP Mapping Summary" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
182 <data name="unpaired" format ="sam" label="BSMAP Unpaired Hits">
b35020882aad Uploaded
eiriche
parents:
diff changeset
183 <filter>(singlePaired['sPaired'] == 'paired')</filter>
b35020882aad Uploaded
eiriche
parents:
diff changeset
184 </data>
b35020882aad Uploaded
eiriche
parents:
diff changeset
185
b35020882aad Uploaded
eiriche
parents:
diff changeset
186 </outputs>
b35020882aad Uploaded
eiriche
parents:
diff changeset
187 <help>
b35020882aad Uploaded
eiriche
parents:
diff changeset
188 **What it does**
b35020882aad Uploaded
eiriche
parents:
diff changeset
189
b35020882aad Uploaded
eiriche
parents:
diff changeset
190 BSMAP is a short reads mapping software for bisulfite sequencing reads. It has the following features:
b35020882aad Uploaded
eiriche
parents:
diff changeset
191
b35020882aad Uploaded
eiriche
parents:
diff changeset
192 - read length up to 144 nt, allow up to 15 mismatches, gap size up to 3 bp.
b35020882aad Uploaded
eiriche
parents:
diff changeset
193
b35020882aad Uploaded
eiriche
parents:
diff changeset
194 - support single end and pair end mapping. support multi-thread mapping.
b35020882aad Uploaded
eiriche
parents:
diff changeset
195
b35020882aad Uploaded
eiriche
parents:
diff changeset
196 - support both "Lister protocol" (sequence 2 forward strands only) and "Cokus protocol" (sequence all 4 bisulfite converted strands)
b35020882aad Uploaded
eiriche
parents:
diff changeset
197
b35020882aad Uploaded
eiriche
parents:
diff changeset
198 - reads are directly mapped to original reference genome sequence, no need to preprocess the reads and reference genome to convert C to T.
b35020882aad Uploaded
eiriche
parents:
diff changeset
199
b35020882aad Uploaded
eiriche
parents:
diff changeset
200 - support both whole genome bisulfite sequencing (WGBS) mode and reduced representation bisulfite sequencing (RRBS) mode, allow changing the digestion site information to support different digestion enzymes for RRBS.
b35020882aad Uploaded
eiriche
parents:
diff changeset
201
b35020882aad Uploaded
eiriche
parents:
diff changeset
202 - allow trimming adapter sequences and low quality nucleotides from the 3'end of reads
b35020882aad Uploaded
eiriche
parents:
diff changeset
203
b35020882aad Uploaded
eiriche
parents:
diff changeset
204 - allow trade off between speed/memory usage/mapping sensitivity. For human genome, the RRBS mode uses ~3GB. In WGBS mode, the typical memory usage is ~9GB, but can be as low as 5GB.
b35020882aad Uploaded
eiriche
parents:
diff changeset
205
b35020882aad Uploaded
eiriche
parents:
diff changeset
206 - allow alignment for other nucleotide transitions, for example, can be set to detect the A=>I(G) transition in RNA editing.
b35020882aad Uploaded
eiriche
parents:
diff changeset
207
b35020882aad Uploaded
eiriche
parents:
diff changeset
208 .. _BSMAP: http://code.google.com/p/bsmap/
b35020882aad Uploaded
eiriche
parents:
diff changeset
209
b35020882aad Uploaded
eiriche
parents:
diff changeset
210 **Input formats**
b35020882aad Uploaded
eiriche
parents:
diff changeset
211
b35020882aad Uploaded
eiriche
parents:
diff changeset
212 BSMAP accepts files in FASTA/FASTQ format.
b35020882aad Uploaded
eiriche
parents:
diff changeset
213
b35020882aad Uploaded
eiriche
parents:
diff changeset
214 **Outputs**
b35020882aad Uploaded
eiriche
parents:
diff changeset
215
b35020882aad Uploaded
eiriche
parents:
diff changeset
216 The output contains the following files:
b35020882aad Uploaded
eiriche
parents:
diff changeset
217
b35020882aad Uploaded
eiriche
parents:
diff changeset
218 - mapped reads in SAM format
b35020882aad Uploaded
eiriche
parents:
diff changeset
219
b35020882aad Uploaded
eiriche
parents:
diff changeset
220 - mapping summary
b35020882aad Uploaded
eiriche
parents:
diff changeset
221
b35020882aad Uploaded
eiriche
parents:
diff changeset
222 - unpaired hits (only for paired-end mapping)
b35020882aad Uploaded
eiriche
parents:
diff changeset
223
b35020882aad Uploaded
eiriche
parents:
diff changeset
224 </help>
b35020882aad Uploaded
eiriche
parents:
diff changeset
225
b35020882aad Uploaded
eiriche
parents:
diff changeset
226 <tests>
b35020882aad Uploaded
eiriche
parents:
diff changeset
227 <test>
b35020882aad Uploaded
eiriche
parents:
diff changeset
228 <param name="sPaired" value="single" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
229 <param name="indices" value="bismark" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
230 <param name="sInput1" ftype="fastq" value="bismark_test_single.fastq" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
231 <param name="sParams" value="preSet" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
232 <output name="mapped" ftype="SAM" file="bismark_result_single_1.SAM" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
233 <ouput name="summary" ftype="txt" file="bismark_result_single_2.txt" />
b35020882aad Uploaded
eiriche
parents:
diff changeset
234 </test>
b35020882aad Uploaded
eiriche
parents:
diff changeset
235 </tests>
b35020882aad Uploaded
eiriche
parents:
diff changeset
236 </tool>
b35020882aad Uploaded
eiriche
parents:
diff changeset
237