e5f2b5cb689c5bd5c71872a400884d70cd4d2413 markd Tue Mar 29 13:05:59 2022 -0700 Added liftover chains for T2T CHM13. Clicking on these currently fails due to chromAlias issues. diff --git src/hg/makeDb/doc/hg38/chm13LiftOver.txt src/hg/makeDb/doc/hg38/chm13LiftOver.txt new file mode 100644 index 0000000..1524269 --- /dev/null +++ src/hg/makeDb/doc/hg38/chm13LiftOver.txt @@ -0,0 +1,38 @@ +############################################################################# +# chm13 liftover alignments (2022-03-29 markd) + +# preliminary CHM13 <-> hg38 liftOver chains until NCBI produces consensus alignments +# provide by Nae-Chyun Chen <naechyun.chen@gmail.com> + + mkdir -p /hive/data/genomes/hg38/bed/chm13LiftOver + cd /hive/data/genomes/hg38/bed/chm13LiftOver + +# Obtain GRCh38 from T2T Globus: team-liftover/v1_nflo/grch38-chm13v2.chain + +# rename to better match UCSC convetions and compress + mv grch38-chm13v2.chain hg38-chm13v2.over.chain + pigz hg38-chm13v2.over.chain + +# build bigChain files: + hgLoadChain -noBin -test none bigChain hg38-chm13v2.over.chain.gz + sed 's/\.000000//' chain.tab | awk 'BEGIN {OFS="\t"} {print $2, $4, $5, $11, 1000, $8, $3, $6, $7, $9, $10, $1}' > bigChainIn.tab + bedToBigBed -type=bed6+6 -as=${HOME}/kent/src/hg/lib/bigChain.as -tab bigChainIn.tab ../../chrom.sizes hg38-chm13v2.over.chain.bb + tawk '{print $1, $2, $3, $5, $4}' link.tab | csort -k1,1 -k2,2n --parallel=64 > bigLinkIn.tab + bedToBigBed -type=bed4+1 -as=${HOME}/kent/src/hg/lib/bigLink.as -tab bigLinkIn.tab ../../chrom.sizes hg38-chm13v2.over.link.bb + + rm *.tab + +# link to gbdb + mkdir -p /gbdb/hg38/bbi/chm13LiftOver + ln -sf $(pwd)/*.bb /gbdb/hg38/bbi/chm13LiftOver/ + +# make downloads, can't add to liftOver directory due to license in that directory + mkdir -p /usr/local/apache/htdocs-hgdownload/goldenPath/hg38/chm13LiftOver/ + ln -sf $(pwd)/hg38-chm13v2.over.chain.gz /usr/local/apache/htdocs-hgdownload/goldenPath/hg38/chm13LiftOver/ + +# push these files + /gbdb/hg38/bbi/chm13LiftOver/hg38-chm13v2.over.chain.bb + /gbdb/hg38/bbi/chm13LiftOver/hg38-chm13v2.over.link.bb + /usr/local/apache/htdocs-hgdownload/goldenPath/hg38/chm13LiftOver/hg38-chm13v2.over.chain.gz + +# NOTES: chromAlias doesn't work on click through yet, so disable for now