[maker-devel] mapping cDNA to updated genome

[Prashant S Hosmani]

Hi All,

I am working on updating a plant genome annotation. I would like to map genes from previous annotation to a new genome build. There is a protocol about this in Campbell et al 2014, current protocols in bioinformatics (basic protocol 4 - Mapping annotations to a new assembly). I followed that protocol exactly with setting est_forward=1. But in output I’m getting large number of genes. My input cDNA fasta contains ~35K genes and after mapping there are ~58K genes.

I’m using maker version 3.0. There are few changes in the genome and I’m not expecting many changes in the mapping previous genes.

Please let me know if there are any other parameters to control mapping of EST’s. I was hoping to get similar number of genes mapped on to new assembly with very few changes.

Thank you for your help in advance.
Prashant


Prashant Hosmani
Sol Genomics Network
Boyce Thompson Institute, Ithaca, NY, USA



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://yandell-lab.org/pipermail/maker-devel_yandell-lab.org/attachments/20160831/1190998f/attachment-0002.html>