<div dir="ltr">Hi Carson,<div><br></div><div>Thank you very much for getting back to me on this. Sorry I was too quick to jump to a suspicion! Thanks for pointing it out and please see attached the EST (face_cDNA.fa) and genome scaffold sequences (albo_scaf3928+5347_for_face_cDNA.fa) I used for Maker. We used blastn to identify these two scaffolds for this cDNA sequence, and exonerate result also showed that this cDNA can be aligned to the scaffolds.</div>
<div><br></div><div>Thanks!</div><div><br></div><div>Xin </div></div><div class="gmail_extra"><br><br><div class="gmail_quote">On Mon, Nov 11, 2013 at 12:10 AM, Carson Holt <span dir="ltr"><<a href="mailto:carsonhh@gmail.com" target="_blank">carsonhh@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="font-size:14px;font-family:Calibri,sans-serif;word-wrap:break-word"><div>Sorry for the slow reply, I was moving from Canada to the US, got sick, and my moving truck got delayed one week. So I fell behind on the maker list.</div>
<div><br></div><div>What you observe is correct behavior. The config is blasted against the database of ESTs/proteins. It really doesn't matter which way you blast, you pick the direction that is most efficient for the question you are asking, and you only have to adjust the Z value for database size if you want alignment scores to be reproducible both ways.</div>
<div><br></div><div>Not getting results means they were filtered out for some reason (you don't even have blastn results in your output). If you can send a test dataset I can take a look and tell you which filter.</div>
<div><br></div><div>Thanks,</div><div>Carson</div><div><br></div><div><br></div><span><div style="border-right:medium none;padding-right:0in;padding-left:0in;padding-top:3pt;text-align:left;font-size:11pt;border-bottom:medium none;font-family:Calibri;border-top:#b5c4df 1pt solid;padding-bottom:0in;border-left:medium none">
<span style="font-weight:bold">From: </span> Xin Huang <<a href="mailto:xh33@georgetown.edu" target="_blank">xh33@georgetown.edu</a>><br><span style="font-weight:bold">Date: </span> Tuesday, November 5, 2013 2:40 PM<br>
<span style="font-weight:bold">To: </span> <<a href="mailto:maker-devel@yandell-lab.org" target="_blank">maker-devel@yandell-lab.org</a>><br><span style="font-weight:bold">Subject: </span> [maker-devel] A suspected strange behavior of Maker...<br>
</div><div><div class="h5"><div><br></div><div dir="ltr">Hi,<div><br></div><div><span style="font-family:arial,sans-serif;font-size:13px">I found a very bizarre thing when I ran Maker several times. The blastn research seems to reverse query and subject the way it's specified in the maker_opts.ctl file. When I put the scaffold sequence after genome= and the cDNA sequence after est=, the blastn result actually used the scaffold as query and the cDNA as subject. I realized that this sounds unreal, so I reversed the places for the two and the blastn result is normal with the cDNA as query and the scaffold as subject, but at the same time, the GFF file is of the cDNA sequence. I tried a few times, and got the same result. I deleted Maker (2.28) and reinstalled the beta version (2.30), still getting the same result.</span><br>
</div><div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">My wild guess was that this might be the reason why the est2genome prediction using the EST sequence didn't go into the GFF3 file, as stated in a previous email to the Maker-devel. I couldn't figure out how to test this, though. I wonder if there's anything dramatically wrong that I've been doing to get such odd results... If you'd like me to send you any output files or control files, please let me know and I'll promptly pass those along.</span></div>
<div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Thank you very much for considering my email.</span></div><div><span style="font-family:arial,sans-serif;font-size:13px"><br>
</span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Sincerely,</span></div><div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Xin Huang</span></div>
</div></div></div>
_______________________________________________
maker-devel mailing list
<a href="mailto:maker-devel@box290.bluehost.com" target="_blank">maker-devel@box290.bluehost.com</a>
<a href="http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org" target="_blank">http://box290.bluehost.com/mailman/listinfo/maker-devel_yandell-lab.org</a>
</span></div>
</blockquote></div><br></div>