moved mum.rb

[biopieces.git] / bp_bin / bwa_seq

diff --git a/bp_bin/bwa_seq b/bp_bin/bwa_seq
index 43637ca5e9677e275c6e1ea037b74378f5cd122c..4fd89a253427fa7ab8f24c52785317cbf1fc1de3 100755 (executable)
--- a/bp_bin/bwa_seq
+++ b/bp_bin/bwa_seq
@@ -43,7 +43,6 @@ casts << {:long=>'cpus',       :short=>'c', :type=>'uint',   :mandatory=>false,
 options = Biopieces.options_parse(ARGV, casts)
 
 tmp_dir = Biopieces.mktmpdir
-tmp_dir = "Sletmig"
 tmp_fq  = File.join(tmp_dir, "bwa.fq")
 tmp_sai = File.join(tmp_dir, "bwa.sai")
 tmp_sam = File.join(tmp_dir, "bwa.sam")
@@ -117,7 +116,7 @@ Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
 
   Sam.open(tmp_sam, 'r') do |io_sam|
     io_sam.each do |entry|
-      output.puts Sam.to_bp(entry)
+      output.puts Sam.to_bp(entry) unless entry[:RNAME] == '*'
     end
   end
 end