diff --git a/java/src/org/broadinstitute/sting/utils/GenomeLocParser.java b/java/src/org/broadinstitute/sting/utils/GenomeLocParser.java index 7f22f7294..18a1e7ffd 100644 --- a/java/src/org/broadinstitute/sting/utils/GenomeLocParser.java +++ b/java/src/org/broadinstitute/sting/utils/GenomeLocParser.java @@ -435,7 +435,7 @@ public class GenomeLocParser { // read is unmapped and not placed anywhere on the genome return GenomeLoc.UNMAPPED; else { - // + // Use Math.max to ensure that end >= start (Picard assigns the end to reads that are entirely within an insertion as start-1) int end = read.getReadUnmappedFlag() ? read.getAlignmentStart() : Math.max(read.getAlignmentEnd(), read.getAlignmentStart()); return createGenomeLoc(read.getReferenceName(), read.getReferenceIndex(), read.getAlignmentStart(), end, false); }