txUnique: Count Number of Bases Unique to Each Transcript
In Rsubread: Mapping, quantification and variant analysis of sequencing data

Description Usage Arguments Details Value Author(s)

For each transcript, number of unique bases and total number of bases will be reported.

txUnique(

    # basic input/output options
    GTF_Annotation_File,
    Feature_Type = "exon",
    Gene_ID_Attribute = "gene_id",
    Transcript_ID_Attribute = "transcript_id")

`GTF_Annotation_File`	a character string giving the name of a GTF file containing the transcript-level annotations.
`Feature_Type`	a character string specifying the type of annotations. Only the annotations of the given type will be loaded from the GTF file. `exon` by default. The third column in the GTF file contains the type of the annotation in each line
`Gene_ID_Attribute`	a character string specifying the attribute name in the 9-th column in the GTF file for the gene identifier. `gene_id` by default.
`Transcript_ID_Attribute`	a character string specifying the attribute name in the 9-th column in the GTF file for the transcript identifier. `transcript_id` by default.

This function compares the transcripts belonging to the same gene and then counts the number of bases unique to each transcript. It also reports the total number of bases for each transcript. When a transcript is found to contain overlapping exons, the overlapping exons will be merged and each overlapping base will be counted only once.

A data matrix in which each row represents a transcript. The matrix includes four columns: gene identifier, transcript identifier, number of bases unique to each transcript and total number of bases each transcript has.

Yang Liao and Wei Shi

Rsubread documentation built on March 17, 2021, 6:01 p.m.