###

test_cigarQNarrow <- function()
{
    cigar <- c("25M4D10M", "6S17M6I3M3S")

    ans <- cigarQNarrow(cigar)
    ans0 <- cigar
    attr(ans0, "rshift") <- c(0L, 0L)
    checkEquals(ans, ans0)

    ans <- cigarQNarrow(cigar, start=3, end=-3)
    ans0 <- c("23M4D8M", "4S17M6I3M1S")
    attr(ans0, "rshift") <- c(2L, 0L)
    checkEquals(ans, ans0)

    ans <- cigarQNarrow(cigar, start=7, end=-4)
    ans0 <- c("19M4D7M", "17M6I3M")
    attr(ans0, "rshift") <- c(6L, 0L)
    checkEquals(ans, ans0)

    ans <- cigarQNarrow(cigar, start=8, end=-5)
    ans0 <- c("18M4D6M", "16M6I2M")
    attr(ans0, "rshift") <- c(7L, 1L)
    checkEquals(ans, ans0)

    ans <- cigarQNarrow(cigar, start=25)
    ans0 <- c("1M4D10M", "5I3M3S")
    attr(ans0, "rshift") <- c(24L, 17L)
    checkEquals(ans, ans0)

    ans <- cigarQNarrow(cigar, start=26)
    ans0 <- c("10M", "4I3M3S")
    attr(ans0, "rshift") <- c(29L, 17L)
    checkEquals(ans, ans0)

    ans <- cigarQNarrow(cigar, start=26, end=-8)
    ans0 <- c("3M", "3I")
    attr(ans0, "rshift") <- c(29L, 17L)
    checkEquals(ans, ans0)

    ans <- cigarQNarrow(cigar, start=26, end=-10)
    ans0 <- c("1M", "1I")
    attr(ans0, "rshift") <- c(29L, 17L)
    checkEquals(ans, ans0)
}

test_cigarNarrow <- function()
{
    cigar <- c("25M4D10M", "6S17M6I3M3S")

    ans <- cigarNarrow(cigar)
    ans0 <- c("25M4D10M", "17M6I3M")
    attr(ans0, "rshift") <- c(0L, 0L)
    checkEquals(ans, ans0)

    ans <- cigarNarrow(cigar, start=3, end=-3)
    ans0 <- c("23M4D8M", "15M6I1M")
    attr(ans0, "rshift") <- c(2L, 2L)
    checkEquals(ans, ans0)

    ans <- cigarNarrow(cigar, start=7, end=-4)
    ans0 <- c("19M4D7M", "11M")
    attr(ans0, "rshift") <- c(6L, 6L)
    checkEquals(ans, ans0)

    ans <- cigarNarrow(cigar, start=8, end=-5)
    ans0 <- c("18M4D6M", "9M")
    attr(ans0, "rshift") <- c(7L, 7L)
    checkEquals(ans, ans0)

    ans <- cigarNarrow(cigar[1], start=26, end=-10)
    ans0 <- "1M"
    attr(ans0, "rshift") <- 29L
    checkEquals(ans, ans0)
}

test_cigarToIRangesListByRName <- function()
{
    cigar <- c("30M5000N10M", "50M", "90M10M5I50M10D40M", "18M10I22M", "99I")
    rname <- c("chr6", "chr6", "chr2", "chr6", "chr2")
    pos <- c(101, 201, 1001,  301, 2001)

    ans <- as.list(cigarToIRangesListByRName(cigar, rname, pos,
                                             merge.ranges=FALSE))
    ir2 <- IRanges(start=c(1001, 1091, 1101, 1151, 1161),
                   end=c(1090, 1100, 1150, 1160, 1200))
    ir6 <- IRanges(start=c(101, 5131, 201, 301, 319),
                   end=c(130, 5140, 250, 318, 340))
    ans0 <- list(chr2=ir2, chr6=ir6)
    checkEquals(ans, ans0)

    ans <- as.list(cigarToIRangesListByRName(cigar, rname, pos,
                                             drop.D.ranges=TRUE,
                                             merge.ranges=FALSE))
    ans0 <- list(chr2=ir2[-4], chr6=ir6)
    checkEquals(ans, ans0)

    ans <- as.list(cigarToIRangesListByRName(cigar, rname, pos,
                                             drop.D.ranges=TRUE))
    ir6 <- IRanges(start=c(101, 5131, 201, 301),
                   end=c(130, 5140, 250, 340))
    ans0 <- list(chr2=reduce(ir2[-4]), chr6=ir6)
    checkEquals(ans, ans0)

    ans <- as.list(cigarToIRangesListByRName(cigar, rname, pos))
    ir6 <- IRanges(start=c(101, 5131, 201, 301),
                   end=c(130, 5140, 250, 340))
    ans0 <- list(chr2=reduce(ir2), chr6=ir6)
    checkEquals(ans, ans0)
}

