I have a set of metatranscriptomics reads. I would like to map them to a given genome. This genome is not an annotated genome (no GFF file available). It is an unannotated genome assembly. With this kind of data at hand, what is the best way (which tools) to find out if there are transcripts from the given genome available within the metatranscriptomics data?
1 answer
It is not that difficult to annotate a prokaryotic genome. Really, it can be done in under 10 minutes. A basic way of doing it would be to simply predict genes. Prodigal will do it in under a minute and create a .gff file. If you want to invest a bit more time, prokka will do the annotation and create many useful files, one of which is a .gff file.
Log in to answer this question.
-
Including unannotated transcripts in ballgown diffexpression analysis
written by gabrielpenun •
-
How to get raw read counts of genes/cds/transcripts from SAM/BAM file?
written by shail.nair05 •
-
Combining prodigal for metagenomics and transdecoder for metatranscriptomics in same dataset?
written by O.rka
-
Denovo RNA seq
written by kanwarjag
-
Generate GTF/GFF file (coordinates) from a FASTA annotated file.
written by marquezg48 •
-
Recent metatranscriptomics workflow
written by dark.lord •
-
De novo assembled transcript mapping to genome scaffold
written by Bioinfonext
-
How to make an annotation file (.gtf/.gff) for a customized reference genome?
written by Shahzad •
-
Rna Seq On Unannotated Genome
written by kanwarjag
-
Finding Novel Transcripts From A De Novo Transcriptome Assembly
written by ilovepython •
