Question

Import A List Of Regions In Dgv And Get How Many Times We Have Deletions/Duplications In These Regions

0

Entering edit mode

13.3 years ago

Omid ▴ 590

Is there a simple way to import a list of regions(hg18) like below

chr1:72794860-72856439
chr1:195005320-195066067
chr6:32605081-32633891

in DGV(database of Genomic Variants) and get how many times we have deletion or duplication in these regions?Thanks

cnv • 1.8k views

ADD COMMENT • link updated 13.3 years ago by Pierre Lindenbaum 164k • written 13.3 years ago by Omid ▴ 590

0

Entering edit mode

'import... in DGV ' or 'fetch ... from DGV' ?

ADD REPLY • link 13.3 years ago by Pierre Lindenbaum 164k

0

Entering edit mode

thanks Pierre can I import a list of regions into DGV and get output?

ADD REPLY • link 13.3 years ago by Omid ▴ 590

score 1 · Answer 1 · 2011-09-07

You can download the data from http://projects.tcag.ca/variation/tableview.asp?table=DGV_Content_Summary.txt and count the regions using, for example, bedtools.

or you can use a simple command line:

curl -s "http://projects.tcag.ca/variation/downloads/variation.hg18.v10.nov.2010.txt" |\
cut -d '       ' -f 3,4,5,6  |\
awk -F '     ' '{if($1=="chr1" &&  !(int($2)>72856439 || int($3)< 72794860 )) { print ; } }' |\
cut -d '      ' -f 4 | sort | uniq -c

  3 CopyNumber