Stringtie --merge and gtf_list.txt issue
0
0
Entering edit mode
5.7 years ago
Morris_Chair ▴ 370

Hello guys,

I’m having problems merging the GTF files generated with Stringtie. I know it must me a simple mistake that I m doing with the cmdline but also googling I couldn’t figure it out.

The command line that I’m using for merging the files is this:

stringtie --merge -G Homo_sapiens.GRCh37.75.gtf -p 8 -o Treated_1.gtf Treated_2.gtf Treated_3.gtf Untreated_1.gtf Untreated_2.gtf Untreated_3.gtf gtf_list.txt

Error: no transcripts were found in input file gtf_list.txt

Beside the error that , it generates a folder with files inside named “Transcript_s000.bam and sam” but I can’t understand why it stops working.

In the file "gtf_list.txt' I listed the file name_gtf

Can you please help me with that?

Thank you

RNA-Seq • 7.3k views
ADD COMMENT
0
Entering edit mode

Paste the contents of gtf_list.txt

your contents should look like this

assembly/ERR188044_chrX.gtf
assembly/ERR188104_chrX.gtf
assembly/ERR188234_chrX.gtf
assembly/ERR188245_chrX.gtf
assembly/ERR188257_chrX.gtf
assembly/ERR188273_chrX.gtf
assembly/ERR188337_chrX.gtf
assembly/ERR188383_chrX.gtf
assembly/ERR188401_chrX.gtf
assembly/ERR188428_chrX.gtf
assembly/ERR188454_chrX.gtf
assembly/ERR204916_chrX.gtf
ADD REPLY
0
Entering edit mode

Hi bioExplorer,

here is the head of my gtf.txt file which is different than yours, did you put also the path beside of each? in my case the txt file is in the same folder of the gtf files so I think I should be ok

[@ws7910 Files.GTF]$ head gtf_list.txt
Treated_1.gtf
Treated_2.gtf
Treated_3.gtf
Untreated_1.gtf
Untreated_2.gtf
Untreated_3.gtf[@ws7910 Files.GTF]$

This gtf come out from this command line,

stringtie file.bam -o file.bam.To.gtf -G Homo_sapiens.GRCh37.75.gtf --rf -A Treated_1.txt

Thank you

ADD REPLY
1
Entering edit mode

Can you provide the absolute paths in the gtf_list.txt like this?

for e.g. say your gtf files are in /home/user/assembly directory, then your gtf_list.txt should look like

/home/user/assembly/ERR188428_chrX.gtf
/home/user/assembly/ERR188454_chrX.gtf
/home/user/assembly/ERR204916_chrX.gtf

AND, I noticed that your command is incorrect, it should be

stringtie --merge -G Homo_sapiens.GRCh37.75.gtf -p 8 -o stringtie_merged.gtf gtf_list.txt

stringtie --mergegenerated single output file called stringtie_merged.gtf.

ADD REPLY
0
Entering edit mode

Hole' ! problem solved , I have this precious file !

Basically I had to specify the path into the txt file

Thank you

PS: my previous command line was used to generate every single file.txt from single BAM files

ADD REPLY

Login before adding your answer.

Traffic: 1793 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6