tests/testthat/test-getInsSeqs.R

In HLindsay/AlignmentTools: Common operations on GenomicAlignments

context("Get insertion sequences from aligned sequences")


test_that("", {


    # wdths <- c(10,10,5,3)
    #gr <- GenomicRanges::GRanges("chr1",
    #                           IRanges::IRanges(start = c(5,7,4,13), width = wdths),
    #                           cigar = sprintf("%sM", wdths), strand = "+",
    #                           flag=c(0,0,0,2048))
    #names(gr) <- c("A", "B", "C","C")
    #seqs <- Biostrings::DNAStringSet(subseq(rep("ACTGACTGAC",
    #                                          length(gr)),1, width(gr)))
    #gr$seq <- seqs
    #galns <- as(gr, "GAlignments"))


})