<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Mike and Carson,
<div class=""><br class="">
</div>
<div class="">Thank you for your help. I used masked genome for aligning cDNAs. And yes, this was due to multiple aligning cDNA’s. I guess you could also filter according genes based on the alignment score from gff. I used GMAP (<a href="http://research-pub.gene.com/gmap/" class="">http://research-pub.gene.com/gmap/</a>)
to align cDNA on to the updated genome. GMAP has parameters to filter based on alignment scores and also can choose best path per cDNA.</div>
<div class=""><br class="">
</div>
<div class="">Regards,</div>
<div class="">Prashant</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<div class="">Prashant Hosmani</div>
<div class="">
<div dir="ltr" style="font-family: arial; font-size: small;" class="">Sol Genomics Network<br class="">
</div>
<div dir="ltr" class="">
<div style="font-family: arial; font-size: small;" class="">Boyce Thompson Institute, Ithaca, NY, USA</div>
</div>
</div>
<div class=""><br class="">
</div>
<br class="Apple-interchange-newline">
</div>
<br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Aug 31, 2016, at 12:12 PM, Carson Holt <<a href="mailto:carsonhh@gmail.com" class="">carsonhh@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Also if you have multiple alignments of the same cDNA, you can use the score column of the mRNA feature to see which aligns best. If they have the same score, you will have to disambiguate manually or just remove all copies.
<div class=""><br class="">
</div>
<div class="">—Carson</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Aug 31, 2016, at 10:10 AM, Michael Campbell <<a href="mailto:michael.s.campbell1@gmail.com" class="">michael.s.campbell1@gmail.com</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Prashant,
<div class=""><br class="">
</div>
<div class="">I’m almost positive that the additional genes are coming from multiply aligning cDNAs. Did you repeat mask your genome before mapping things forward?</div>
<div class=""><br class="">
</div>
<div class="">Another thought, what kind of whole genome duplications has your plant been through. it may be that the multiple alignments are to pseudogenes is some stage of decay. If that is the case it would probably be safe to keep the the gene from longest/best
aligned cDNA.</div>
<div class=""><br class="">
</div>
<div class="">Thanks,</div>
<div class="">Mike<br class="">
<div class="">
<blockquote type="cite" class="">
<div class="">On Aug 31, 2016, at 10:35 AM, Prashant S Hosmani <<a href="mailto:psh65@cornell.edu" class="">psh65@cornell.edu</a>> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi All,
<div class=""><br class="">
</div>
<div class="">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. </div>
<div class=""><br class="">
</div>
<div class="">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.</div>
<div class=""><br class="">
</div>
<div class="">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.</div>
<div class=""><br class="">
</div>
<div class="">Thank you for your help in advance.</div>
<div class="">Prashant<br class="">
<div class=""><br class="">
</div>
<div class=""><br class="">
<div class="">
<div class="">Prashant Hosmani</div>
<div class="">
<div dir="ltr" style="font-family: arial; font-size: small;" class="">Sol Genomics Network<br class="">
</div>
<div dir="ltr" class="">
<div style="font-family: arial; font-size: small;" class="">Boyce Thompson Institute, Ithaca, NY, USA</div>
</div>
</div>
<div class=""><br class="">
</div>
<br class="Apple-interchange-newline">
</div>
<br class="">
</div>
</div>
</div>
_______________________________________________<br class="">
maker-devel mailing list<br class="">
<a href="mailto:maker-devel@box290.bluehost.com" class="">maker-devel@box290.bluehost.com</a><br class="">
<a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" class="">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a><br class="">
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
_______________________________________________<br class="">
maker-devel mailing list<br class="">
<a href="mailto:maker-devel@box290.bluehost.com" class="">maker-devel@box290.bluehost.com</a><br class="">
<a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" class="">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a><br class="">
</div>
</blockquote>
</div>
<br class="">
</div>
</div>
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>