<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; color: rgb(0, 0, 0); font-size: 14px; font-family: Calibri, sans-serif; ">
<div>
<div>
<div>
<blockquote type="cite" style="font-family: Calibri; font-size: medium; ">
<div>2. Is there any way to track the maker gene to its corresponding query proteins through setting up some parameters in maker? Or is it possible to directly let maker name the gene model based on its query protein?</div>
</blockquote>
</div>
</div>
</div>
<div>
<div>Not directly. Proteins can align to multiple locations, so you will get more than one protein overlapping a gene. A reciprocal BLAST would probably be the best thing to do here. If the MAKER derived annotation, and a gene are both each others best hit
in separate blast jobs, then you could usually safely call them the same gene. You would have to do that outside of MAKER.</div>
</div>
<div><br>
</div>
<div>Thanks,</div>
<div>Carson</div>
<div><br>
</div>
<div><br>
</div>
<span id="OLK_SRC_BODY_SECTION">
<div style="font-family:Calibri; font-size:11pt; text-align:left; color:black; BORDER-BOTTOM: medium none; BORDER-LEFT: medium none; PADDING-BOTTOM: 0in; PADDING-LEFT: 0in; PADDING-RIGHT: 0in; BORDER-TOP: #b5c4df 1pt solid; BORDER-RIGHT: medium none; PADDING-TOP: 3pt">
<span style="font-weight:bold">From: </span>Barry Moore <<a href="mailto:barry.moore@genetics.utah.edu">barry.moore@genetics.utah.edu</a>><br>
<span style="font-weight:bold">Date: </span>Thursday, 6 September, 2012 12:35 PM<br>
<span style="font-weight:bold">To: </span>Zhou Qi <<a href="mailto:muntjaczhou@gmail.com">muntjaczhou@gmail.com</a>><br>
<span style="font-weight:bold">Cc: </span><<a href="mailto:maker-devel@yandell-lab.org">maker-devel@yandell-lab.org</a>><br>
<span style="font-weight:bold">Subject: </span>Re: [maker-devel] gene duplicates and track the query protein<br>
</div>
<div><br>
</div>
<div>
<div style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">
<br>
<div>
<div>On Sep 5, 2012, at 8:26 PM, Zhou Qi wrote:</div>
<br class="Apple-interchange-newline">
<blockquote type="cite">
<div>Dear maker users and developers:<br>
<br>
I'm new to maker and trying to use it to annotate a Drosophila genome given protein sequences of its closely related species. I have two questions:<br>
<br>
1. How maker decide which one to keep if one query protein has two similarly best hits in the target genome, e.g., gene duplication?<br>
</div>
</blockquote>
<div><br>
</div>
<div>MAKER uses the aligned proteins as evidence to support gene predictions (i.e. from SNAP). If the supporting protein can align within the specified parameters (as defined in maker_bopts.ctl) then all of the valid alignments could be used as evidence for
gene predictions at those loci.</div>
<br>
<blockquote type="cite">
<div><br>
2. Is there any way to track the maker gene to its corresponding query proteins through setting up some parameters in maker? Or is it possible to directly let maker name the gene model based on its query protein?<br>
<br>
</div>
</blockquote>
<div><br>
</div>
<div>I don't know of a way to do this by setting parameters in the MAKER control file - maybe Carson has some idea for you here?</div>
<div><br>
</div>
<div><br>
</div>
B</div>
<div><br>
<blockquote type="cite">
<div>I looked through the previous threads, one possible way seems to use model_gff and map_forward option combined. I guess people more experienced than me using maker might know other ways?<br>
<br>
Many thanks!<br>
<br>
Qi<br>
_______________________________________________<br>
maker-devel mailing list<br>
<a href="mailto:maker-devel@box290.bluehost.com">maker-devel@box290.bluehost.com</a><br>
<a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a><br>
</div>
</blockquote>
</div>
<br>
<div><span class="Apple-style-span" style="border-collapse: separate; color: rgb(0, 0, 0); font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-align: auto; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; ">
<div><span class="Apple-style-span" style="font-family: Arial; font-size: 12px; ">
<div>Barry Moore</div>
<div>Research Scientist</div>
<div>Dept. of Human Genetics</div>
<div>University of Utah</div>
<div>Salt Lake City, UT 84112</div>
<div>--------------------------------------------</div>
<div>(801) 585-3543</div>
<div><br class="khtml-block-placeholder">
</div>
</span></div>
<div><br>
</div>
</span><br class="Apple-interchange-newline">
</div>
<br>
</div>
</div>
_______________________________________________ maker-devel mailing list <a href="mailto:maker-devel@box290.bluehost.com">
maker-devel@box290.bluehost.com</a> <a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org">
http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a> </span>
</body>
</html>