Has anyone tried "star --readFilesManifest manifest.tsv" and then "featureCounts --byReadGroup" to quickly get gene quantification for multiple samples?
0
1
Entering edit mode
3.2 years ago
Xianjun ▴ 310

One typical workflow for RNA-seq is to run STAR on the individual sample and then call htseq-count or featureCount to get the read count per gene. Recently I noticed that STAR has an option --readFilesManifest manifest.tsv that you can map multiple samples in one run and generate one BAM file with reads group tag @RG to mark individual samples, then you can run featureCounts −−byReadGroup on the BAM file to separate the multiple samples into individual read count profiles. I am not sure if this is a faster and more efficient way to run on multiple samples (than looping thru individual files). Has anyone tried it? Love to hear feedback. Thanks.

RNA-Seq STAR featureCounts • 1.6k views
ADD COMMENT

Login before adding your answer.

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