This function calculates the position and length of the two sequenced segments from a single RNA template, potentially after merging when both segments overlap.

phsc.bam.get.length.and.pos.of.mergedreads(bam.file.name, error.strict = TRUE)

Arguments

bam.file.name

full path name to bam file.

Value

data.table with columns QNAME (template query ID), POS (leftmost position of read), LEN (length of read)