]> git.donarmstrong.com Git - biopieces.git/commitdiff
fixed paths in assemble_seq_idba + velvet
authormartinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Wed, 18 May 2011 15:18:59 +0000 (15:18 +0000)
committermartinahansen <martinahansen@74ccb610-7750-0410-82ae-013aeee3265d>
Wed, 18 May 2011 15:18:59 +0000 (15:18 +0000)
git-svn-id: http://biopieces.googlecode.com/svn/trunk@1409 74ccb610-7750-0410-82ae-013aeee3265d

bp_bin/assemble_seq_idba
bp_bin/assemble_seq_velvet

index da2f4dab4d4a0b03e56bce0ee8669954f86eb24b..ef35973f6ea627f1fa280b01d182d5437e2cbd60 100755 (executable)
@@ -48,7 +48,7 @@ options = bp.parse(ARGV, casts)
 
 Dir.mkdir(options[:directory]) unless Dir.exists?(options[:directory])
 
-file_fasta = [options[:directory], "IDBA"].join(File::SEPARATOR) + ".fna"
+file_fasta = File.join(options[:directory], "IDBA") + ".fna"
 
 Fasta.open(file_fasta, mode="w") do |fasta_io|
   bp.each_record do |record|
@@ -57,7 +57,7 @@ Fasta.open(file_fasta, mode="w") do |fasta_io|
 end
 
 unless File.size(file_fasta) == 0
-  output = [options[:directory], "IDBA"].join(File::SEPARATOR)
+  output = File.join(options[:directory], "IDBA")
 
   commands = []
   commands << "nice -n 19"
@@ -82,7 +82,7 @@ unless File.size(file_fasta) == 0
     $stderr.puts "Command failed: #{command}"
   end
 
-  file_contig = [options[:directory], "IDBA"].join(File::SEPARATOR) + "-contig.fa"
+  file_contig = File.join(options[:directory], "IDBA") + "-contig.fa"
 
   Fasta.open(file_contig, mode="r") do |fasta_io|
     fasta_io.each do |entry|
index c3e25f35e35119b582ceaef72baa8bdee1142280..08629852c921b42ed75011879489e6a771520a9e 100755 (executable)
@@ -45,7 +45,7 @@ class Velvet
     kmer = @kmer_min
 
     while kmer <= @kmer_max
-      dir_velveth = [@directory, "Kmer_#{kmer}"].join(File::SEPARATOR)
+      dir_velveth = File.join(@directory, "Kmer_#{kmer}")
 
       Dir.mkdir(dir_velveth)
 
@@ -68,7 +68,7 @@ class Velvet
       files_velveth = Dir.glob("#{dir_velveth}/*")
 
       cov_cutoffs.each do |cov_cutoff|
-        dir_velvetg = [dir_velveth, "Cov_cutoff_#{cov_cutoff}"].join(File::SEPARATOR)
+        dir_velvetg = File.join(dir_velveth, "Cov_cutoff_#{cov_cutoff}")
 
         Dir.mkdir(dir_velvetg)
         FileUtils.cp_r files_velveth, dir_velvetg
@@ -160,7 +160,7 @@ raise ArgumentError, "kmer_min >= kmer_max: #{options[:kmer_min]} >= #{options[:
 
 Dir.mkdir(options[:directory]) unless Dir.exists?(options[:directory])
 
-file_fasta = [options[:directory], "sequence_in.fna"].join(File::SEPARATOR)
+file_fasta = File.join(options[:directory], "sequence_in.fna")
 
 Fasta.open(file_fasta, mode="w") do |fasta_io|
   bp.each_record do |record|