gravatar for Xianjun

2 hours ago by

Great Boston Area

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.


modified 1 hour ago

2 hours ago


Source link