A few examples

See the ‘examples’ subdirectory for complete examples.

STAMPS data set

The ‘stamps’ data set is a fake metagenome-like data set containing two species, mixed at a 10:1 ratio. The source genomes are in ‘data/stamps-genomes.fa’. The reads file is in ‘data/stamps-reads.fa.gz’, and consists of 100-base reads with a 1% error rate.

The example shows how to construct k-mer abundance histograms, as well as the effect of digital normalization and partitioning on the k-mer abundance distribution.

See the script for running everything and the IPython Notebook.

For an overall discussion and some slides to explain what’s going on, visit the Web site for a 2013 HMP metagenome assembly webinar that Titus Brown gave.

comments powered by Disqus

Table Of Contents

Previous topic

Installing and running khmer

Next topic

An assembly handbook for khmer - rough draft

This Page

Edit this document!

This file can be edited directly through the Web. Anyone can update and fix errors in this document with few clicks -- no downloads needed.

  1. Go to A few examples on GitHub.
  2. Edit files using GitHub's text editor in your web browser (see the 'Edit' tab on the top right of the file)
  3. Fill in the Commit message text box at the bottom of the page describing why you made the changes. Press the Propose file change button next to it when done.
  4. Then click Send a pull request.
  5. Your changes are now queued for review under the project's Pull requests tab on GitHub!

For an introduction to the documentation format please see the reST primer.