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,83 +1,85 @@ #!/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 = ""; 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; my $geneCount = `bigBedInfo $ensGeneBbi | egrep "itemCount:|basesCovered:" | xargs echo | sed -e 's/itemCount/Gene count/; s/ basesCovered/; Bases covered/;'`; chomp $geneCount; print <<_EOF_

Description

-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.

These gene predictions were generated by Ensembl.

$geneCount

Methods

For a description of the methods used in Ensembl gene predictions, please refer to Hubbard et al. (2002), also listed in the References section below.

Credits

We would like to thank Ensembl for providing this annotation. For more information, please see: Ensembl's genome annotation page.

References

Hubbard T, Barker D, Birney E, Cameron G, Chen Y, Clark L, Cox T, Cuff J, Curwen V, Down T et al. The Ensembl genome database project. Nucleic Acids Res. 2002 Jan 1;30(1):38-41. PMID: 11752248; PMC: PMC99161

_EOF_ ;