This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Rna-Seq With Reference Genome But No Gene Annotations

I have an RNAseq data from a fungi and the public database for that particular fungi has whole genome shotgun sequences. But there are no gene annotations. Is it possible to analyze RNAseq data for differential expression by aligning to the shotgun reference genome.

Thanks

rnaseq seq

2 answers

Yes, for example you can use cufflinks to assemble trnascript without genome annotation: read the documentation here http://cufflinks.cbcb.umd.edu/

On idea is to use a well annotated related species (in your case another fungi ) and to map the gene sequences ( use cDNA sequences, so UTR+exons ) with gmap ( it can output a gff file) to have an basic annotation.

Log in to answer this question.