SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Dictionary size problem in Picard HsMetrics cursecatcher Bioinformatics 1 05-06-2018 03:49 PM
Filtering in Picard HSMetrics? sgroenewold Bioinformatics 2 01-27-2014 11:30 AM
picard HSmetrics cmccabe Bioinformatics 28 07-16-2013 01:31 PM
Picard: HsMetrics Seq84 Bioinformatics 0 04-18-2011 02:24 AM
Picard HsMetrics NEAR_BAIT_BASES interval Bruins Bioinformatics 0 01-25-2011 02:46 AM

Reply
 
Thread Tools
Old 02-11-2020, 12:12 PM   #1
floydian_slip
Junior Member
 
Location: RTP

Join Date: Jun 2012
Posts: 6
Default Confusing Picard HsMetrics ON_BAIT_VS_SELECTED for near=250 and near=0

Hi, for an assay developed by my vendor, I ran latest Picard HsMetrics (2.21.8) with two different options for the same set of BAM files: one with a NEAR_DISTANCE=250 and one with NEAR_DISTANCE=0.
The bait and Target interval file was the same.
With the default value of NEAR_DISTANCE=250, I got ON_BAIT_VS_SELECTED value around 0.45 for my 12 samples and with NEAR_DISTANCE=0 I got it around 0.5.

ON_BAIT_VS_SELECTED is defined as the fraction of bases on or near baits that are covered by baits, ON_BAIT_BASES/(ON_BAIT_BASES + NEAR_BAIT_BASES).
In the 2nd case, since my near_distance = 0 it means that all the bases are either on the bait or outside the near distance so near_bait_bases = 0 which means that I should see ON_BAIT_VS_SELECTED = 1 but I see it around 0.5.
Am I misunderstanding anything?


The other question I have is: if both my Target and Bait files are the same, why is my mean target coverage around 2000X where as mean bait coverage is 3000X.

Thank you for your help in advance!
floydian_slip is offline   Reply With Quote
Reply

Tags
assay development, bait, hsmetrics, picard, target capture

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


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