99856cacdbaa78c1e10650c00ee8cff19bc9c0da braney Tue Oct 11 10:37:37 2022 -0700 update clinVarSub track diff --git src/hg/makeDb/doc/hg38/clinvarSubLolly.txt src/hg/makeDb/doc/hg38/clinvarSubLolly.txt index 4054a20..ed72158 100644 --- src/hg/makeDb/doc/hg38/clinvarSubLolly.txt +++ src/hg/makeDb/doc/hg38/clinvarSubLolly.txt @@ -1,18 +1,18 @@ ############################################################################# -# clinvarSubLolly track DONE BRANEY 12/14/2020 +# clinvarSubLolly track DONE BRANEY 10/8/2022 mkdir /cluster/data/hg38/bed/clinvarSubLolly cd /cluster/data/hg38/bed/clinvarSubLolly bigBedToBed /gbdb/hg38/bbi/clinvar/clinvarMain.bb stdout | tawk '{print $40, $1,$2,$2+1,$4}' | sort -S 40g > sort.main.bed hgsql hg38 -Ne "select varId,clinSign,scv from clinvarSub" | sort -S 40g > clinvarSubSub.txt join -t $'\t' sort.main.bed clinvarSubSub.txt | tawk '{print $2,$3,$4,$5,$6,$1, $7}' | sort -S 40g -k1,1 -k2,2n -k5,5 | tawk -f makeFranklin | tawk -f assignColors > tmp1 # add the line break after v409 #tawk '{print $1":"$2 + 1"-"$3" <br>Variants (submissions):"$11}' tmp1 > tmp2 tawk '{print $1":"$2 + 1"-"$3" Variants (submissions):"$11}' tmp1 > tmp2 paste tmp1 tmp2 > bigBedInput.bed bedToBigBed -as=$HOME/kent/src/hg/lib/clinvarSubLolly.as -type=bed9+5 -tab bigBedInput.bed /cluster/data/hg38/chrom.sizes clinvarSubLolly.bb mkdir -p /gbdb/hg38/clinvarSubLolly ln -s `pwd`/clinvarSubLolly.bb /gbdb/hg38/clinvarSubLolly/clinvarSubLolly.bb