Question: NA for gff compare class code
0
gravatar for zoegward
2.9 years ago by
zoegward100
zoegward100 wrote:

Hi,

I have used gffcompare to look at the transcripts in my data set and some of them are coming out with 'NA' for a class code. eg. the transcript is present when I create a transcript count matrix but when I pull out all of the class codes for my transcripts from gffcompare and merge this to the count matrix a few do not have any class code affiliated to them.

Could anyone explain how this would happen?

rna-seq • 1.3k views
ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by zoegward100

Please paste your code and any warning or error messages.

ADD REPLYlink written 2.9 years ago by Kevin Blighe68k

There aren't any errors or warning messages. My understanding is that you can use gffcompare to compare your transcripts with a reference annotation and each transcript is given a class code depending where it resides with respect to the reference annotation. I don't quite understand why some are getting no class code at all and wondered if anyone else had seen this and had a theory as to why it would happen. One explanation could be that the transcript doesn't satisfy any class but when looking at all of the class codes it looks as though every scenario is covered? class code descriptions

ADD REPLYlink written 2.9 years ago by zoegward100

Thanks for posting more information!

  • Were the GTFs generated in the same way?
  • If you're also suppyling a reference GTF / GFF, is it the same as you used for aligning the files?
  • Are you using gffcompare from the old TopHat2 / Cufflinks packages or the upgraded HiSat2 / StringTie packages? - looks like StringTie.

I would imagine that gffcompare would assign a 'NA', also, if the contig IDs did not match between your GTF and your counts file / aligned BAM, for whatever reason. This is usually down to 'chr1' versus '1', for example.

ADD REPLYlink modified 2.9 years ago • written 2.9 years ago by Kevin Blighe68k
0
gravatar for zoegward
2.9 years ago by
zoegward100
zoegward100 wrote:

Should anyone else have the same problem - it is a known bug in gffcompare.... gffcompare_issues

ADD COMMENTlink modified 2.9 years ago • written 2.9 years ago by zoegward100

Thanks for following up. Can you mark your answer as accepted just for others. Thanks also for providing the helpful link.

ADD REPLYlink written 2.9 years ago by Kevin Blighe68k
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: 1018 users visited in the last hour