annotate bsmap_meth_caller.xml @ 7:be88d0f3f6f2 draft

Uploaded
author eiriche
date Thu, 29 Nov 2012 10:14:57 -0500
parents e5025e0378d0
children
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
4
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
1 <tool id="bsmap_meth_caller" name="BSMAP Methylation Caller">
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
2 <requirements>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
3 <requirement type='package'>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
4 bsmap
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
5 </requirement>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
6 </requirements>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
7 <requirements>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
8 <requirement type='package'>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
9 samtools
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
10 </requirement>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
11 </requirements>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
12 <command interpreter="bash">
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
13 bsmap_meth_caller.sh
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
14 input=$bsmap_sam
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
15 unique=$unique
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
16 properly=$properly
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
17 zero_meth = $zero_meth
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
18 rem_dup = $rem_dup
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
19 combine_cpg = $combine_cpg
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
20 trimN = $trimN
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
21 depth = $depth
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
22 output=$output
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
23 tempdir=$output.files_path
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
24 ref="${ filter( lambda x: str( x[1] ) == str( $bsmap_sam.metadata.dbkey ), $__app__.tool_data_tables['bsmap_fasta'].get_fields() )[0][3] }"
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
25 </command>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
26 <inputs>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
27 <param name="bsmap_sam" format="sam" type="data" label="BSMAP mapping output file" help="Must be in SAM format" />
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
28 <param name="unique" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Process only unique mappings/pairs" />
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
29 <param name="properly" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Process only properly paired mappings" />
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
30 <param name="zero_meth" type="boolean" truevalue="true" falsevalue="false" checked="True" label="report loci with zero methylation ratios" />
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
31 <param name="rem_dup" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Remove duplicated reads" />
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
32 <param name="combine_cpg" type="boolean" truevalue="true" falsevalue="false" checked="False" label="Combine CpG methylaion ratios on both strands" />
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
33 <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" />
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
34 <param name="depth" type="integer" value="1" label="Minimum sequencing depth to report loci" />
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
35 </inputs>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
36 <outputs>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
37 <data name="output" format ="bed" label="BSMAP methylation output" />
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
38 </outputs>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
39 <help>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
40 **What it does**
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
41
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
42 This methylation caller parses the BSMAP SAM output file into bed format.
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
43
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
44
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
45 **Output format** ::
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
46
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
47
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
48 Column Description
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
49 ---------------------- --------------------------------------
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
50 1 chr chromosome
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
51 2 pos position
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
52 3 strand strand
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
53 4 context context (CHH,CHG,CpG)
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
54 5 coverage totally sequenced Cs at that position
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
55 6 methylated methylated Cs at that position
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
56 7 percentage methylated percentage of 6
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
57 </help>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
58 </tool>
e5025e0378d0 Uploaded
eiriche
parents:
diff changeset
59