<div dir="ltr"><div><div>Thank you very much for the quick responses and tips, I&#39;ll battle that guide out tomorrow, another dumb question I have is is there any easy way to make a graph out of the histograms khmer makes to see if I have neat separate peaks of khmer distributions? Here is the sort of thing I&#39;m looking for : <a href="http://ged.msu.edu/angus/_images/stamps-reads.png" target="_blank">http://ged.msu.edu/angus/_images/stamps-reads.png</a><br><br></div>Thanks in advance!<br></div>Diane<br><div><br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature"><div dir="ltr"><pre cols="72">-- 
Dr Diane Hatziioanou
Greek Mobile: <a value="+306909403373">(+30)6909403373</a><br>UK Mobile: <a value="+447779516625">(+44)7779516625</a>
<a href="http://www.linkedin.com/in/dhatziioanou/" target="_blank">www.linkedin.com/in/dhatziioanou/</a>
<a href="https://twitter.com/DianeHIoanou" target="_blank">https://twitter.com/DianeHIoanou</a></pre></div></div></div>
<br><div class="gmail_quote">On 3 June 2015 at 16:50, C. Titus Brown <span dir="ltr">&lt;<a href="mailto:ctbrown@ucdavis.edu" target="_blank">ctbrown@ucdavis.edu</a>&gt;</span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Yes, but what I really wanted to do was point Diane at a somewhat more<br>
friendly set of instructions --<br>
<br>
so, here!<br>
<br>
<a href="https://khmer-protocols.readthedocs.org/en/latest/metagenomics/index.html" target="_blank">https://khmer-protocols.readthedocs.org/en/latest/metagenomics/index.html</a><br>
<br>
More soon...<br>
<br>
best,<br>
--titus<br>
<br>
On Wed, Jun 03, 2015 at 09:47:57AM -0400, Ramakrishnan Srinivasan wrote:<br>
&gt; Titus,<br>
&gt;<br>
&gt; Would adding a False to check_space() help?<br>
&gt;<br>
&gt; --<br>
&gt; Ram<br>
&gt;<br>
&gt; On Wed, Jun 3, 2015 at 9:43 AM, C. Titus Brown &lt;<a href="mailto:ctbrown@ucdavis.edu">ctbrown@ucdavis.edu</a>&gt; wrote:<br>
&gt;<br>
&gt; &gt; Hi Diane, sorry, that&#39;s a bug in our code :(.<br>
&gt; &gt;<br>
&gt; &gt; More in a bit --<br>
&gt; &gt;<br>
&gt; &gt; best,<br>
&gt; &gt; --titus<br>
&gt; &gt;<br>
&gt; &gt; On Wed, Jun 03, 2015 at 04:40:36PM +0300, Diane Hatziioanou wrote:<br>
&gt; &gt; &gt;  Hello,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I&#39;m using khmer for the first time and have encountered a problem.<br>
&gt; &gt; &gt; I&#39;m running khmer version 1.4 and also using screed version 0.8.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I have a dataset which is metagenomic-like with excessive coverage which<br>
&gt; &gt; I<br>
&gt; &gt; &gt; have been following the basic metagenomic pipleline from the khmer page:<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;    1. Apply single-pass digital normalization: normalize-by-median to<br>
&gt; &gt; C=20<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;       2. Run filter-below-abund with C=100<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;    1. Partition reads with load-graph, etc. etc.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I did the load-graph.py, make-initial-stoptags.py and partition-graph.py<br>
&gt; &gt; &gt; successfully following the guide at<br>
&gt; &gt; &gt; <a href="http://khmer.readthedocs.org/en/latest/user/partitioning-big-data.html" target="_blank">http://khmer.readthedocs.org/en/latest/user/partitioning-big-data.html</a><br>
&gt; &gt; and<br>
&gt; &gt; &gt; no matter what I do I can not make the find-knots.py work.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; The command I give is:<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; find-knots.py -N 4 -x 8e9 C100Graph2 &amp;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I&#39;ve also tried:<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; find-knots.py C100Graph2 &amp;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; and:<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; find-knots.py C100Graph2 C100Graph2 &amp;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; the message I always get is:<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Traceback (most recent call last):<br>
&gt; &gt; &gt;   File &quot;/usr/local/bin/find-knots.py&quot;, line 154, in &lt;module&gt;<br>
&gt; &gt; &gt;     main()<br>
&gt; &gt; &gt;   File &quot;/usr/local/bin/find-knots.py&quot;, line 95, in main<br>
&gt; &gt; &gt;     check_space(infiles)<br>
&gt; &gt; &gt; TypeError: check_space() takes at least 2 arguments (1 given)<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; I can&#39;t find any help anywhere which described the two arguments needed,<br>
&gt; &gt; &gt; all the examples i&#39;ve seen only have one! Please let me know if you can<br>
&gt; &gt; see<br>
&gt; &gt; &gt; what I&#39;m doing wrongor if this is a known problem.<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Yours sincerely,<br>
&gt; &gt; &gt;<br>
&gt; &gt; &gt; Artemis<br>
&gt; &gt; &gt; --<br>
&gt; &gt; &gt; Dr Diane Hatziioanou<br>
&gt; &gt; &gt; Greek Mobile: (+30)6909403373<br>
&gt; &gt; &gt; UK Mobile: (+44)7779516625<br>
&gt; &gt; &gt; <a href="http://www.linkedin.com/in/dhatziioanou/" target="_blank">www.linkedin.com/in/dhatziioanou/</a><br>
&gt; &gt; &gt; <a href="https://twitter.com/DianeHIoanou" target="_blank">https://twitter.com/DianeHIoanou</a><br>
&gt; &gt;<br>
&gt; &gt; &gt; _______________________________________________<br>
&gt; &gt; &gt; khmer mailing list<br>
&gt; &gt; &gt; <a href="mailto:khmer@lists.idyll.org">khmer@lists.idyll.org</a><br>
&gt; &gt; &gt; <a href="http://lists.idyll.org/listinfo/khmer" target="_blank">http://lists.idyll.org/listinfo/khmer</a><br>
&gt; &gt;<br>
&gt; &gt;<br>
&gt; &gt; --<br>
&gt; &gt; C. Titus Brown, <a href="mailto:ctbrown@ucdavis.edu">ctbrown@ucdavis.edu</a><br>
&gt; &gt;<br>
&gt; &gt; _______________________________________________<br>
&gt; &gt; khmer mailing list<br>
&gt; &gt; <a href="mailto:khmer@lists.idyll.org">khmer@lists.idyll.org</a><br>
&gt; &gt; <a href="http://lists.idyll.org/listinfo/khmer" target="_blank">http://lists.idyll.org/listinfo/khmer</a><br>
&gt; &gt;<br>
<span class="HOEnZb"><font color="#888888"><br>
--<br>
C. Titus Brown, <a href="mailto:ctbrown@ucdavis.edu">ctbrown@ucdavis.edu</a><br>
</font></span></blockquote></div><br></div>