3442f799e4ddde42dfc0257e514ea1831a514f97 markd Sun Apr 24 21:03:28 2022 -0700 add score to T2T supplied liftOver chains diff --git src/hg/makeDb/doc/hg19.txt src/hg/makeDb/doc/hg19.txt index 6b2a553..8a433d9 100644 --- src/hg/makeDb/doc/hg19.txt +++ src/hg/makeDb/doc/hg19.txt @@ -34770,30 +34770,31 @@ bedGraphToBigWig jarvis.bed ../chrom.sizes jarvis.bw # Go to Brian's cube, push "that was easy" button ############################################################################# # chm13 liftover alignments (2022-03-29 markd) # preliminary CHM13 <-> hg19 liftOver chains until NCBI produces consensus alignments # provide by Nae-Chyun Chen <naechyun.chen@gmail.com> mkdir -p /hive/data/genomes/hg19/bed/chm13LiftOver cd /hive/data/genomes/hg19/bed/chm13LiftOver # originally Obtain GRCh37 from T2T Globus: team-liftover/v1_nflo/hg19-chm13v2.chain # 2022-04-09 updated due to discovery that chrM was missing +# 2022-04-19 it was discover that chains render oddly due to the lack of chain ids. # obtain updated chains from cp /hive/data/genomes/asmHubs/genbankBuild/GCA/009/914/755/GCA_009914755.4_T2T-CHM13v2.0/trackData/hgLiftOver/hg19-chm13v2.over.chain.gz . # make NCBI query names as well, since chromAlias doesn't work yet for click-through # note doc is wrong in chromToPsl PSL query is 10 chainToPslBasic hg19-chm13v2.over.chain.gz stdout | chromToUcsc -k 10 -a /hive/data/genomes/asmHubs/genbankBuild/GCA/009/914/755/GCA_009914755.4_T2T-CHM13v2.0/GCA_009914755.4_T2T-CHM13v2.0.chromAlias.txt | pslToChain stdin stdout | pigz -c > hg19-chm13v2.ncbi-qnames.over.chain.gz # build bigChain files: hgLoadChain -noBin -test none bigChain hg19-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 hg19-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 hg19-chm13v2.over.link.bb