68e60a3ccfe0c3aaa48f2fa85f86ca700d681e81
hiram
Tue Mar 31 12:08:22 2020 -0700
correct hub.txt not testHub.txt refs #23891
diff --git src/hg/makeDb/doc/asmHubs/mkAsmStats.pl src/hg/makeDb/doc/asmHubs/mkAsmStats.pl
index 1b713fa..5d1f5b0 100755
--- src/hg/makeDb/doc/asmHubs/mkAsmStats.pl
+++ src/hg/makeDb/doc/asmHubs/mkAsmStats.pl
@@ -1,328 +1,329 @@
#!/usr/bin/env perl
use strict;
use warnings;
use File::stat;
my $argc = scalar(@ARGV);
if ($argc != 2) {
printf STDERR "mkAsmStats Name asmName\n";
printf STDERR "e.g.: mkAsmStats Mammals mammals\n";
exit 255;
}
my $Name = shift;
my $asmHubName = shift;
my $home = $ENV{'HOME'};
my $toolsDir = "$home/kent/src/hg/makeDb/doc/asmHubs";
my $commonNameList = "$asmHubName.asmId.commonName.tsv";
my $commonNameOrder = "$asmHubName.commonName.asmId.orderList.tsv";
my @orderList; # asmId of the assemblies in order from the *.list files
# the order to read the different .list files:
my %betterName; # key is asmId, value is better common name than found in
# assembly_report
my $assemblyTotal = 0; # complete list of assemblies in this group
my $asmCount = 0; # count of assemblies completed and in the table
my $overallNucleotides = 0;
my $overallSeqCount = 0;
my $overallGapSize = 0;
my $overallGapCount = 0;
##############################################################################
# from Perl Cookbook Recipe 2.17, print out large numbers with comma delimiters:
##############################################################################
sub commify($) {
my $text = reverse $_[0];
$text =~ s/(\d\d\d)(?=\d)(?!\d*\.)/$1,/g;
return scalar reverse $text
}
##############################################################################
### start the HTML output
##############################################################################
sub startHtml() {
my $timeStamp = `date "+%F"`;
chomp $timeStamp;
my $subSetMessage = "subset of $asmHubName only";
if ($asmHubName eq "vertebrate") {
$subSetMessage = "subset of other ${asmHubName}s only";
}
print <<"END"
$Name Genomes assembly hubs, assembly statistics
Assemblies from NCBI/Genbank/Refseq sources, $subSetMessage.
Data resource links
NOTE: Click on the column headers to sort the table by that column
The link to genome browser will attach only that single assembly to
the genome browser.
END
}
##############################################################################
### start the table output
##############################################################################
sub startTable() {
print <<"END"
| count |
common name
link to genome browser |
scientific name
and data download |
NCBI assembly |
sequence
count | genome size
nucleotides |
gap
count | unknown bases
(gap size sum) | masking
percent |
|---|
END
}
##############################################################################
### end the table output
##############################################################################
sub endTable() {
my $commaNuc = commify($overallNucleotides);
my $commaSeqCount = commify($overallSeqCount);
my $commaGapSize = commify($overallGapSize);
my $commaGapCount = commify($overallGapCount);
my $percentDone = 100.0 * $asmCount / $assemblyTotal;
my $doneMsg = "";
if ($asmCount < $assemblyTotal) {
$doneMsg = sprintf(" (%d build completed, %.2f %% finished)", $asmCount, $percentDone);
}
if ($assemblyTotal > 1) {
print <<"END"
| TOTALS: | total assembly count ${assemblyTotal}${doneMsg} |
$commaSeqCount |
$commaNuc |
$commaGapCount |
$commaGapSize |
|
|---|
END
} else {
print <<"END"
END
}
} # sub endTable()
##############################################################################
### end the HTML output
##############################################################################
sub endHtml() {
if ($asmHubName ne "viral") {
- printf "\n
\n";
- printf "| Assembly hubs index pages: | \n";
+ printf "\n
\n\n";
}
print <<"END"