how to evaluate coverage at a regions of interest and generate specificity, sensitivity and PPV for that region?
0
0
Entering edit mode
5.9 years ago

Hi guys,

We have exome sequenced samples at a minimum of 15X coverage using capture kits.

  • We would like to
    • identify genes that may have low coverage because no probes are available
    • identify gene that falls within a region which is difficult to sequence by the capture kit
  • how to evaluate coverage at a regions of interest and generate specificity, sensitivity and PPV for that region?
  • how to develop statistics based on the existing samples of how well the samples provide on-target coverage, uniformity and ability to detect single nucleotide variants (SNVs), insertions and deletions (indels) and copy number variants (CNVs) in the region of interest.
  • Is there any statistics which is already developed?

Kindly experts throw some light on the above queries.

SNP DNAseq specificity sensitivity PPV • 1.1k views
ADD COMMENT
1
Entering edit mode

Do you have a bed file with exome probes? If you do, you can use bedtools to achieve several of your goals:

identify genes that may have low coverage because no probes are available

Use bedtools intersect -v with the genome annotation to get the genes / exons not covered by the probes.

how to evaluate coverage at a regions of interest

Use bedtools coverage.

ADD REPLY

Login before adding your answer.

Traffic: 2622 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