Question: Stringtie error with gtf file
0
gravatar for williamsbrian5064
8 months ago by
williamsbrian5064100 wrote:

Hi! I am getting this error when I try to run stringtie on some RNA-seq data. This is what I am getting.

WARNING: no reference transcripts were found for the genomic sequences where reads were mapped!

Please make sure the -G annotation file uses the same naming convention for the genome sequences.

I am using a gtf file from flybase. I am wondering if it is just not the correct file? I am working on a project with Drosophila. Does anyone know anything about this or could help?

The index that I used for hisat2 was just the premade Drosophilia index they have on their site. Could it be an issue of the index is a older version than the gtf file?

Thanks in advance

ADD COMMENTlink written 8 months ago by williamsbrian5064100

It seems that the header of reference fasta file and naming used in gtf/gff (1st column ) is not matching. Have you downloaded reference fasta & gtf/gff file from same DB and is it of similar version?

ADD REPLYlink written 8 months ago by MSM5570

I believe a similar version. Do you think that is the issue?

ADD REPLYlink written 8 months ago by williamsbrian5064100

Maybe I need to build my own index?

ADD REPLYlink written 8 months ago by williamsbrian5064100

Where did you get the indexes from? Can you post the output of following commands (use real file names that you have)

$ grep "^>" genome.fa 
$ head -10 your_file.gtf
ADD REPLYlink written 8 months ago by genomax55k
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: 618 users visited in the last hour