Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
How to calculate coverage arendon Bioinformatics 53 08-20-2015 07:23 AM
How do I calculate the probability of sequence coverage at particular window? vkpilla Bioinformatics 3 10-22-2013 03:13 AM
Hot to calculate coverage on a 454 targeted region sequenced ? Giorgio C Bioinformatics 9 10-23-2012 01:06 PM
how to calculate the genome heterozygosity luckylove De novo discovery 0 06-19-2012 05:16 PM
How to calculate the sequencing coverage from bioscope result coonya SOLiD 0 12-28-2010 10:14 PM

Thread Tools
Old 08-08-2012, 08:35 AM   #1
Location: New Hampshire

Join Date: Sep 2011
Posts: 20
Default How do calculate % of genome with >20x coverage?

Hi all,

I've searched through the forums for my question but I haven't found anything for my specific stat that I'm looking for. How do you figure out the % of your genome that had >20x coverage? I'm new to bioinformatics so scripting is not something I can do but I can run any program. Is there a program out there that can help me?


Last edited by SeqVicious; 08-08-2012 at 09:04 AM.
SeqVicious is offline   Reply With Quote
Old 08-08-2012, 09:46 AM   #2
Senior Member
Location: Aberdeen, Scotland

Join Date: Jan 2010
Posts: 388

Could use Picard HSmetrics:

I know I use it for exome data, but feed it the right interval list and I'm sure it will work with whole genome. It will give you a column in the output:

PCT_TARGET_BASES_20X: The percentage of ALL target bases acheiving 20X or greater coverage.
Bukowski is offline   Reply With Quote
Old 08-08-2012, 10:13 AM   #3
Location: California

Join Date: Sep 2008
Posts: 45

GATK DepthOfCoverage will do this.
aaronh is offline   Reply With Quote
Old 08-08-2012, 10:47 AM   #4
Location: New Hampshire

Join Date: Sep 2011
Posts: 20

Thank you. I am going to give GATK a try. Luckily I already have this installed.

SeqVicious is offline   Reply With Quote
Old 08-08-2012, 09:48 PM   #5
Location: China

Join Date: Jul 2012
Posts: 14

You can use SOAP2 or BWA for mapping the sequencing reads with the reference ,no matter DNA or RNA data .and then use SOAP.coverage to calculate the coverage and the depth.
Tong.W is offline   Reply With Quote
Old 08-09-2012, 03:56 AM   #6
Junior Member
Location: Liege, Belgium

Join Date: Jun 2011
Posts: 4

You could also use the genomeCoverageBed command from the BEDTools suite.
stephwen 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 09:40 AM.

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