Tool that estimates coverage (and genome size) of dna sequence from reads.

This model is further extended with models for polymorphism in https://github.com/wernerkrampl/covest.

• python 3.4+
• python3-dev
• gcc

We suggest to install CovEst in python3 virtual environment.

pip install covest

pip install -e . from the project directory

type covest --help for the usage.

covest histogram -m model -k K -r read_length

• You can specify the read file using -s reads.fa parameter for more precise genome size computation.
• default K is 21
• default read length is 100
• currently, the supported models are:
  • basic: for simple genomes without repeats
  • repeat: for genomes with repetitive sequences

The input histogram can be generated from the read data using jellyfish.

• jellyfish count -m K -C reads.fa -o table.jf
• jellyfish histo table.jf -o reads.hist

The format of the histogram is just list of lines. Each lines contains an index and value separated by space.

CovEst outputs it's results in simple subset of YAML format for best human readability and possibility of machine processing.

The output are lines containing key: value. The most important keys are coverage and genome_size (or genome_size_reads if reads size was specified).

• geset.py tool for estimation genome size from reads size and known coverage
• reads_size.py tool for computation of the total reads size
• kmer_hist.py custom khmer histogram computation, it is much slower than other tools, so use it only if you have no other option.
• read_sampler.py script for subsampling reads, useful if you have very high coverage data and want to make it smaller.
• fasta_length.py get total length of all sequences in fasta file.

CovEst is licenced under GNU GPLv3 license.

CovEst is research software, so you should cite us when you use it in scientific publications!

Hozza, M., Vinař, T., & Brejová, B. (2015, September). How Big is that Genome? Estimating Genome Size and Coverage from k-mer Abundance Spectra. In String Processing and Information Retrieval (pp. 199-209). Springer International Publishing.