<html><body><div style="font-family: arial, helvetica, sans-serif; font-size: 12pt; color: #000000"><div><br></div><div>Hi !</div><div><br data-mce-bogus="1"></div><div>Thank you <span style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; background-color: #ffffff; display: inline !important; float: none;" data-mce-style="color: #000000; font-family: arial, helvetica, sans-serif; font-size: 16px; font-style: normal; font-variant-ligatures: normal; font-variant-caps: normal; font-weight: normal; letter-spacing: normal; orphans: 2; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-text-stroke-width: 0px; background-color: #ffffff; display: inline !important; float: none;">Prashant </span> for sharing your experience. Indeed using the same blast version 2.2.29 for makeblastdb seems to resolve the problem. It is looking to work fine for all the sequences except one as I have the message above:</div><div><br></div><div>Fatal exception (source file ../../easel/esl_sqio_ascii.c, line 2001):<br>Failed to fetch subsequence residues -- corrupt coords?<br>sh: line 1: 46520 Aborted (core dumped) /hmmer-3.1b2-linux-intel-x86_64/binaries/esl-sfetch -c 1242..19031 all-te.lib rnd-4_family-1731#DNA >> blastx_results-all-te.txt.fnolowm50seq<br></div><div><br data-mce-bogus="1"></div><div>Did you encounter this problem before?</div><div><br></div><div>Thank you for your help.</div><div><br data-mce-bogus="1"></div><div>Amine</div><div><br data-mce-bogus="1"></div><hr id="zwchr" data-marker="__DIVIDER__"><div data-marker="__HEADERS__"><b>De: </b>jiangn@msu.edu<br><b>À: </b>"Prashant S Hosmani" <psh65@cornell.edu><br><b>Cc: </b>"Michael Campbell" <michael.s.campbell1@gmail.com>, "Mohamed Amine Chebbi" <mohamed.amine.chebbi@univ-poitiers.fr><br><b>Envoyé: </b>Jeudi 3 Novembre 2016 23:54:05<br><b>Objet: </b>Re: [maker-devel] ProtExcluder1.2 Error<br></div><br><div data-marker="__QUOTED_TEXT__">
<title></title>
<p>
Hi Prashant,<br>
<br>
Thank you so much for sharing your experience. It is important to keep everything in the same version. I will remind users about this when we update it and I may need to bother you then.<br>
<br>
Best regards,<br>
<br>
Ning<br>
<br>
Quoting Prashant S Hosmani <psh65@cornell.edu>:<br>
<br>
> Hi Amine,<br>
><br>
> I was getting similar error. You need to be careful with the blast<br>
> versions. Try using the same blast version for makeblastdb. I was<br>
> using BLAST 2.2.29+. After recreating new blast database with same<br>
> version, it worked for me.<br>
><br>
> Hope this helps.<br>
> Prashant<br>
><br>
><br>
> Prashant Hosmani<br>
> Sol Genomics Network<br>
> Boyce Thompson Institute, Ithaca, NY, USA<br>
><br>
><br>
><br>
> On Nov 3, 2016, at 1:57 PM, Michael Campbell<br>
> <michael.s.campbell1@gmail.com<mailto:michael.s.campbell1@gmail.com>><br>
> wrote:<br>
><br>
> Hi Amine,<br>
><br>
> That script is maintained by Ning Jiang and Kevin Childs. They know<br>
> best what this script is expecting. I?ve ccd them on this email in<br>
> the hope that they can provide some direction.<br>
><br>
> Thanks,<br>
> Mike<br>
> On Nov 2, 2016, at 2:09 PM, Mohamed Amine Chebbi<br>
> <mohamed.amine.chebbi@univ-poitiers.fr<mailto:mohamed.amine.chebbi@univ-poitiers.fr>><br>
> wrote:<br>
><br>
> Hi!<br>
><br>
> I am working on creating a custom repeat library and I want to use<br>
> ProtExcluder1.2 to trim potential genes from my repeat sequences.<br>
> My blast version is BLAST 2.2.30+<br>
><br>
> I get this message error :<br>
><br>
> Can not open the seqfile test.lib_blast_results.txt.fnolowm50seq<br>
> mergeunmatchedregion.pl seqfile<br>
> Illegal division by zero at ProtExcluder1.2/GCcontent.pl line 122.<br>
><br>
> I wonder if you can help me to fix this.<br>
><br>
> Thank you.<br>
><br>
> Amine<br>
> _______________________________________________<br>
> maker-devel mailing list<br>
> maker-devel@box290.bluehost.com<mailto:maker-devel@box290.bluehost.com><br>
> <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><br>
><br>
> _______________________________________________<br>
> maker-devel mailing list<br>
> maker-devel@box290.bluehost.com<mailto:maker-devel@box290.bluehost.com><br>
> <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><br>
><br>
><br>
</p><br></div></div></body></html>