From 5ab2ffba8c37dd9d620c1b069b289d3bd494f159 Mon Sep 17 00:00:00 2001 From: agouru55 <79283474+agouru55@users.noreply.github.com> Date: Mon, 5 Feb 2024 12:38:54 -0500 Subject: [PATCH] Update Usage.md --- doc/Usage.md | 21 +++++++++++++++++++++ 1 file changed, 21 insertions(+) diff --git a/doc/Usage.md b/doc/Usage.md index 452b5b3..822b6e9 100644 --- a/doc/Usage.md +++ b/doc/Usage.md @@ -9,6 +9,8 @@ liqa -task : quantify: quantify isoform expression diff: detect differential splicing gene/isoform + + novel: detect noval isoform ``` @@ -122,3 +124,22 @@ Output: • Column 3: Hellinger distance which measures the splicing difference between two conditions for a gene in terms of isoform relative abundances ``` + +## Novel isoform detection +For novel isoform detection, the user can use the following LIQA command and parameters. +``` +liqa -task novel + -refgene + -bam + -out + -num_cover <# bp coverage> + -num_support_read <# support reads> +``` +Here is an example of the command you can run: + +``` + liqa -task novel -refgene example.refgene -bam example.bam -out noveltest -num_cover 20 -num_support_read 10 +``` + + +The output is a GTF file with a list of genes/isoforms and their relative abundance.