Hi everyone,
I am trying to understand rmats pipeline for my project.iam trying to run the pipeline with test files provided in the tutorial. I am successful in running the pipeline using bam files .I have problem running pipeline using fastq files, I am confused about star index folder, where do I get star index folder from? The star index folder link mentioned tutorial has many files in it. I am new to NGS tutorial.
python rMATS-turbo-xxx-UCSx/rmats.py --s1 s1.txt --s2 s2.txt --gtf gtf/Homo_sapiens.Ensembl.GRCh37.72.gtf --bi ~/STARindex/hg19 --od out_test -t paired --nthread 6 --readLength 50 --tophatAnchor 8 --cstat 0.0001 --tstat 6
Any help would be appreciated
Thanks in advance
•
written
by
manasa22.reddy •
1 answer
As indicated in the documentation under the "Installation" tab at http://rnaseq-mats.sourceforge.net/user_guide.htm#install ,
Obtain STAR genome index for genome by either of the following two ways
- Download pre-built STAR indexes if using Human (hg38, hg19) or Mouse (mm10)
- Build your own STAR index following STAR manual from genome fasta sequence
You can download these indexes from http://labshare.cshl.edu/shares/gingeraslab/www-data/dobin/STAR/STARgenomes/
Log in to answer this question.
-
rMATS reports zero evenets?
written by PK •
-
rMats Run Does Not Generate More Than One Output Row Per File
written by Y •
-
How robust is rMATS when swapping samples?
written by Apex92 •
-
error in rmats when mapping
written by katie326 •
-
weird error while running rmats on fastq files using hg38 genome
written by Sara
-
Getting errors trying to run rmats
written by Omnia •
-
Detection of splicing defects in RNA-Seq data of Arabidopsis thaliana
written by julius •
-
rMATS 4.0.2 (turbo) error when running the statistical part
written by Ivó Hernández Hernández •
-
I got some diffculties when run rMATs
written by younglin113
-
input for rMATS
written by bisht20diksha •
I added code markup to your post for increased readability. You can do this by selecting the text and clicking the 101010 button. When you compose or edit a post that button is in your toolbar, see image below:
Furthermore, you posted this as a "Tool", a post type which is intended to announce new tools. I have converted it to a Question.
I up this topic because I have a question relative to this
I have done an analysis of alternative splicing using rMATS, with fastq files as input (so the pipeline execute by itself the mapping with STAR)
But recently I have heard about the 2 pass-mode (http://labshare.cshl.edu/shares/gingeraslab/www-data/dobin/STAR/STAR.posix/doc/STARmanual.pdf) that allow to detect novel splicing junctions
MY question is did people compare these two approaches? Is it worth it to restart my analysis with this 2-pass approach, and if so, how to find in the output file predicted splices de novo?
Thank's by advance