Extract the sequences for a given region

Extract the sequences for a given region

extract_sequence(bam, index, reg)

Arguments

bam	the cram/bam/sam file
index	the index of the cram/bam/sam file
reg	the region of interest, typically in format of chr1:start-begin

Value

a character vector with the sequences in the given region

Examples

# NOT RUN {
count_kmer(bam, index, "chr1:10001-100050")
# }