Picard HsMetrics on Whole Genome (no bedfile) need "HS_Penalty_20X & HS_LIBRARY_SIZE"
0
0
Entering edit mode
6.7 years ago
Jon17 ▴ 20

Is there a way to get these metrics on a whole genome scale?

I might do a hack where I create a bed & interval files covering the entire genome.

  • HS_Penalty_20X - "hybrid selection penalty" incurred to get 80% of target bases to 20X.
  • HS_LIBRARY_SIZE - estimated number of unique molecules in the selected part of the library.

I'm guessing for HS_LIBRARY_SIZE this is # of reads left post duplicates removed for whole genome sequencing?

But for HS_PENALTY_20X should I just calculate this manually? Using the formula in the documentation:

  • PF_ALIGNED_BASES / ($REFERENCE_SIZE * 20) = HS_PENALTY_20X.

That doesn't seem like it will work... :-p

picard gatk • 1.3k views
ADD COMMENT

Login before adding your answer.

Traffic: 2689 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6