zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

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
This commit is contained in:
hanna 2011-03-01 17:38:30 +00:00
parent dd30ad751a
commit 3c7ae0d1a6
1 changed files with 10 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

10
java/src/org/broadinstitute/sting/gatk/datasources/reads/LowMemoryIntervalSharder.java
View File

@ -34,6 +34,7 @@ import org.broadinstitute.sting.utils.GenomeLocSortedSet;
import java.util.ArrayList; import java.util.ArrayList;
import java.util.Arrays; import java.util.Arrays;
import java.util.Collections;
import java.util.Iterator; import java.util.Iterator;
import java.util.List; import java.util.List;
import java.util.NoSuchElementException; import java.util.NoSuchElementException;
@ -85,6 +86,15 @@ public class LowMemoryIntervalSharder implements Iterator<FilePointer> {
nextFilePointer = null; nextFilePointer = null;
while(nextFilePointer == null && currentLocus != 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()); nextFilePointer = new FilePointer(currentLocus.getContig());
int coveredRegionStart = 1; int coveredRegionStart = 1;