get_splice_sites_info2: Get all splice-sites and read-end positions of a BAM files

In vladimirsouza/lrRNAseqBenchmark: Figure out how to call varaints from Iso-Seq data

Get all splice-sites and read-end positions of a BAM files

Description

Similar to get_splice_sites_info, but also count read ends.

Usage

get_splice_sites_info2(input_bam, threads)

Arguments

input_bam	The input BAM file to extract splice site positions from.
threads	Number of threads.

Details

This function may take several hours to run.

Value

A 2-element list ('splice_site', and 'transcript_end'). Element 'splice_site' should contain the same returned by get_splice_sites_info. Element 'transcript_end' counts of read ends (start or final/end site of the read).

vladimirsouza/lrRNAseqBenchmark documentation built on March 23, 2023, 7:32 a.m.