From 86acfd47f12db220e9aa2809357c6feb9b81738d Mon Sep 17 00:00:00 2001 From: gaoyubang <1489582340@qq.com> Date: Sat, 24 Dec 2022 19:49:45 +0800 Subject: [PATCH] Update README.md description --- README.md | 2 ++ 1 file changed, 2 insertions(+) diff --git a/README.md b/README.md index 754cbba..97f6009 100644 --- a/README.md +++ b/README.md @@ -222,6 +222,8 @@ EHMT1 NM_001354612.2 NM_001354611.2 NM_001145527.2 NM_001354259.2 NM_001354263.2 ``` (4) the --support parameter: The minimum number of DRS reads supporting a modified m6A site in genomic coordinates from one million DRS reads. The default is 10. Due to the low sequencing depth for DRS reads, quantification of m6A modification in low abundance gene is difficult. With this option, the pipeline will attempt to normalize library using this formula: Total number of DRS reads/1,000, 000 to generate \'per million scaling factor'. Then the 'per million scaling factor' multiply reads from -r option to generate the cuttoff for the number of modified transcripts for each modified m6A site. For example, the option (-r = 10, total DRS reads=2, 000, 000) will generate (2000000/1000000)*10=20 as cuttoff. Than means that modified A base supported by at least 20 modified transcripts will be identified as modified m6A sites in genomic coordinates. +In the test data, the parameter was set 30442 dut to small library size。But in the actual data, the parameter was set 10 in default。 + The main output is the ratio.x.tsv in the output dir. The header of ratio.x.tsv.