STAR mapping two pass mode for multiple samples with split steps or one step ?
0
0
Entering edit mode
7 months ago
octpus616 ▴ 100

Hi,

I was recently using STAR to compare RNAseq data, and I noticed that it has a two pass mode with different way:

STAR can run in a 2-pass mode in two ways,

1) By re-generating the genome indices from splice junctions obtained

2) from 1-pass By using splice junctions directly during the mapping step

In the first method, we should run the STAR and generated splice junctions (SJ.out.tab), you can use these splice junctions for genome index re-generation, and run STAR with re-generated genome index again. But in latter method, we only need a params --twopassMode Basic.

According to the manual, there seems to be some advantage to using the former, especially when I have a lot of samples, but I've noticed that a lot of the literature seems to use the latter setup, why is that? How should this parameter be selected?

Best Regrads

Zhang

STAR RNA-seq • 759 views
ADD COMMENT
0
Entering edit mode

Which is your experimental setting? Would you expect to see splicing alterations sample-specific or rather more batch/experiment specific?

ADD REPLY
0
Entering edit mode

I am runing a typical RNA-seq pipeline to quantify gene expresstion, in fact. But I'd like the got more precise alignment to support some custom analysis start from aligned reads. I think the two step 2-pass mode is more precise, but is seem like less using in literature

ADD REPLY

Login before adding your answer.

Traffic: 2334 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