SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Picard gwilymh General 6 12-12-2014 12:42 PM
About Picard Jin Pan Bioinformatics 8 07-15-2014 05:55 AM
picard HSmetrics cmccabe Bioinformatics 28 07-16-2013 02:31 PM
Picard: HsMetrics Seq84 Bioinformatics 0 04-18-2011 03:24 AM
Picard HsMetrics NEAR_BAIT_BASES interval Bruins Bioinformatics 0 01-25-2011 03:46 AM

Reply
 
Thread Tools
Old 10-25-2013, 12:58 PM   #1
sgroenewold
Junior Member
 
Location: Iowa, USA

Join Date: Jan 2013
Posts: 4
Default Filtering in Picard HSMetrics?

I've been using HSMetrics for a while to get 2X/10X/20X/30X coverage, as well as average probe and target coverage.

I've noticed lately that the 2X/10X/20X/30X numbers start to skew lower for alignments with high duplicate levels. In fact, if I run MarkDuplicates and then re-run HSMetrics, the 2X/10X/20X/30X numbers actually increase, which seems counterintuitive. On top of that, coverage counts that I generate using BEDTools coverageBed don't match the 2X/10X/20X/30X numbers or the average target coverage depth counts I get from HSMetrics.

I've searched high and low and not seen a good explanation of the inner workings of HSMetrics. Is it applying some of its own filtering? Is it perhaps sampling only a portion of the data?

My next steps are to try both HSMetrics and BEDTools with a small target set to see if I can replicate and explain the behavior. After that, I guess I'll try to walk through the Picard code.

But if someone out there knows what might be going on, I'm all ears.
sgroenewold is offline   Reply With Quote
Old 10-28-2013, 03:22 PM   #2
sgroenewold
Junior Member
 
Location: Iowa, USA

Join Date: Jan 2013
Posts: 4
Default One bump for another shot at a response.

Can anyone shed a little light?

Thanks.
sgroenewold is offline   Reply With Quote
Old 01-27-2014, 12:30 PM   #3
fashar1
Junior Member
 
Location: Baltimore, MD

Join Date: Jan 2013
Posts: 2
Default

This is very interesting. What is the Percent_Duplication on your samples? I just ran it with a couple of my samples and don't see a significant difference.
fashar1 is offline   Reply With Quote
Reply

Tags
coverage depth, picard

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 11:02 AM.


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