Mercurial > repos > cmonjeau > mapsembler2
diff mapsembler2.py @ 2:5aea5b993ae8
fix .dat unknown extension for kissreads
author | cmonjeau |
---|---|
date | Fri, 11 Sep 2015 10:04:40 +0000 |
parents | f905f982ee3d |
children | c62c5fd56ebf |
line wrap: on
line diff
--- a/mapsembler2.py Mon Jul 06 15:53:45 2015 +0000 +++ b/mapsembler2.py Fri Sep 11 10:04:40 2015 +0000 @@ -12,7 +12,7 @@ http://www.irisa.fr/symbiose/people/ppeterlongo/mapsembler2_2.2.3.zip -or with the galaxy_mapsembler2 package in the GenOuest toolshed +or with the package_mapsembler2_2_2_3 package in the GUGGO toolshed and main toolshed """ @@ -49,8 +49,13 @@ #inputs cmd_line.append("-r") - #cmd_line.append(inputs) - cmd_line.append(' '.join(options.input_files.split(","))) + inputs_tab = [] + + for input in options.input_files.split(","): + os.symlink(input, os.path.basename(input)+'.'+options.extension_format) + inputs_tab.append(os.path.basename(input)+'.'+options.extension_format) + + cmd_line.append(' '.join(inputs_tab)) # add parameters into the command line cmd_line.extend(["-t", options.output_extension, "-k", options.kmer, "-c", options.coverage, "-d", options.substitutions, "-g", options.genome_size, "-f", options.process_search, "-x", options.max_length, "-y", options.max_depth]) @@ -59,15 +64,8 @@ log = open(options.output, "w") log.write("[COMMAND LINE] "+' '.join(cmd_line)) - process=subprocess.Popen(cmd_line, - stdout=subprocess.PIPE, stderr=subprocess.PIPE) - - stdoutput, stderror = process.communicate() + process=subprocess.call(cmd_line) - # results recuperation - log.write(stdoutput) - log.write(stderror) - # close log file log.close()