1313ab4f97f5d5c1f2edd63533dce23ad8426ab7 braney Wed Oct 19 10:50:11 2022 -0700 draw zoom indicator box in the right place and make it transparent so you can see the underlying pixel diff --git src/hg/utils/otto/clinvar/clinvarSubLolly src/hg/utils/otto/clinvar/clinvarSubLolly new file mode 100755 index 0000000..a31eef6 --- /dev/null +++ src/hg/utils/otto/clinvar/clinvarSubLolly @@ -0,0 +1,20 @@ + +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"
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 + +bigBedToBed /gbdb/hg38/bbi/clinvar/clinvarMain.bb stdout | tawk '{print $40, $1,$2,$2+1,$4,$13,$15,$18,$19}' | sort -S 40g > sort.main.bed +hgsql hg38 -Ne "select * from clinvarSub" | sort -S 40g > clinvarSubSub.txt +join -t $'\t' sort.main.bed clinvarSubSub.txt | tawk '{print $2,$3,$4,$5,0,"+",0,0,"0,0,0",$6,$20,$8, $9,$1,$10,$7,$11,$12,$13,$14,$15,$16,$17,$18,$19,$21}' | sort -S 40g -k1,1 -k2,2n | tawk -f assignScore > bigBedInput.bed + +bedToBigBed -as=clinvarSubBB.as -type=bed9+11 -tab bigBedInput.bed /cluster/data/hg38/chrom.sizes clinvarSub.bb +ln -s `pwd`/clinvarSub.bb /gbdb/hg38/clinvarSubLolly/clinvarSub.bb