How do calculate % of genome with >20x coverage?

[SeqVicious]

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?

Thanks.

[Bukowski]

Could use Picard HSmetrics:

http://picard.sourceforge.net/picard...html#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.

[aaronh]

GATK DepthOfCoverage will do this. http://www.broadinstitute.org/gsa/ga...ageWalker.html

[SeqVicious]

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

Cheers.

[Tong.W]

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.

[stephwen]

You could also use the genomeCoverageBed command from the BEDTools suite.