diff --git a/java/src/org/broadinstitute/sting/playground/utils/GenomicMap.java b/java/src/org/broadinstitute/sting/playground/utils/GenomicMap.java index efbb067e1..46ef74be5 100644 --- a/java/src/org/broadinstitute/sting/playground/utils/GenomicMap.java +++ b/java/src/org/broadinstitute/sting/playground/utils/GenomicMap.java @@ -146,7 +146,13 @@ public class GenomicMap implements Iterable 0 && + segments.get(segments.size()-1).getStop()+1 == newSegment.getStart() && + segments.get(segments.size()-1).getContigIndex() == newSegment.getContigIndex()) + System.out.println("WARNING: strictly adjacent segments found in custom contig "+name); + + segments.add(newSegment); p1 = p2+1; // set p1 after the comma while ( p1 < line.length() && Character.isWhitespace(line.charAt(p1))) p1++; // skip whitespaces @@ -183,8 +189,8 @@ public class GenomicMap implements Iterable currStop + 1 ) { + // add "deletion" w/respect to the master ref over the region between adjacent segments + // (and do not add anything if segments are strictly adjacent, i.e. refPos == currStop+1): newCigar.add(new CigarElement((int)(refPos-currStop-1),CigarOperator.D)); } currStop = gl.getStop();