stringtie Error: input file cannot be found
2
0
Entering edit mode
5.6 years ago

Hello,

I want to use StringTie to assemble transcripts for each sample, but when I run the stringtie, I have an error:

Error: input file –l cannot be found!


here is my command:

stringtie -p 8 -G /annotation/Ghirsutum_458_v1.1.gene.gff3.gz -o NMHT-80-23.gtf –l NMHT-80-23 NMHT-80-23.bam


Please kindly let me know what should I do?

rna-seq RNA-Seq alignment sequence Assembly • 5.9k views
0
Entering edit mode

Have you tried giving a different name for your output transcripts for the -l option?

0
Entering edit mode

yes, I did but i still have that error. what format supposed to be the output file after -l ?

0
Entering edit mode

For POC, try eliminating the -l, and let it use the default STRG. You can do a sed later, sed -i 's/STRG/NMHT-80-23/g' NMHT-80-23.gtf if it works.

The user manual also as the bam file in front of the options:

stringtie <aligned_reads.bam> [options]*

5
Entering edit mode
5.4 years ago
Ellango R ▴ 50

This error may be due to strings error in the script options. I do suffer the same problem while i copied the script directly from article and paste it into terminal. (*Name of output perfix –l =error (which i copied from article); -l = working)

1
Entering edit mode

That's right. As Pertea pointed, in the paper code we can found '–' instead of '-'; the first one is a character not recognized by stringtie. The problem is fixed just using the small '-'.

0
Entering edit mode

That's exactly what happened to me. Thanks :D

0
Entering edit mode
4.2 years ago
spriyar ▴ 10

Have you tried using the full path for the reference files. I had the same problem and adding the path to the genome reference solved it! Good luck!