Stranded setting in RSEM and RNASeQC
KangchengHou opened this issue · comments
Kangcheng Hou commented
According to the instruction at https://github.com/broadinstitute/gtex-pipeline/tree/master/rnaseq, when running RSEM, there will be a strand-specific setting. However, in the usage in the RNASeQC, we don't specify that. Could you explain why? Thanks very much!
François Aguet commented
These settings are specific to GTEx samples, which were sequenced with an unstranded protocol. Both tools have options for stranded protocols, and the latest wrappers (currently on the v9 branch) include them: https://github.com/broadinstitute/gtex-pipeline/blob/v9/rnaseq/src/run_rnaseqc.py