# HG changeset patch
# User bgruening
# Date 1408005972 14400
# Node ID 30caca800c9b17e17c7e2445ffd9fac0e29af33c
# Parent 862fb59a9a256f5139c5cf3c115c8293e881d012
Uploaded
diff -r 862fb59a9a25 -r 30caca800c9b bismark_bowtie2_wrapper.xml
--- a/bismark_bowtie2_wrapper.xml Mon Apr 14 16:42:38 2014 -0400
+++ b/bismark_bowtie2_wrapper.xml Thu Aug 14 04:46:12 2014 -0400
@@ -7,14 +7,19 @@
samtools
bowtie2
-
+
+
+
+
+
+
bismark_wrapper.py
## Change this to accommodate the number of threads you have available.
--num-threads "\${GALAXY_SLOTS:-24}"
- --bismark_path \$SCRIPT_PATH
+ ##--bismark_path \$SCRIPT_PATH
--bowtie2
@@ -34,7 +39,6 @@
## Input parameters
##
-
#if $singlePaired.sPaired == "single":
--single-paired $singlePaired.input_singles
@@ -58,14 +62,17 @@
--mate1 #echo ','.join($mate1)
--mate2 #echo ','.join($mate2)
- #if $singlePaired.mate_list[0].input_mate1.ext == "fastqillumina":
- --phred64-quals
- --fastq
- #elif $singlePaired.mate_list[0].input_mate1.ext == "fastqsanger":
- --fastq
- #elif $singlePaired.mate_list[0].input_mate1.ext == "fasta":
- --fasta
- #end if
+ #for $mate_pair in $singlePaired.mate_list:
+ #if $mate_pair.input_mate1.ext == "fastqillumina":
+ --phred64-quals
+ --fastq
+ #elif $mate_pair.input_mate1.ext == "fastqsanger":
+ --fastq
+ #elif $mate_pair.input_mate1.ext == "fasta":
+ --fasta
+ #end if
+ #break
+ #end for
-I $singlePaired.minInsert
-X $singlePaired.maxInsert
@@ -275,18 +282,22 @@
-
+
- singlePaired['sPaired'] == "paired"
- params['settingsType'] == "custom"
- params['supressed_read_file'] is True
+
+ ((
+ singlePaired['sPaired'] == "paired" and
+ params['settingsType'] == "custom" and
+ params['suppressed_read_file'] is True
+ ))
+
@@ -295,9 +306,9 @@
-
+
@@ -319,18 +330,22 @@
-
+
- singlePaired['sPaired'] == "paired"
- params['settingsType'] == "custom"
- params['unmapped_read_file'] is True
+
+ ((
+ singlePaired['sPaired'] == "paired" and
+ params['settingsType'] == "custom" and
+ params['unmapped_read_file'] is True
+ ))
+
@@ -339,9 +354,9 @@
-
+
@@ -641,4 +656,7 @@
always used and its default value is set to 10.
+
+ 10.1093/bioinformatics/btr167
+
diff -r 862fb59a9a25 -r 30caca800c9b bismark_bowtie_wrapper.xml
--- a/bismark_bowtie_wrapper.xml Mon Apr 14 16:42:38 2014 -0400
+++ b/bismark_bowtie_wrapper.xml Thu Aug 14 04:46:12 2014 -0400
@@ -7,11 +7,16 @@
samtools
bowtie
-
+
+
+
+
+
+
bismark_wrapper.py
- --bismark_path \$SCRIPT_PATH
+ ##--bismark_path \$SCRIPT_PATH
##
## Bismark Genome Preparation, if desired.
@@ -53,14 +58,17 @@
--mate1 #echo ','.join($mate1)
--mate2 #echo ','.join($mate2)
- #if $singlePaired.mate_list[0].input_mate1.ext == "fastqillumina":
- --phred64-quals
- --fastq
- #elif $singlePaired.mate_list[0].input_mate1.ext == "fastqsanger":
- --fastq
- #elif $singlePaired.mate_list[0].input_mate1.ext == "fasta":
- --fasta
- #end if
+ #for $mate_pair in $singlePaired.mate_list:
+ #if $mate_pair.input_mate1.ext == "fastqillumina":
+ --phred64-quals
+ --fastq
+ #elif $mate_pair.input_mate1.ext == "fastqsanger":
+ --fastq
+ #elif $mate_pair.input_mate1.ext == "fasta":
+ --fasta
+ #end if
+ #break
+ #end for
-I $singlePaired.minInsert
-X $singlePaired.maxInsert
@@ -251,9 +259,9 @@
-
+
@@ -271,9 +279,9 @@
-
+
@@ -295,9 +303,9 @@
-
+
@@ -314,15 +322,13 @@
-
+
-
-
@@ -547,4 +553,7 @@
temporary folder can be either relative or absolute.
+
+ 10.1093/bioinformatics/btr167
+
diff -r 862fb59a9a25 -r 30caca800c9b bismark_methylation_extractor.xml
--- a/bismark_methylation_extractor.xml Mon Apr 14 16:42:38 2014 -0400
+++ b/bismark_methylation_extractor.xml Thu Aug 14 04:46:12 2014 -0400
@@ -13,7 +13,7 @@
--infile $input
- --bismark_path \$SCRIPT_PATH
+ #--bismark_path \$SCRIPT_PATH
#if $singlePaired.sPaired == "single":
--single-end
diff -r 862fb59a9a25 -r 30caca800c9b bismark_wrapper.py
--- a/bismark_wrapper.py Mon Apr 14 16:42:38 2014 -0400
+++ b/bismark_wrapper.py Thu Aug 14 04:46:12 2014 -0400
@@ -17,7 +17,6 @@
def __main__():
- print 'tempfile_location',tempfile.gettempdir()
#Parse Command Line
parser = argparse.ArgumentParser(description='Wrapper for the bismark bisulfite mapper.')
parser.add_argument( '-p', '--num-threads', dest='num_threads',
@@ -258,9 +257,6 @@
# Final bismark command:
cmd = cmd % arguments
- print 'bismark_cmd:', cmd
- #sys.stderr.write( cmd )
- #sys.exit(1)
# Run
try:
tmp_out = tempfile.NamedTemporaryFile().name
@@ -321,8 +317,8 @@
"""
#tmp_out = tempfile.NamedTemporaryFile( dir=output_dir ).name
tmp_stdout = open( tmp_out, 'wab' )
- tmp_err = tempfile.NamedTemporaryFile( dir=output_dir ).name
- tmp_stderr = open( tmp_err, 'wb' )
+ #tmp_err = tempfile.NamedTemporaryFile( dir=output_dir ).name
+ tmp_stderr = open( tmp_err, 'wab' )
tmp_res = tempfile.NamedTemporaryFile( dir= output_dir).name
@@ -338,18 +334,22 @@
if returncode != 0:
raise Exception, open( tmp_stderr.name ).read()
else:
- tmp_res = bam_files[0]
+ tmp_res = bam_files[0]
bam_path = "%s" % tmp_res
if os.path.exists( bam_path ):
- if args.sort_bam:
- cmd = 'samtools sort -@ %s %s %s' % (args.num_threads, bam_path, args.output)
- else:
- shutil.copy( bam_path, args.output )
+ if args.sort_bam:
+ cmd = 'samtools sort -@ %s %s sorted_bam' % (args.num_threads, bam_path)
+ proc = subprocess.Popen( args=shlex.split( cmd ) )
+ returncode = proc.wait()
+ if returncode != 0:
+ raise Exception("Error during '%s'" % cmd)
+ shutil.move( 'sorted_bam.bam', args.output )
+ else:
+ shutil.move( bam_path, args.output )
else:
stop_err( 'BAM file no found:\n' + str( bam_path ) )
-
# TODO: look for errors in program output.
diff -r 862fb59a9a25 -r 30caca800c9b tool_dependencies.xml
--- a/tool_dependencies.xml Mon Apr 14 16:42:38 2014 -0400
+++ b/tool_dependencies.xml Thu Aug 14 04:46:12 2014 -0400
@@ -1,7 +1,7 @@
-
+
$REPOSITORY_INSTALL_DIR