We've moved!
This site is now read-only. You can find our new documentation site and support forum for posting questions here.
Be sure to read our welcome blog!

Picard CollectWgsMetrics for WES?

Hello,

I want to get a coverage report for WES using Picard CollectWgsMetrics, since the DepthofCoverage function in GATK is relatively slow. But the coverage results differ a lot from the 2 programs.

Settings:
* CollectWgsMetrics and DepthofCoverage wereprovided with the interval list file from the WES.
* ran with default settings.

Results:
CollectWgsMetrics reports mean coverage = 33.50, while DepthofCoverage reports 58.33. So I tried to find out why such big difference. I found that over 29.4% aligned bases were excluded in PCT_EXC_DUPE based on CollectWgsMetrics. However, only 8.9% of reads are filtered by DuplicateReadFilter in DepthofCoverage log file.

Is the lower coverage in CollectWgsMetrics due to the larger number of duplicated bases?

DoesCollectWgsMetrics and DepthofCoverage count duplicated reads differently?

Thanks a lot,
JL

Answers

Sign In or Register to comment.