Calculate the GC content for a region

Calculate the GC content for a region

gc_content(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 dataframe with the sequnce reads, counts of GC bases, and proportion of GC per read

Examples

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