92509db0231e9b1ae83aa0c85182883dd8eabc9b hiram Wed Sep 14 13:59:35 2022 -0700 fixup chainNet trackDb and now sending files to dynablat and arranging hub files to place it in use refs #29811 diff --git src/hg/utils/automation/asmHubChainNet.pl src/hg/utils/automation/asmHubChainNet.pl index 7c6f771..a168518 100755 --- src/hg/utils/automation/asmHubChainNet.pl +++ src/hg/utils/automation/asmHubChainNet.pl @@ -1,43 +1,42 @@ #!/usr/bin/env perl use strict; use warnings; use FindBin qw($Bin); use lib "$Bin"; use AsmHub; use HgAutomate; use File::Basename; my $argc = scalar(@ARGV); if ($argc != 4) { - printf STDERR "usage: asmHubChainNet.pl asmId asmId.names.tab queryId hubPath > asmId.chainNet.html\n"; + printf STDERR "usage: asmHubChainNet.pl asmId ncbiAsmId asmId.names.tab queryId > asmId.chainNet.html\n"; printf STDERR "where asmId is the assembly identifier,\n"; printf STDERR "and asmId.names.tab is naming file for this assembly,\n"; printf STDERR "and queryId is the asmId or db of the other organism,\n"; - printf STDERR "and hubPath is the path to this assembly directory in .../hubs/.\n"; exit 255; } # specific to UCSC environment my $dbHost = "hgwdev"; my $asmId = shift; +my $ncbiAsmId = shift; my $namesFile = shift; my $queryId = shift; -my $hubUrl = shift; # if assembly hub, need to find the real full assembly ID if ($queryId =~ m/^GC/) { my $gcX = substr($queryId,0,3); my $d0 = substr($queryId,4,3); my $d1 = substr($queryId,7,3); my $d2 = substr($queryId,10,3); my $hubBuildDir = "refseqBuild"; $hubBuildDir = "genbankBuild" if ($gcX eq "GCA"); $queryId = `ls -d /hive/data/genomes/asmHubs/$hubBuildDir/$gcX/$d0/$d1/$d2/${queryId}*`; chomp $queryId; $queryId =~ s#.*/##; } my $ncbiAssemblyId = `grep -v "^#" $namesFile | cut -f10`;