<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Indeed, this stems from incorrect reading of the fasta file of the Genbank record (GBK). "lseq" has value e.g. 2 somehow, when it should be ~ 4.5 Million.<br>
</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> seqinr-forum-bounces@lists.r-forge.r-project.org <seqinr-forum-bounces@lists.r-forge.r-project.org> on behalf of Matthew Mark Parks <map2085@med.cornell.edu><br>
<b>Sent:</b> Monday, July 24, 2017 8:01:40 PM<br>
<b>To:</b> seqinr-forum@lists.r-forge.r-project.org<br>
<b>Subject:</b> [Seqinr-forum] Error in oriloc for many sequences</font>
<div> </div>
</div>
<div>
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p>Dear SeqinR developers,</p>
<p><br>
</p>
<p>The function "oriloc" returns the following error:</p>
<p><br>
</p>
<blockquote>
<p>ERROR : wrong sign in 'by' argument<br>
</p>
</blockquote>
<p><br>
</p>
<p>For many records from Genbank, even for very common, well documented organisms.</p>
<p><br>
</p>
<p>For example:</p>
<blockquote>
<p>library(seqinr)</p>
<p>oriloc(gbk = "<span>U00096.gb" )</span><br>
</p>
<p><span>ERROR : wrong sign in 'by' argument</span></p>
</blockquote>
<p>This occurs for many other Genbank records as well, including: <span>CP007592.1</span>,
<span>AMKD01000101.1, <span>JXMX01000001.1</span></span><br>
</p>
<p><br>
</p>
<p>Thank you,</p>
<p>Matthew<br>
</p>
<p><br>
</p>
<p><br>
</p>
</div>
</div>
</body>
</html>