Hi everyone,
I am looking for a data analysis pipeline specifically for RNA-Seq data obtained from 3' mRNA-Seq assays (i.e. only the 3' end of the fragmented transcript is sequenced). Particularly, I would like to know which STAR (for the alignment step) and HTSeq (for the quantification step) parameters should be used and to what values should they be set in this case?
Thanks a lot for your help!
•
written
by
Ld_60 •
1 answer
Hi Ld_60,
I guess it depends on the library prep that was used. If it targets the actual TTS, you'll might need rather an end-to-end alignment rather than splice detection.
If your reads are not directly on the TTS and rather randomly distributed, you'll go very well with the 'ENCODE Options' according to the manual.
I'd check the alignments with RSEQC'S read_distribution, infer_experinent, and geneBodyCoverage analysis tools. There you'll see if your annotation is "too short", which strandedness your libraries have (for HTSEQ), and how your reads are distributed over the genes.
Cheers
Michael
Log in to answer this question.
-
Current best practice for alignment and quantification for downstream analysis in DESeq2
written by yura.grabovska •
-
Downstream analysis from STAR Alignment at transcript level
written by bassanio •
-
How to get fpkm of TCGA data
written by Maryam •
-
STAR ALIGNMENT PARAMETERS FOR RNA-SEQ OF THE PLANT SAMPLE ?
written by Meeran •
-
Why does TCGA RNA-seq pipeline starts from a BAM file?
written by シン •
-
STAR for mouse RNA-seq alignment, which parameters should I use besides the basic ones in practice?
written by FantasticAI •
-
Normalisation for single nuclei RNA-seq
written by MutationalMeltdown
-
STAR aligner - multimaps vs. strandedness
written by cjag •
-
Recommendations for miRNA-Seq quantification with Salmon
written by labadorf •
-
gdc RNA-Seq Pipeline gene expression Quantification parameters
written by young yu •
