bdc5df573ffaa879405a4e13f10c205fed01c4d1 hiram Fri May 1 12:55:06 2026 -0700 explicit /cluster/bin/ path names for kent binaries to allow these scripts to function in cron jobs refs #31811 diff --git src/hg/utils/automation/asmHubAllGaps.pl src/hg/utils/automation/asmHubAllGaps.pl index f0976c83574..290d096a31f 100755 --- src/hg/utils/automation/asmHubAllGaps.pl +++ src/hg/utils/automation/asmHubAllGaps.pl @@ -1,82 +1,82 @@ #!/usr/bin/env perl use strict; use warnings; use FindBin qw($Bin); use lib "$Bin"; use AsmHub; use File::Basename; my $argc = scalar(@ARGV); if ($argc != 5) { printf STDERR "usage: asmHubAllGaps.pl asmId asmId.names.tab asmId.agp.gz hubUrl bbi/asmId > asmId.allGaps.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 asmId.agp.gz is AGP file for this assembly,\n"; printf STDERR "and hubUrl is the URL to this assembly directory for the AGP file.\n"; printf STDERR "and bbi/asmId is the path prefix to .allGaps.bb.\n"; exit 255; } my $asmId = shift; my $namesFile = shift; my $agpFile = shift; my $agpFileBase = basename($agpFile); my $hubUrl = shift; my $bbiPrefix = shift; my $allGapsBbi = "$bbiPrefix.allGaps.bb"; if ( ! -s $agpFile ) { printf STDERR "ERROR: can not find AGP file:\n\t'%s'\n", $agpFile; exit 255; } if ( ! -s $allGapsBbi ) { printf STDERR "ERROR: can not find allGaps bbi file:\n\t'%s'\n", $allGapsBbi; exit 255; } my $em = ""; my $noEm = ""; my $assemblyDate = `grep -v "^#" $namesFile | cut -f9`; chomp $assemblyDate; my $ncbiAssemblyId = `grep -v "^#" $namesFile | cut -f10`; chomp $ncbiAssemblyId; my $organism = `grep -v "^#" $namesFile | cut -f5`; chomp $organism; print <<_EOF_

Description

This track shows all gaps (any sequence of N's) in the $assemblyDate $em${organism}$noEm/$asmId genome assembly. Not all gaps were annotated in the AGP file, this track includes all sequences of N's in the assembly.

Genome assembly procedures are covered in the NCBI assembly documentation.
NCBI also provides specific information about this assembly.

Any sequence of N's in the assembly is marked as a gap in this track. The standard gap track only shows the gaps as annotated in the AGP file: $asmId.agp.gz
The NCBI document AGP Specification describes the format of the AGP file.

Gaps are represented as black boxes in this track. There is no information in this track about order or orientation of the contigs on either side of the gap.

_EOF_ ; -my $gapCount = `bigBedInfo $allGapsBbi | egrep "itemCount:|basesCovered:" | xargs echo | sed -e 's/itemCount/gap count/; s/basesCovered/bases covered/;'`; +my $gapCount = `/cluster/bin/x86_64/bigBedInfo $allGapsBbi | egrep "itemCount:|basesCovered:" | xargs echo | sed -e 's/itemCount/gap count/; s/basesCovered/bases covered/;'`; chomp $gapCount; printf "Gap count and coverage: %s\n", $gapCount; printf "

\n";