This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Stringtie not find the merged gtf reference

Hi there I am trying to use stringtie to count reads using the merged gtf from bam files but the command that I find it to do the job give me back an error.

$ stringtie –e –B -p 8 -G genome/TAIR10/annotation/stringtie_merged.gtf -o ballgown/Athaliana_root_1/Athaliana_root_1.gtf sam_bam/Athaliana_root_1.bam

Error: input file –e cannot be found!

I try to copy the merged file to different directories but not seems to work! I got this command from a tutorial. Do you have any clue or suggestion why this isn't working?
Thanks

rna-seq
• 2,371 views • link
written by schlogl schlogl

Stupid suggestion but have you looked in the manual/help section of the stringtie software? ;)

it seems like the first argument it expects is a file and not a parameter or other option. (== try to restructure you cmdline)

• 0 views • link
written by lieven.sterck lieven.sterck

indeed it is mentioned like this in this tutorial.

So either it's just wrong or it might be a version difference?

• 0 views • link
written by lieven.sterck lieven.sterck

2 answers

The syntax is wrong. From the manual: stringtie <aligned_reads.bam> [options]*

It means the input file is the first argument after stringtie which currently is -e, which is why the tool interpretes it as input.

Must be:

stringtie sam_bam/Athaliana_root_1.bam –e –B -p 8 -G genome/TAIR10/annotation/stringtie_merged.gtf -o ballgown/Athaliana_root_1/Athaliana_root_1.gtf
• 0 views • link
written by ATpoint ATpoint

Thanks ATpoint, generally I check the manual... but you right

stringtie <aligned_reads.bam> [options]*
• 0 views • link
written by schlogl schlogl

I got it done...

$ stringtie -e -B -p 4 sam_bam/Athaliana_root_1.bam -G genome/TAIR10/annotation/stringtie_merged.gtf -o ballgown/Athaliana_root_2/athaliana_root_1.count -A ballgown/Athaliana_root_1/Athaliana_root_1.out

Thanks

• 0 views • link
written by schlogl schlogl

Following ATpoint instructions this should/will not work , as the input bam is still not the first argument after the stringtie cmd.

• 0 views • link
written by lieven.sterck lieven.sterck

Log in to answer this question.

More posts like this