Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
Overlap Graph navin_elango Bioinformatics 2 05-20-2013 11:56 AM
how to calculate genome coverage correcting for multiple genome hits alice.lunardon Bioinformatics 0 04-30-2013 03:03 AM
The way to normalize RNA-SEQ coverage data for multiple samples? puggie RNA Sequencing 8 04-09-2013 04:23 AM
overlap between multiple de novo assembled samples boetsie Bioinformatics 1 12-13-2011 09:08 AM
compare distribution of coverage in solexa data scuellar Genomic Resequencing 3 03-10-2010 10:36 AM

Thread Tools
Old 09-26-2013, 09:58 AM   #1
Junior Member
Location: Bloomington, IN

Join Date: Sep 2013
Posts: 2
Default Need to compare multiple mpileups for coverage overlap

I have 11 mpileups generated by novoalign from whole genome sequencing of Daphnia, which has a genome size of ~200MB. Each is the result of alignment of about 130M paired-end illumina reads. I find that all of the mpileups have about 110M lines, implying that we have info for about 55% of the reference. I want to quantify the degree of overlap in coverage among the mpileups. Is there an easy way to do this with samtools or something similar? Also, is there a way to construct an mpileup file that has a line for every site in the reference, rather than just those with read coverage?
kspitze is offline   Reply With Quote
Old 09-26-2013, 10:50 AM   #2
Senior Member
Location: Birmingham, UK

Join Date: Jul 2009
Posts: 356

Look into the bedtools package for easy ways to do things like this:
nickloman is offline   Reply With Quote

coverage, mpileups, novoalign, samtools

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:17 AM.

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