annotate bsmap.xml @ 15:fc7ffd15b17f draft

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