annotate bsmap.xml @ 33:91a4092971bc draft default tip

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