Question: StringTie missing header
gravatar for SarRo
12 months ago by
SarRo0 wrote:


I am a newbie to RNA seq and differential gene expression and I'm struggling with stringtie.

After merging all my stringtie outputs into a reference .gtf I've got same error messages for every sample.

samopen] no @SQ lines in the header.
[sam_read1] missing header? Abort!

I've used the following command:

stringtie -e -B -p 4 -G /merged_gtf_location/stringie_merged.gtf -o stringtie_out.gtf sample1.gtf --rf

When I converted my hisat2 outputs into .bam files I checked for headers and also running stringtie the first time didn't give me that error. I guess it's just a very stupid mistake but I really can't figure it out on my own and I'd be really happy for any help.

rna-seq assembly • 496 views
ADD COMMENTlink modified 12 months ago by lakhujanivijay4.8k • written 12 months ago by SarRo0

Apparently I didn't. Wow, that's embarrassing! I've checked my command so many times and it really didn't occur to me... Thank you!

ADD REPLYlink written 12 months ago by SarRo0

If an answer was helpful, you should upvote it; if the answer resolved your question, you should mark it as accepted. You can accept more than one if they work.

ADD REPLYlink written 12 months ago by genomax80k
gravatar for swbarnes2
12 months ago by
United States
swbarnes27.5k wrote:

Umm, where in that command line are you telling StringTie where your bam is?

ADD COMMENTlink written 12 months ago by swbarnes27.5k
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: 1521 users visited in the last hour