src/hg/utils/automation/makeGenomeDb.pl 1.28
1.28 2010/02/19 23:15:51 angie
Added new required setting assemblyShortLabel, an abbreviated form of the sequencing center's ID for the assembly (which is often more familiar to users), e.g. NCBI36 for hg18.
Index: src/hg/utils/automation/makeGenomeDb.pl
===================================================================
RCS file: /projects/compbio/cvsroot/kent/src/hg/utils/automation/makeGenomeDb.pl,v
retrieving revision 1.27
retrieving revision 1.28
diff -b -B -U 4 -r1.27 -r1.28
--- src/hg/utils/automation/makeGenomeDb.pl 5 Feb 2010 23:52:29 -0000 1.27
+++ src/hg/utils/automation/makeGenomeDb.pl 19 Feb 2010 23:15:51 -0000 1.28
@@ -78,9 +78,15 @@
- The month and year in which this assembly was released by the sequencing
center.
assemblyLabel XXXXXX
- - The sequencing center's label or version identifier for this release.
+ - The detailed/long form of the sequencing center's label or version
+ identifier for this release (e.g. 'Genome Center at Washington University,
+ St. Louis, Genus_species 1.2.3').
+
+assemblyShortLabel XXXXXX
+ - The abbreviated form of the sequencing center's label or version identifier
+ for this release (e.g. 'WUGSC 1.2.3').
orderKey NN
- A priority number (for the central database's dbDb.orderKey column)
that will determine db's relative position in the assembly menu.
@@ -178,9 +184,9 @@
$opt_noGoldGapSplit
/;
# Required config parameters:
-my ($db, $scientificName, $assemblyDate, $assemblyLabel, $orderKey,
+my ($db, $scientificName, $assemblyDate, $assemblyLabel, $assemblyShortLabel, $orderKey,
$mitoAcc, $fastaFiles, $dbDbSpeciesDir, $taxId);
# Conditionally required config parameters:
my ($fakeAgpMinContigGap, $fakeAgpMinScaffoldGap,
$clade, $genomeCladePriority);
@@ -247,8 +253,9 @@
$db = &requireVar('db', \%config);
$scientificName = &requireVar('scientificName', \%config);
$assemblyDate = &requireVar('assemblyDate', \%config);
$assemblyLabel = &requireVar('assemblyLabel', \%config);
+ $assemblyShortLabel = &requireVar('assemblyShortLabel', \%config);
$orderKey = &requireVar('orderKey', \%config);
$mitoAcc = &requireVar('mitoAcc', \%config);
$fastaFiles = &requireVar('fastaFiles', \%config);
$dbDbSpeciesDir = &requireVar('dbDbSpeciesDir', \%config);
@@ -389,9 +396,9 @@
my $acat = "cat";
my $fcat = "cat";
my $sli = "";
- if ($subsetLittleIds eq "Y") {
+ if (defined $subsetLittleIds && $subsetLittleIds eq "Y") {
$sli = "-1 ";
}
foreach my $file (`ls $fastaFiles 2> /dev/null`) {
if ($file =~ m/\.gz$/) {
@@ -839,9 +846,9 @@
(name, description, nibPath, organism,
defaultPos, active, orderKey, genome, scientificName,
htmlPath, hgNearOk, hgPbOk, sourceName, taxId)
VALUES
- ("$db", "$assemblyDate", "$HgAutomate::gbdb/$db", "$genome",
+ ("$db", "$assemblyDate ($assemblyShortLabel/$db)", "$HgAutomate::gbdb/$db", "$genome",
"$defaultPos", 1, $orderKey, "$genome", "$scientificName",
"$HgAutomate::gbdb/$db/html/description.html", 0, 0, "$assemblyLabel",
$taxId);
_EOF_
@@ -1072,12 +1079,15 @@
}
close($fh);
$fh = &HgAutomate::mustOpen(">$topDir/html/gap.html");
+ my $em = $commonName ? "" : "<em>";
+ my $noEm = $commonName ? "" : "</em>";
if ($gotAgp) {
print $fh <<_EOF_
<H2>Description</H2>
-This track depicts gaps in the assembly.
+This track depicts gaps in the draft assembly ($assemblyDate, $assemblyLabel)
+of the $em\$organism$noEm genome.
*** Developer: remove this statement if no future assemblies are expected:
Many of these gaps — with the
@@ -1142,9 +1152,10 @@
;
} else {
print $fh <<_EOF_
<H2>Description</H2>
-This track depicts gaps in the assembly.
+This track depicts gaps in the draft assembly ($assemblyDate, $assemblyLabel)
+of the $em\$organism$noEm genome.
*** Developer: remove this statement if no future assemblies are expected:
Many of these gaps — with the
@@ -1183,15 +1194,14 @@
}
close($fh);
$fh = &HgAutomate::mustOpen(">$topDir/html/gold.html");
- my $em = $commonName ? "" : "<em>";
- my $noEm = $commonName ? "" : "</em>";
if ($gotAgp) {
print $fh <<_EOF_
<H2>Description</H2>
<P>
-This track shows the draft assembly of the $em\$organism$noEm genome.
+This track shows the draft assembly ($assemblyDate, $assemblyLabel)
+of the $em\$organism$noEm genome.
*** Developer: check if this is accurate:
Whole-genome shotgun reads were assembled into contigs. When possible,
@@ -1215,9 +1225,10 @@
} else {
print $fh <<_EOF_
<H2>Description</H2>
<P>
-This track shows the draft assembly of the $em\$organism$noEm genome.
+This track shows the draft assembly ($assemblyDate, $assemblyLabel)
+of the $em\$organism$noEm genome.
*** Developer: check if this is accurate:
Whole-genome shotgun reads were assembled into contigs. When possible,