This is a test version of Biostars. For the public version, visit https://www.biostars.org.
RNA-Seq DEGenes Without GTF

I am working a specific species where we have a reference genome, and I was able to align our RNA-seq data to it. However, I am looking to create a count table to funnel our dataset to edgeR for DE analysis. I don’t have a GTF to use with something like HTSeq. Anybody have an idea of how to proceed?

rna-seq

It's a species that has a reference genome but that is not well annotated in terms of gene expression? If that is the case, I would have done:

  1. HISAT2 / StringTie to produce a reference transcriptome GTF
  2. Use HTseq with the new GTF to derive raw counts
  3. Feed raw counts into EdgeR
  1. ...
  2. ...
  3. ...
  4. extract the unannotated gtf features and blast them to get an annotation.

0 answers

No answers yet.

Log in to answer this question.