Question: Stringtie and prepDE
0
gravatar for panacotaforcota
12 months ago by
Russia
panacotaforcota0 wrote:

Hello!

I try to use Stringtie to generate read coverage tables for prepDE.py. So, I do

./stringtie-2.0.3/stringtie file_sorted.bam -B -e -G merged.gtf > file.gtf

Is it correct that this command return gtf file or not?

Then I do

./prepDE.py -i file.gtf -g ./gene.counts.csv -t ./transc.counts.csv

and get it

Error: Text file with sample ID and path invalid (19 StringTie transcript 3033990 3034589 1000 . . gene_id "MSTRG.1"; transcript_id "MSTRG.1.1"; cov "3.420559"; FPKM "69.257027"; TPM "153.416428";)

Why?

stringtie • 578 views
ADD COMMENTlink modified 12 months ago by ATpoint40k • written 12 months ago by panacotaforcota0
1

I don't know by heart what stringtie writes to the stdout. I'd rather use the output option for the GTF than writing all output to a file.

The error message shows a line without strand information (+ or -). I'd check where that got lost.

ADD REPLYlink written 12 months ago by michael.ante3.6k

Not helping with this problem but do you really need stringtie? If you work with well-annotated transcriptomes there is rarely a reason to assemble the transcriptome. THere are other tools to conveniently get a count matrix for DEG like featureCounts or quantifiers such as salmon or kallisto that you might want to check out.

Also, please use the code option to hightlight code and error messages. I did it for you this time. enter image description here

ADD REPLYlink modified 12 months ago • written 12 months ago by ATpoint40k
Please log in to add an answer.

Help
Access

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