Detect the most common type of alignment errors in full-length RNA-seq reads aligned by Minimap2
detect_alignment_errors(long_reads, min_read_support = 2)
GRangesList object which is identical to the input object, except for the extra logical column
GRangesList object is returned by the
The Minimap2 aligner is prone to under-split the full-length RNA-seq reads. In certain reads, two adjacent exons may appear as a single subalignment. Assuming that most of the reads are still aligned correctly, one can detect alignment errors by comparing each subalignment to the set of constitutive exons. The constitutive exons are defined here as exonic intervals supported by more than 50% of reads in given locus (but not less than
Subalignment with an alternative 5'- and/or 3-border (relative to a constitutive exon) are considered valid alternative exons, only if the next constitutive exon is also present in the read.
Otherwise, the alignment is marked as a possible alignment error.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.