X-Git-Url: https://git.donarmstrong.com/?a=blobdiff_plain;f=bp_bin%2Ffind_orphans;h=20cc9c27a5dc2ff74e193c1e18551567547bd886;hb=b1fadc0ce9224b6a13f8deae332bb902a54c65af;hp=d8f2e4e6d8f80784618b4af1d8761c67d24447ae;hpb=cb41f472b532613a00a69d2cb03afc88d4f4811d;p=biopieces.git

diff --git a/bp_bin/find_orphans b/bp_bin/find_orphans
index d8f2e4e..20cc9c2 100755
--- a/bp_bin/find_orphans
+++ b/bp_bin/find_orphans
@@ -24,7 +24,7 @@
 
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> DESCRIPTION <<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 
-# Identify paired end orphan in records with sequence data.
+# Find paired end orphan in records with sequence data.
 
 # >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>><<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<
 
@@ -48,22 +48,24 @@ records = []
 
 Biopieces.open(options[:stream_in], options[:stream_out]) do |input, output|
   input.each_record do |record|
-    records << record
+    if record[:SEQ_NAME]
+      records << record
 
-    if records.size == 2
-      if Illumina.name_match?(records.first[:SEQ_NAME], records.last[:SEQ_NAME])
-        records.first[:TYPE] = "paired"
-        records.last[:TYPE]  = "paired"
+      if records.size == 2
+        if Illumina.name_match?(records.first[:SEQ_NAME], records.last[:SEQ_NAME])
+          records.first[:TYPE] = "paired"
+          records.last[:TYPE]  = "paired"
 
-        output.puts records.first
-        output.puts records.last
+          output.puts records.first
+          output.puts records.last
 
-        records = []
-      else
-        records.first[:TYPE] = "orphan"
-        output.puts records.first
+          records = []
+        else
+          records.first[:TYPE] = "orphan"
+          output.puts records.first
 
-        records.shift
+          records.shift
+        end
       end
     end
   end