convenient tools to collect target mapping metrics
Entering edit mode
2.7 years ago
J.F.Jiang ▴ 880

Hi all,

I want to collect some metrics for my NGS panel.

To evaluate the capture efficiency of the probes, I want to calculate the on target rate from reads and base.

For reads, bedtools is applied to get the rate, however, reads can be calculate twice if they span two probes. Any method to avoid this?

To calculate the rate from base, I first get the depth file from bam using samtools. The sum of all these depth is listed as M. Then probe base depth was extracted from the depth file among those targets., which is listed as N. Base on target rate is calculated as N/M. Any quick guide for more efficient tool?

Thanks. Junfeng

target mapping • 617 views
Entering edit mode

Picard tools would be one place to start.

Entering edit mode

I've tried PICARD to get the HSMetrics from bam file. However, it is very interesting that it filter those reads not uniq mapped to the target, as PF_BASES_ALIGNED.

However, ON_BAIT_BASES and ON_TARGET_BASES are quite different, although same interval file is offered.

For example, 90% bases are uniquely mapped to the genome, 60% bases are mapped to the BAIT, 12% bases are mapped near to the BAIT.

But only 36.8% bases are mapped to the target.

I am wondering why so different results are obtained since the same interval file is provided.


Login before adding your answer.

Traffic: 1306 users visited in the last hour
Help About
Access RSS

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

Powered by the version 2.3.6