046e02c14c273f26be1707e52a192323c6dfcedb hiram Mon May 17 12:49:32 2021 -0700 correct printf string no redmine diff --git src/hg/utils/automation/asmHubNcbiRefSeq.pl src/hg/utils/automation/asmHubNcbiRefSeq.pl index 60097d0..64a2be6 100755 --- src/hg/utils/automation/asmHubNcbiRefSeq.pl +++ src/hg/utils/automation/asmHubNcbiRefSeq.pl @@ -196,75 +196,75 @@ _EOF_ ; if ( -s "$trackDataDir/ncbiRefSeq/${asmId}.ncbiRefSeqCurated.stats.txt" ) { $geneStats = `cat $trackDataDir/ncbiRefSeq/${asmId}.ncbiRefSeqCurated.stats.txt | awk '{printf "%d\\n", \$2}' | xargs echo`; chomp $geneStats; ($itemCount, $basesCovered) = split('\s+', $geneStats); $percentCoverage = sprintf("%.3f", 100.0 * $basesCovered / $totalBases); $itemCount = commify($itemCount); $basesCovered = commify($basesCovered); printf <<_EOF_ <p> <b>Curated gene count: </b>$itemCount<br> <b>Bases in curated genes: </b>$basesCovered<br> -<b>Percent genome coverage: </b>% $percentCoverage<br> +<b>Percent genome coverage: </b>%% $percentCoverage<br> </p> _EOF_ } else { printf <<_EOF_ <p> <b>There are no curated gene annotations.</b> </p> _EOF_ } if ( -s "$trackDataDir/ncbiRefSeq/${asmId}.ncbiRefSeqPredicted.stats.txt" ) { $geneStats = `cat $trackDataDir/ncbiRefSeq/${asmId}.ncbiRefSeqPredicted.stats.txt | awk '{printf "%d\\n", \$2}' | xargs echo`; chomp $geneStats; ($itemCount, $basesCovered) = split('\s+', $geneStats); $percentCoverage = sprintf("%.3f", 100.0 * $basesCovered / $totalBases); $itemCount = commify($itemCount); $basesCovered = commify($basesCovered); printf <<_EOF_ <p> <b>Predicted gene count: </b>$itemCount<br> <b>Bases in genes: </b>$basesCovered<br> -<b>Percent genome coverage: </b>% $percentCoverage<br> +<b>Percent genome coverage: </b>%% $percentCoverage<br> </p> _EOF_ } else { printf <<_EOF_ <p> <b>there are no predicted gene annotations</b> </p> _EOF_ } if ( -s "$trackDataDir/ncbiRefSeq/${asmId}.ncbiRefSeqOther.stats.txt" ) { $geneStats = `cat $trackDataDir/ncbiRefSeq/${asmId}.ncbiRefSeqOther.stats.txt | awk '{printf "%d\\n", \$2}' | xargs echo`; chomp $geneStats; ($itemCount, $basesCovered) = split('\s+', $geneStats); $percentCoverage = sprintf("%.3f", 100.0 * $basesCovered / $totalBases); $itemCount = commify($itemCount); $basesCovered = commify($basesCovered); printf <<_EOF_ <p> <b>Other annotation count: </b>$itemCount<br> <b>Bases in other annotations: </b>$basesCovered<br> -<b>Percent genome coverage: </b>% $percentCoverage<br> +<b>Percent genome coverage: </b>%% $percentCoverage<br> </p> _EOF_ } printf <<_EOF_ <h2>Credits</h2> <p> This track was produced at UCSC from data generated by scientists worldwide and curated by the NCBI RefSeq project. </p> <h2>References</h2> <p> Kent WJ. <a href="https://genome.cshlp.org/content/12/4/656.full" target="_blank">BLAT - the BLAST-like alignment tool</a>.