Question: 'gene_name' is missing in StringTie output file 't_data.ctab'
gravatar for swarajgaikwad15
3 months ago by
Tezpur University, Assam, India
swarajgaikwad150 wrote:

Hello, I used New Tuxedo protocol in which StringTie is used for quantification step using:

stringtie -e -B -p 8 -G merged_gtf -o SRRXXX.gtf SRRXXX.bam

gives the output files as following:

e2t.ctab  e_data.ctab  i2t.ctab  i_data.ctab  SRRXXX.gtf  t_data.ctab

t_data.ctab columns are used for making countdata for DESeq2 using command

I tried to import t_data.ctab for DESeq2 with the help of tximport manual

   tx2gene <- tmp[, c("t_name", "gene_name")]

but my t_data.ctab contains '.' in 'gene_name' column, which is inappropriate for creation of countdata. Therefore I can't proceed my differential expression of genes. My question is Can I use 'gene_id' column instead of 'gene_name' from t_data.ctab. Or am I supposed to directly switch the quantification tool itself, if yes then which tool will be better as compared to StringTie?

ADD COMMENTlink modified 3 months ago by Viswanathan0 • written 3 months ago by swarajgaikwad150
gravatar for Viswanathan
3 months ago by
Viswanathan0 wrote:

To extract read count information, you can use the script provided by the StringTie authors.

ADD COMMENTlink written 3 months ago by Viswanathan0
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2090 users visited in the last hour