1b17e173bde7713fe7080825f27efb8e952419f4 hiram Thu Aug 1 14:51:39 2019 -0700 updating text after QA comments refs #23734 diff --git src/hg/utils/automation/asmHubGc5Percent.pl src/hg/utils/automation/asmHubGc5Percent.pl index c43c8cf..f123400 100755 --- src/hg/utils/automation/asmHubGc5Percent.pl +++ src/hg/utils/automation/asmHubGc5Percent.pl @@ -1,64 +1,64 @@ #!/usr/bin/env perl use strict; use warnings; use FindBin qw($Bin); use lib "$Bin"; use AsmHub; my $argc = scalar(@ARGV); if ($argc != 3) { printf STDERR "usage: asmHubGc5Percent.pl asmId asmId.names.tab buildDir\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 buildDir is the directory with bbi/asmId.gc5Base.bw.\n"; exit 255; } my $asmId = shift; my $namesFile = shift; my $buildDir = shift; my $gc5Bw = "$buildDir/bbi/$asmId.gc5Base.bw"; if ( ! -s $gc5Bw ) { printf STDERR "ERROR: can not find CpG masked file:\n\t'%s'\n", $gc5Bw; 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 $averageGC = `bigWigInfo $gc5Bw | egrep "mean:" | sed -e 's/mean: //;'`; chomp $averageGC; $averageGC = sprintf("%.2f", $averageGC); print <<_EOF_ <h2>Description</h2> <p> The GC percent track shows the percentage of G (guanine) and C (cytosine) bases in 5-base windows on the $assemblyDate $em${organism}$noEm/$asmId/$ncbiAssemblyId genome assembly. High GC content is typically associated with gene-rich areas. The average -overall GC percent for the entire assembly is % $averageGC; +overall GC percent for the entire assembly is % $averageGC. </p> <p> This track may be configured in a variety of ways to highlight different aspects of the displayed information. Click the "Graph configuration help" link for an explanation of the configuration options. </p> <h2>Credits</h2> <p> The data and presentation of this graph were prepared by <a href="mailto:hiram@soe .ucsc.edu">Hiram Clawson</a>. </p> _EOF_ ;