|
24
|
1 <tool id="bsmap_meth_caller" name="BSMAP Methylation Caller">
|
|
|
2 <requirements>
|
|
|
3 <requirement type='package' version="2.6">bsmap</requirement>
|
|
|
4 <requirement type='package'>samtools</requirement>
|
|
|
5 <requirements>
|
|
|
6 <command interpreter="bash">
|
|
|
7 bsmap_meth_caller.sh
|
|
|
8 input=$bsmap_sam
|
|
|
9 unique=$unique
|
|
|
10 properly=$properly
|
|
|
11 zero_meth = $zero_meth
|
|
|
12 rem_dup = $rem_dup
|
|
|
13 combine_cpg = $combine_cpg
|
|
|
14 trimN = $trimN
|
|
|
15 depth = $depth
|
|
|
16 output=$output
|
|
|
17 tempdir=$output.files_path
|
|
|
18 #if $refGenomeSource.genomeSource == "history":
|
|
|
19 ref="${refGenomeSource.myFile}"
|
|
|
20 #else
|
|
|
21 ref="${refGenomeSource.builtinFile.fields.path}"
|
|
|
22 #end if
|
|
|
23
|
|
|
24 </command>
|
|
|
25 <inputs>
|
|
|
26 <conditional name="refGenomeSource">
|
|
|
27 <param name="genomeSource" type="select" label="Will you select a reference genome from your history or use a built-in reference?" help="Must be the same reference as used for the mapping">
|
|
|
28 <option value="builtin">Use a built-in reference</option>
|
|
|
29 <option value="history">Use one from the history</option>
|
|
|
30 </param>
|
|
|
31 <when value="builtin">
|
|
|
32 <param name="builtinFile" type="select" label="Select the reference genome">
|
|
|
33 <options from_data_table="bsmap_fasta">
|
|
|
34 <filter type="sort_by" column="2" />
|
|
|
35 <validator type="no_options" message="No reference genomes are available" />
|
|
|
36 </options>
|
|
|
37 </param>
|
|
|
38 </when>
|
|
|
39 <when value="history">
|
|
|
40 <param name="myFile" type="data" format="fasta" label="Select the reference genome" />
|
|
|
41 </when>
|
|
|
42 </conditional>
|
|
|
43
|
|
|
44 <param name="bsmap_sam" format="sam" type="data" label="BSMAP mapping output file" help="Must be in SAM format" />
|
|
|
45 <param name="unique" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Process only unique mappings/pairs" />
|
|
|
46 <param name="properly" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Process only properly paired mappings" />
|
|
|
47 <param name="zero_meth" type="boolean" truevalue="true" falsevalue="false" checked="True" label="report loci with zero methylation ratios" />
|
|
|
48 <param name="rem_dup" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Remove duplicated reads" />
|
|
|
49 <param name="combine_cpg" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Combine CpG methylaion ratios on both strands" />
|
|
|
50 <param name="trimN" type="integer" value="2" label="Trim N fill-in nucleotides in DNA fragment end-repairing" help="This option is only for pair-end mapping. For RRBS, N could be detetmined by the distance between cuttings sites on forward and reverse strands. For WGBS, N is usually between 0~3" />
|
|
|
51 <param name="depth" type="integer" value="1" label="Minimum sequencing depth to report loci" />
|
|
|
52 </inputs>
|
|
|
53 <outputs>
|
|
|
54 <data name="output" format ="bed" label="BSMAP methylation output" />
|
|
|
55 </outputs>
|
|
|
56 <help>
|
|
|
57 **What it does**
|
|
|
58
|
|
|
59 This methylation caller parses the BSMAP SAM output file into bed format.
|
|
|
60
|
|
|
61
|
|
|
62 **Output format** ::
|
|
|
63
|
|
|
64
|
|
|
65 Column Description
|
|
|
66 ---------------------- --------------------------------------
|
|
|
67 1 chr chromosome
|
|
|
68 2 pos position
|
|
|
69 3 strand strand
|
|
|
70 4 context context (CHH,CHG,CpG)
|
|
|
71 5 coverage totally sequenced Cs at that position
|
|
|
72 6 methylated methylated Cs at that position
|
|
|
73 7 percentage methylated percentage of 6
|
|
|
74 </help>
|
|
|
75 </tool>
|
|
|
76
|
