# HG changeset patch # User yutaka-saito # Date 1429450766 -32400 # Node ID 20930a8f700bd62d011995052c4a941b3741750f # Parent 06f8460885ff41ca1065c21a03103ade72b40bf1 rename some files diff -r 06f8460885ff -r 20930a8f700b bsf-call_wrapper.pl --- a/bsf-call_wrapper.pl Sun Apr 19 20:51:13 2015 +0900 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,51 +0,0 @@ -#!/usr/bin/perl - -use strict; -use warnings; - -use FindBin; - -print STDOUT "The tool script is called with:\n", join(" ", ($0, @ARGV)), "\n\n"; - -my ($idx, $in) = ("", ""); -my $default_option = "-o bsf-call.out -W bsfwork"; - -my $tooldir = shift(@ARGV); -$tooldir = $FindBin::Bin; -$ENV{PATH} = "$tooldir/bin:" . $ENV{PATH}; -my $reference_source = shift(@ARGV); -my $read_end = shift(@ARGV); -my $gslot = shift(@ARGV); -#$idx = "$tooldir/data/chrX.sub.fa"; - -if ($reference_source eq "indexed") { - $idx = shift(@ARGV); -} -elsif ($reference_source eq "history") { - my $own = shift(@ARGV); - $idx = "reference.fa"; - &invoke_command("ln -s $own reference.fa"); -} -else { - die "never reach here\n"; -} - -if ($read_end eq "single-end") { - $in = shift(@ARGV); - &invoke_command("$tooldir/bin/bsf-call $default_option -p $gslot $idx $in"); -} -elsif ($read_end eq "paired-end") { - my $in1 = shift(@ARGV); - my $in2 = shift(@ARGV); - $in = $in1 . "," . $in2; - &invoke_command("$tooldir/bin/bsf-call $default_option -p $gslot $idx $in"); -} -else { - die "never reach here\n"; -} - -sub invoke_command { - my ($command) = @_; - print "invoking: $command\n"; - system($command); -} diff -r 06f8460885ff -r 20930a8f700b bsf-call_wrapper.xml --- a/bsf-call_wrapper.xml Sun Apr 19 20:51:13 2015 +0900 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,116 +0,0 @@ - - Mapping bisulfite-seq reads and calling methylated cytosines - - - - TOOLDIR - - - - bsf-call_wrapper.pl TOOLDIR $reference.source $read.end \${GALAXY_SLOTS:-1} - - #if $reference.source=="indexed": - $reference.index.fields.path - #else if $reference.source=="history": - $reference.own_file - #else - - #end if - - #if $read.end=="single-end" - $in - #else if $read.end=="paired-end" - $in1 $in2 - #else - - #end if - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -**bsf-call** - -Mapping bisulfite-seq reads and calling methylated cytosines - ------- - -**Input format** - -Inputs are bisulfite-seq reads in fastqsanger format (single-end or paired-end), and a reference genome index (built-in or constructed from your fasta file). - ------- - -**Output format** - -Output is a six-column tab-delimited file:: - - Col.| Description - ----+-------------------------------------- - 1 | chromosome label (e.g. chr1) - 2 | genomic position (0-based) - 3 | strand (+,-) - 4 | mC context (CG, CHG, CHH) - 5 | mC rate (float) - 6 | read coverage - ------- - -**Contact** - -Toutai Mituyama - -mituyama-toutai AT aist.go.jp - - - - 10.1093/nar/gkt1373 - - - - - diff -r 06f8460885ff -r 20930a8f700b bsfcall_wrapper.pl --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/bsfcall_wrapper.pl Sun Apr 19 22:39:26 2015 +0900 @@ -0,0 +1,51 @@ +#!/usr/bin/perl + +use strict; +use warnings; + +use FindBin; + +print STDOUT "The tool script is called with:\n", join(" ", ($0, @ARGV)), "\n\n"; + +my ($idx, $in) = ("", ""); +my $default_option = "-o bsf-call.out -W bsfwork"; + +my $tooldir = shift(@ARGV); +$tooldir = $FindBin::Bin; +$ENV{PATH} = "$tooldir/bin:" . $ENV{PATH}; +my $reference_source = shift(@ARGV); +my $read_end = shift(@ARGV); +my $gslot = shift(@ARGV); +#$idx = "$tooldir/data/chrX.sub.fa"; + +if ($reference_source eq "indexed") { + $idx = shift(@ARGV); +} +elsif ($reference_source eq "history") { + my $own = shift(@ARGV); + $idx = "reference.fa"; + &invoke_command("ln -s $own reference.fa"); +} +else { + die "never reach here\n"; +} + +if ($read_end eq "single-end") { + $in = shift(@ARGV); + &invoke_command("$tooldir/bin/bsf-call $default_option -p $gslot $idx $in"); +} +elsif ($read_end eq "paired-end") { + my $in1 = shift(@ARGV); + my $in2 = shift(@ARGV); + $in = $in1 . "," . $in2; + &invoke_command("$tooldir/bin/bsf-call $default_option -p $gslot $idx $in"); +} +else { + die "never reach here\n"; +} + +sub invoke_command { + my ($command) = @_; + print "invoking: $command\n"; + system($command); +} diff -r 06f8460885ff -r 20930a8f700b bsfcall_wrapper.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/bsfcall_wrapper.xml Sun Apr 19 22:39:26 2015 +0900 @@ -0,0 +1,116 @@ + + Mapping bisulfite-seq reads and calling methylated cytosines + + + + TOOLDIR + + + + bsfcall_wrapper.pl TOOLDIR $reference.source $read.end \${GALAXY_SLOTS:-1} + + #if $reference.source=="indexed": + $reference.index.fields.path + #else if $reference.source=="history": + $reference.own_file + #else + + #end if + + #if $read.end=="single-end" + $in + #else if $read.end=="paired-end" + $in1 $in2 + #else + + #end if + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +**bsf-call** + +Mapping bisulfite-seq reads and calling methylated cytosines + +------ + +**Input format** + +Inputs are bisulfite-seq reads in fastqsanger format (single-end or paired-end), and a reference genome index (built-in or constructed from your fasta file). + +------ + +**Output format** + +Output is a six-column tab-delimited file:: + + Col.| Description + ----+-------------------------------------- + 1 | chromosome label (e.g. chr1) + 2 | genomic position (0-based) + 3 | strand (+,-) + 4 | mC context (CG, CHG, CHH) + 5 | mC rate (float) + 6 | read coverage + +------ + +**Contact** + +Toutai Mituyama + +mituyama-toutai AT aist.go.jp + + + + 10.1093/nar/gkt1373 + + + + + diff -r 06f8460885ff -r 20930a8f700b tool-data/bsf-call_indexes.loc.sample --- a/tool-data/bsf-call_indexes.loc.sample Sun Apr 19 20:51:13 2015 +0900 +++ /dev/null Thu Jan 01 00:00:00 1970 +0000 @@ -1,2 +0,0 @@ -hg19 hg19 hg19 /disk/reference/Bisulfighter/Homo_sapiens/UCSC/hg19/Sequence/BsfcallIndex/genome.fa -test test test /disk/reference/Bisulfighter/Homo_sapiens/UCSC/hg19/Sequence/BsfcallIndex/chrX.sub.fa diff -r 06f8460885ff -r 20930a8f700b tool-data/bsfcall_indexes.loc.sample --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/tool-data/bsfcall_indexes.loc.sample Sun Apr 19 22:39:26 2015 +0900 @@ -0,0 +1,2 @@ +hg19 hg19 hg19 /disk/reference/Bisulfighter/Homo_sapiens/UCSC/hg19/Sequence/BsfcallIndex/genome.fa +test test test /disk/reference/Bisulfighter/Homo_sapiens/UCSC/hg19/Sequence/BsfcallIndex/chrX.sub.fa diff -r 06f8460885ff -r 20930a8f700b tool_data_table_conf.xml.sample --- a/tool_data_table_conf.xml.sample Sun Apr 19 20:51:13 2015 +0900 +++ b/tool_data_table_conf.xml.sample Sun Apr 19 22:39:26 2015 +0900 @@ -1,6 +1,6 @@ - +
value, dbkey, name, path - +