SEQanswers

SEQanswers (http://seqanswers.com/forums/index.php)
-   Genomic Resequencing (http://seqanswers.com/forums/forumdisplay.php?f=28)
-   -   Panel Validation using NA12878 DNA (http://seqanswers.com/forums/showthread.php?t=75867)

NielQC 05-05-2017 03:36 AM

Panel Validation using NA12878 DNA
 
Hi there,

we have made a panel that comprises several (190) exonic alterations and now we want to test how capable are we to call these variants starting from a DNA sample. This post is for retrieve suggestions about how can we do this and/or to know if my approach is not correct.

Basically what I have been thinking about is purchasing NA12878 DNA from Coriell Institute, sequencing it and running a variant calling pipeline (without going into details, mapping vs GRCh37 and call variants with mpileup). Once I have the variants, check them vs the high-confidence calls set focusing on my 190 exonic alterations and if the alleles matches at these positions would mean that we call these variants correctly.

Is this a correct way?

Previously, aiming to evaluate just the bioinformatic pipeline and not the laboratory part, I did the same analysis with the GIAB NA12878 dataset (HiSeq x300) and checked the output variants with this callset ftp://ftp.ncbi.nlm.nih.gov/giab/ftp/...ransfer.vcf.gz available at the GIAB ftp site. Previously I normalized the variants representation of both calls sets using *vcflib vcfallelicprimitives*.

And a last question. For my validation purpose, is there any difference between purchasing the NIST reference DNA and purchasing the Coriell Institute DNA Sample? I think it's the same, but look the difference in prices...

Any suggestion will be welcome


All times are GMT -8. The time now is 05:22 AM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.