282d648d8896a7bb45af20bbc306873f6a5ffb87 angie Wed Jan 17 11:03:32 2024 -0800 Added chrUn_KI270752v1 to hg38.{p11,p12,p13}.chromAlias.txt files after user pointed out in MLQ #32874 that it was missing. It is present in hg38.p14.chromAlias.txt and all *chromAlias.bb files. diff --git src/hg/makeDb/doc/hg38/patchUpdate.13.txt src/hg/makeDb/doc/hg38/patchUpdate.13.txt index b3379c0..65d7a18 100644 --- src/hg/makeDb/doc/hg38/patchUpdate.13.txt +++ src/hg/makeDb/doc/hg38/patchUpdate.13.txt @@ -880,15 +880,36 @@ # TODO: OMIM Genes needs liftOver to new alts and fixes (or redo from ncbiRefSeq). # OMIM Phenotypes needs liftOvers to all alts and fixes. Sometimes it spans a region larger # than an alt/fix, so maybe lower the percentage that has to map? ############################################################################## # GRC Incident Database (DONE 2021-10-12 galt) # Wait until the updated hg38 files have been pushed to RR because GRC Incident update is # automated. Then update the file used to map GRC's RefSeq accessions to our names: hgsql hg38 -NBe 'select alias,chrom from chromAlias where source = "refseq" order by alias;' \ > /hive/data/outside/grc/incidentDb/GRCh38/refSeq.chromNames.tab ############################################################################# +# Update hg38.p13.chromAlias.txt (DONE 2024-01-17 Angie) + + # In MLQ#32874, the user reported that chrUn_KI270752v1 is missing from hg38.p12.chromAlias.txt. + # That's because when the hg38.{p11,p12,p13}.chromAlias.txt files were initially created + # above ("Correctly versioned hg38.chromAlias.txt files in downloads"), + # chrUn_KI270752v1 was omitted because it had been removed from the RefSeq assembly as + # contamination. However, it is confusing for users to have a sequence in the db with no + # aliases despite it having Assembly, Ensembl and INSDC aliases. So add it back. + hgsql hg38 -NBe 'select * from chromAlias where chrom = "chrUn_KI270752v1"' +#HSCHRUN_RANDOM_CTG29 chrUn_KI270752v1 assembly +#KI270752.1 chrUn_KI270752v1 ensembl,genbank + cd /hive/data/genomes/hg38/goldenPath/bigZips + echo -e "chrUn_KI270752v1\tHSCHRUN_RANDOM_CTG29\tKI270752.1\t" >> hg38.p13.chromAlias.txt + # p13/hg38.p13.chromAlias.txt is a symlink to the one in this directory. + # p13/hg38.p13.chromAlias.bb is a file in p13/ . -- But it already has chrUn_KI270752v1 + # so it does not need to be updated, great. + bigBedToBed -chrom=chrUn_KI270752v1 p13/hg38.p13.chromAlias.bb stdout +#chrUn_KI270752v1 0 27745 chrUn_KI270752v1 HSCHRUN_RANDOM_CTG29 KI270752.1 KI270752.1 + + +#############################################################################