# HG changeset patch # User eiriche # Date 1354520498 18000 # Node ID b5a8c2ecef59e1013b7b2bb78afd00d97e2dec53 # Parent e3e36b45c1dac831d7c7af64049d90dbdd57e19c Deleted selected files diff -r e3e36b45c1da -r b5a8c2ecef59 bsmap_meth_caller.sh --- a/bsmap_meth_caller.sh Mon Dec 03 02:38:26 2012 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,32 +0,0 @@ -#!/bin/bash -# -# Galaxy wrapper for BSMAP Methylation Caller -# - -set -e - -#get parameters - -until [ $# -eq 0 ] -do - case $1 in - input=*) - input=${1#input=} - ;; - method=*) - method=${1#method=} - ;; - output=*) - output=${1#output=} - ;; - tempdir=*) - tempdir=${1#tempdir=} - ;; - ref=*) - ref=${1#ref=} - ;; - esac - shift -done - -methratio.py -o $output -d $ref -q $input \ No newline at end of file diff -r e3e36b45c1da -r b5a8c2ecef59 bsmap_meth_caller.xml --- a/bsmap_meth_caller.xml Mon Dec 03 02:38:26 2012 -0500 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,75 +0,0 @@ - - - bsmap - samtools - - bsmap_meth_caller.sh - input=$bsmap_sam - unique=$unique - properly=$properly - zero_meth = $zero_meth - rem_dup = $rem_dup - combine_cpg = $combine_cpg - trimN = $trimN - depth = $depth - output=$output - tempdir=$output.files_path - #if $refGenomeSource.genomeSource == "history": - ref="${refGenomeSource.myFile}" - #else - ref="${refGenomeSource.builtinFile.fields.path}" - #end if - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -**What it does** - -This methylation caller parses the BSMAP SAM output file into bed format. - - -**Output format** :: - - - Column Description - ---------------------- -------------------------------------- - 1 chr chromosome - 2 pos position - 3 strand strand - 4 context context (CHH,CHG,CpG) - 5 coverage totally sequenced Cs at that position - 6 methylated methylated Cs at that position - 7 percentage methylated percentage of 6 - - -