53d682e17e5ae17bd10a2f240178d05b4da9d0f4 hiram Mon Feb 17 16:39:15 2020 -0800 correct behavior for trackDb for ensGene track refs #23891 diff --git src/hg/utils/automation/asmHubEnsGene.pl src/hg/utils/automation/asmHubEnsGene.pl index 09302c1..4a4b82a 100755 --- src/hg/utils/automation/asmHubEnsGene.pl +++ src/hg/utils/automation/asmHubEnsGene.pl @@ -1,59 +1,61 @@ #!/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 != 3) { - printf STDERR "usage: asmHubEnsGene.pl asmId asmId.names.tab bbi/asmId\n"; +if ($argc != 4) { + printf STDERR "usage: asmHubEnsGene.pl asmId asmId.names.tab bbi/asmId ensVersion\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 bbi/asmId is the path prefix to .ensGene.bb.\n"; + printf STDERR "the ensVersion is from trackData/ensGene/version.txt\n"; exit 255; } my $asmId = shift; my $namesFile = shift; my $bbiPrefix = shift; +my $ensVersion = shift; my $ensGeneBbi = "$bbiPrefix.ensGene.bb"; if ( ! -s $ensGeneBbi ) { printf STDERR "ERROR: can not find ensGene bbi file:\n\t'%s'\n", $ensGeneBbi; exit 255; } my $em = "<em>"; my $noEm = "</em>"; 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; my $geneCount = `bigBedInfo $ensGeneBbi | egrep "itemCount:|basesCovered:" | xargs echo | sed -e 's/itemCount/Gene count/; s/ basesCovered/; Bases covered/;'`; chomp $geneCount; print <<_EOF_ <h2>Description</h2> <p> -This track shows the Ensembl gene, version 86, annotations on +This track shows the Ensembl gene, $ensVersion, annotations on the $assemblyDate $em${organism}$noEm/$asmId genome assembly.<br> <br> These gene predictions were generated by <a href="http://www.ensembl.org/index.html" target="_blank">Ensembl</a>.<br> <br> $geneCount </p> <h2>Methods</h2> <p> For a description of the methods used in Ensembl gene predictions, please refer to <a href="https://academic.oup.com/nar/article/30/1/38/1332872/The-Ensembl-genome-database-project" target="_blank">Hubbard <em>et al</em>. (2002)</a>,