Mercurial > repos > eiriche > bsmap
comparison bsmap_meth_caller.xml @ 26:701d6bb0e28d draft
Uploaded
| author | eiriche |
|---|---|
| date | Mon, 03 Dec 2012 02:44:43 -0500 |
| parents | |
| children |
comparison
equal
deleted
inserted
replaced
| 25:4c6d4be21daf | 26:701d6bb0e28d |
|---|---|
| 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 |