87f551dfb14925f9e7e1d6f23d798fef313fd73c
hiram
Mon Sep 1 23:25:33 2025 -0700
new chain net construction and update genePred names
diff --git src/hg/utils/automation/asmHubChainNetTrackDb.pl src/hg/utils/automation/asmHubChainNetTrackDb.pl
index 9d66641e9c1..ef461ab02de 100755
--- src/hg/utils/automation/asmHubChainNetTrackDb.pl
+++ src/hg/utils/automation/asmHubChainNetTrackDb.pl
@@ -37,30 +37,31 @@
$commonName{$gcX} = $comName;
}
close ($CN);
open ($CN, "-|", "hgsql -N -e 'select name,organism from dbDb;' hgcentraltest") or die "can not hgsql -N -e 'select name,organism from dbDb;'";
while (my $line = <$CN>) {
chomp $line;
my ($gcX, $comName) = split('\t', $line);
$comName =~ s/\s\(.*//;
$commonName{$gcX} = "$comName/${gcX}";
}
close ($CN);
`mkdir -p $buildDir/bbi`;
`mkdir -p $buildDir/liftOver`;
+`mkdir -p $buildDir/quickLift`;
open (DL, "ls -d $buildDir/trackData/lastz.*|") or die "can not list $buildDir/trackData/lastz.*";
while (my $lastzDir =
) {
chomp $lastzDir;
my $queryDb = basename($lastzDir);
$queryDb =~ s/lastz.//;
my $Qdb = ucfirst($queryDb);
# push @queryList, $queryDb;
$queryPrio{$queryDb} = 100;
my $fbTxt = `ls $buildDir/trackData/lastz.${queryDb}/fb.${targetAccession}.chain${Qdb}Link.txt 2> /dev/null`;
chomp $fbTxt;
if (-s "${fbTxt}") {
my $prio = `cut -d' ' -f5 $fbTxt | tr -d '()%'`;
chomp $prio;
# $queryPrio{$queryDb} = sprintf("%d", int((100 - $prio)+0.5));
@@ -261,34 +262,45 @@
", $QueryDb, $targetDb, $N, $comName, $comName, $queryDb, $queryAsmName, $queryDate, $queryDb, $targetDb,
$QueryDb, $targetDb, $QueryDb, $queryDb, $targetDb, $queryPrio{$queryDb};
}
$N++;
}
$N = 0;
$headerOut = 0;
foreach my $queryDb (@queryList) {
$comName = $queryDb;
$comName = $commonName{$queryDb} if (defined($commonName{$queryDb}));
$QueryDb = ucfirst($queryDb);
if ( -s "$buildDir/trackData/lastz.$queryDb/axtChain/chainLiftOver${QueryDb}.bb" ) {
+ printf STDERR "### $targetDb.chainLiftOver${QueryDb}.bb\n";
`rm -f $buildDir/bbi/$targetDb.chainLiftOver${QueryDb}.bb`;
`rm -f $buildDir/bbi/$targetDb.chainLiftOver${QueryDb}Link.bb`;
`ln -s ../trackData/lastz.$queryDb/axtChain/chainLiftOver${QueryDb}.bb $buildDir/bbi/$targetDb.chainLiftOver${QueryDb}.bb`;
`ln -s ../trackData/lastz.$queryDb/axtChain/chainLiftOver${QueryDb}Link.bb $buildDir/bbi/$targetDb.chainLiftOver${QueryDb}Link.bb`;
+
+ if ( -s "$buildDir/trackData/lastz.$queryDb/axtChain/${targetAccession}.${queryDb}.quick.bb" ) {
+ printf STDERR "### $targetDb.chainQuickLift${QueryDb}.bb\n";
+ `rm -f $buildDir/quickLift/${queryDb}.bb`;
+ `rm -f $buildDir/quickLift/${queryDb}.link.bb`;
+ `ln -s ../trackData/lastz.$queryDb/axtChain/${targetAccession}.${queryDb}.quick.bb $buildDir/quickLift/${queryDb}.bb`;
+ `ln -s ../trackData/lastz.$queryDb/axtChain/${targetAccession}.${queryDb}.quickLink.bb $buildDir/quickLift/${queryDb}.link.bb`;
+ } else {
+ printf STDERR "### NOT FOUND: $buildDir/trackData/lastz.$queryDb/axtChain/${targetAccession}.${queryDb}.quick.bb\n";
+ }
if ( $queryDb !~ m/^GC/ ) {
$queryDate = `hgsql -N -e 'select description from dbDb where name="$queryDb"' hgcentraltest | sed -e 's/ (.*//;'`;
chomp $queryDate;
} else {
($queryDate, $queryAsmName) = &HgAutomate::hubDateName($queryDb);
}
if (0 == $headerOut) {
printf "
track %sChainNetViewLiftOver
shortLabel Lift over
view liftover
visibility hide
parent %sChainNet
spectrum on