From 3c7ae0d1a6079f732c1603c32435497170d5d628 Mon Sep 17 00:00:00 2001
From: hanna <hanna@348d0f76-0448-11de-a6fe-93d51630548a>
Date: Tue, 1 Mar 2011 17:38:30 +0000
Subject: [PATCH] Special case handling of unmapped region in low memory
 sharder.

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@5346 348d0f76-0448-11de-a6fe-93d51630548a
---
 .../datasources/reads/LowMemoryIntervalSharder.java    | 10 ++++++++++
 1 file changed, 10 insertions(+)

diff --git a/java/src/org/broadinstitute/sting/gatk/datasources/reads/LowMemoryIntervalSharder.java b/java/src/org/broadinstitute/sting/gatk/datasources/reads/LowMemoryIntervalSharder.java
index 40c2cb1af..10d73216e 100644
--- a/java/src/org/broadinstitute/sting/gatk/datasources/reads/LowMemoryIntervalSharder.java
+++ b/java/src/org/broadinstitute/sting/gatk/datasources/reads/LowMemoryIntervalSharder.java
@@ -34,6 +34,7 @@ import org.broadinstitute.sting.utils.GenomeLocSortedSet;
 
 import java.util.ArrayList;
 import java.util.Arrays;
+import java.util.Collections;
 import java.util.Iterator;
 import java.util.List;
 import java.util.NoSuchElementException;
@@ -85,6 +86,15 @@ public class LowMemoryIntervalSharder implements Iterator<FilePointer> {
 
         nextFilePointer = null;
         while(nextFilePointer == null && currentLocus != null) {
+            // special case handling of the unmapped shard.
+            if(currentLocus == GenomeLoc.UNMAPPED) {
+                nextFilePointer = new FilePointer(GenomeLoc.UNMAPPED);
+                for(SAMReaderID id: dataSource.getReaderIDs())
+                    nextFilePointer.addFileSpans(id,null);
+                currentLocus = null;
+                continue;
+            }
+
             nextFilePointer = new FilePointer(currentLocus.getContig());
 
             int coveredRegionStart = 1;