[khmer] abundance-dist for paired end data

[Fields, Christopher J]

Personally, I interleave the reads and run them together for all the analyses, then split them later (after filtering, etc).  But it would be interesting to hear of other options…

chris

On Oct 16, 2013, at 3:54 PM, Nacho Caballero <nachocab at gmail.com<mailto:nachocab at gmail.com>>
 wrote:


I’m following the STAMPS tutorial<http://ged.msu.edu/angus/2013-hmp-assembly-webinar/exploring-stamps-data.html> and I created the bloom filter for my two paired end files

python load-into-counting.py -x 1e8 -k 20 sample.kh<http://sample.kh/> sample_1.fastq sample_2.fastq

Now I’m trying to generate the distribution, but there is no -p option, how should I do it?

python abundance-dist.py sample.kh<http://sample.kh/> ??.fastq sample.hist

_______________________________________________
khmer mailing list
khmer at lists.idyll.org<mailto:khmer at lists.idyll.org>
http://lists.idyll.org/listinfo/khmer

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.idyll.org/pipermail/khmer/attachments/20131016/2f99df33/attachment-0002.htm>