Latest
Open
RNA-Seq
ChIP-Seq
SNP
Assembly
Tutorials
Tools
Jobs
Forum
Planet
All »
Home
Log In
Question: Number of SNPs by target sequencing
0
gravatar for MAPK
2.8 years ago by
MAPK • 1.4k
United States
MAPK • 1.4k wrote:

I am trying to find the number of SNPs captured by each of these sequencing panel listed here: http://www.illumina.com/products/trusight-panels.html . Can someone please tell me how I can determine the number of SNPs captured by each of these targeted method?

snp target panel sequencing • 958 views
ADD COMMENT • link •
modified 2.8 years ago by Nicola Casiraghi • 440 • written 2.8 years ago by MAPK • 1.4k
2
gravatar for Nicola Casiraghi
2.8 years ago by
Nicola Casiraghi • 440
Trento, IT
Nicola Casiraghi • 440 wrote:

First, you should retrieve the BED file with targeted regions for each panel: for example for panel TruSight One Sequencing Panel Kits you can find the link to BED file Access the gene list and other product files ยป in section Product Configurations. The BED file for this panel is named TruSight One File for UCSC Browser (BED Format)

Once you retrieved your panel in BED format, download the preferred version of the dbSNP catalogue (VCF or BED format) for the Human genome at ftp://ftp.ncbi.nih.gov/snp/organisms/

Finally you can exploit bedtools to "determine number of SNPs captured by each of these targeted method" : 

bedtools intersect -a dbSNP  -b mypanel.bed
ADD COMMENT • link written 2.8 years ago by Nicola Casiraghi • 440
1

Since a lot of the TruSight panels are designed for cancer, you might not be interested in all SNPs, but only cancer-related SNPs. If this is the case, you might want to use COSMIC instead:

  1. Download the Illumina TruSight bed as Nicola mentions above.
  2. Download the COSMIC coding mutations VCF from http://cancer.sanger.ac.uk/cosmic/download (requires registration)
  3. (optional) A few of the TruSight panels cover non-coding regions as well. If this is the case for your panel, you may want to also download the non-coding COSMIC VCF, and merge the two together.
  4. Run bedtools intersect with the COSMIC VCF in place of dbSNP above.
ADD REPLY • link modified 2.8 years ago • written 2.8 years ago by Keith Callenberg • 890

Correct, or refer to the ClinVar dataset from NCBI.

ADD REPLY • link modified 2.8 years ago • written 2.8 years ago by Nicola Casiraghi • 440
Please log in to add an answer.

Similar posts • Search »

Content
Help
Access
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2572 users visited in the last hour