Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
N50 explained maasha Bioinformatics 20 09-12-2021 02:45 PM
Contig length, k-mer coverage, and differential expression nbogard General 3 09-10-2013 12:30 PM
k-mer size impacts coverage distribution (animated gif inside!) seb567 Bioinformatics 0 11-06-2010 06:20 PM
ABySS contig coverage? jgibbons1 Bioinformatics 2 09-24-2010 07:56 AM
PubMed: Digital PCR provides sensitive and absolute calibration for high throughput s Newsbot! Literature Watch 0 03-21-2009 06:03 AM

Thread Tools
Old 12-19-2010, 01:45 PM   #1
Junior Member
Location: Germany

Join Date: Feb 2008
Posts: 5
Default absolute k-mer coverage explained (Abyss)

The (very limited) Abyss manual explains the fasta headers in the contigs assembly outfile like this:

>n iii jjjj
Where n is the numeric contig ID, iii is the contig length in nucleotides, and jjjj is the absolute k-mer coverage.

I don't understand what "absolute k-mer coverage" means. Can someone maybe explain? But what I would really like to know is which basepair in a contig are represented by how many reads. IS there a way to get this from the k-mer coverage?

harrb is offline   Reply With Quote
Old 12-21-2010, 02:01 PM   #2
Senior Member
Location: NL, Leiden

Join Date: Feb 2010
Posts: 245

Hi harrb,

have a look at this older thread about the header;

About the coverage (per position on the contig i think you mean?) you can probably run ABYSS with the option


This option (as stated at records the k-mer coverage histogram in FILE. Though, I'm not 100% sure, since i can not test it @tm .

Otherwise, you can map your reads to the contigs (Bowtie, BWA), producing a .sam file. Convert the .sam file to .bam file and do;

samtools pileup -f your_ref.fa your.bam

Have a look at this thread;

Good luck,
boetsie is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 12:40 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2021, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO