Bovine genome is already in ENSEMBL 61 and was annotated using the whole ENSEMBL annotation pipeline. Problem is, version of the genome used by ENSEMBL (lest call it A) is dated, and newer assembly (B) using exactly the same genomic DNA is available. While B is also annotated using Gnomon pipeline I want to be able to compare annotations and use ENSEMBL gene ids.
Since it will be months until B may be incorporated in ENSEMBL, I want to use B but with annotations without rerunning the whole annotation pipeline.
So my question is:
what is the best (or a t least sensible...) way to get gene annotations from one assembly mapped to another?
edit: changed assembly symbols (A=old, B=new) for clarity.
•
written
by
Darked89
1 answer
I would suggest you look at MAKER.
It has functionality specifically for migrating annotations from one assembly to another.
•
written
by
Scott Cain •
Log in to answer this question.
-
Which SPAdes output should be used for genome annotation: contigs.fasta or scaffolds.fasta?
written by Abhijit •
-
Protein secuence extraction from gpff file for a genes list
written by reza
-
Conversion between Ensembl identifiers and gene symbols obtained from different annotation files
written by lu.ne
-
How to quantify gene expression when stringtie considers two adjacent genes as single gene?
written by kousi31 •
-
Ensembl assembly choice for RNA-Seq reference mapping
written by eoin •
-
Best strategy for gene annotation for de novo genome assembly without RNA-seq data
written by soniabedi.07 •
-
Deprecated identifiers in Ensembl?
written by Solowars
-
some questions on GFF3
written by prasundutta87
-
Genome annotation using EST assembly
written by psiwach29 •
-
Genome Specific Database (Gbrowse / Ensembl Type)
written by Darked89
Do you mean "annotation" as in "finding proteins" or as in "assigning names to proteins"?
"annotation" as "describing chunks of DNA sequence as exons of a gene". Also "if gen X has following coordinates on genome assembly A, what are gene X coordinates on genome assembly B". Lastly, at this stage I am not interested in improving existing annotation for gene X. If it does miss an exon or two on assembly A, so be it, even if this region is present in assembly B.