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

Bug fix 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@434 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
aaron 2009-04-15 19:12:18 +00:00
parent 1edfe48194
commit 67ea66c866
1 changed files with 6 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
java/src/org/broadinstitute/sting/gatk/dataSources/shards/LocusShardStrategy.java
View File

@ -96,8 +96,11 @@ public abstract class LocusShardStrategy implements ShardStrategy {
this.dic = dic; this.dic = dic;
this.intervals = intervals; this.intervals = intervals;
this.currentInterval = 0; this.currentInterval = 0;
// set the starting point to the beginning interval
mLoc = new GenomeLoc(0, 0, 0); if (intervals.size() < 1) {
throw new IllegalArgumentException("Interval files must contain at least one interval");
}
mLoc = new GenomeLoc(intervals.get(0).getContig(),intervals.get(0).getStart()-1,intervals.get(0).getStart()-1);
if (dic.getSequences().size() > 0) { if (dic.getSequences().size() > 0) {
nextContig = true; nextContig = true;
} }
@ -162,7 +165,7 @@ public abstract class LocusShardStrategy implements ShardStrategy {
lastGenomeLocSize = proposedSize; lastGenomeLocSize = proposedSize;
// the next sequence should start at the begining of the next contig // the next sequence should start at the begining of the next contig
Shard ret = LocusShard.toShard(new GenomeLoc(intervals.get(currentInterval).getContigIndex(), nextStart, nextStart + proposedSize - 1)); Shard ret = LocusShard.toShard(new GenomeLoc(intervals.get(currentInterval).getContigIndex(), nextStart, nextStart + proposedSize));
++currentInterval; ++currentInterval;
if (intervals.size() > currentInterval) { if (intervals.size() > currentInterval) {