Analysis of MiRDeep2 and DESeq2
0
0
Entering edit mode
4.5 years ago
guq4602 • 0

I want a pipeline for expression profiling and will proceed with DEanalysis later. I found the pipeline while i was searching for miRNA-seq data analysis.

using miRDeep2

mapper.pl input(fastq) -e -h -j -m -s output

quantifier.pl -p (Ref_hairpin from miRBase) -m (Ref_mature from miRBase) -r mapper.pl output) -t hsa -g 2 -e 2 -f 5

It was a pipeline that made a table by combining the ABC files obtained in this way and proceeded with DEanalysis using DESeq2.

miRNA           read_count  precursor       total   seq     seq(norm)
hsa-let-7a-5p   597882      hsa-let-7a-1    597882  597882  47424.19
hsa-let-7a-3p   29          hsa-let-7a-1    29      29      2.3

I think.. use the read_count value to run DESeq2......

I would like to get help because it is not mentioned why there is no mapping process in reference genome. And I'm wondering if it's right to analyze it like this pipeline.

RNA-Seq rna-seq R alignment sequencing • 1.5k views
ADD COMMENT
0
Entering edit mode

and I wonder if the reason for collapse is just a process to reduce computing power.

ADD REPLY

Login before adding your answer.

Traffic: 1984 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6