a0c8e5e673c15848e6f9ab33e7e77f32c641d995
angie
Mon Sep 18 14:38:57 2017 -0700
Big search & replace: use https instead of http for NCBI URLs. refs #17793
diff --git src/utils/geneReviews/addGeneReviewToBed.pl src/utils/geneReviews/addGeneReviewToBed.pl
index d1f8946..6377721 100755
--- src/utils/geneReviews/addGeneReviewToBed.pl
+++ src/utils/geneReviews/addGeneReviewToBed.pl
@@ -1,53 +1,53 @@
#!/usr/bin/perl
use warnings;
use strict;
sub usage() {
print "usage: ./addGRtoBed.pl dbName > outputFile\n";
}
my $argc=scalar(@ARGV);
if ($argc != 1) {
usage; die "ERROR: Please supply a database name for results.\n";
}
#get the list of (unique) of gene symbols form geneReviews table
my @geneReviews = split('\n',
`hgsql -N -e "select chrom, chromStart, chromEnd, name from geneReviews;" $ARGV[0]`);
my $grBed;
my $clickMsg = " (Click links below to search GeneReviews or GeneTests)";
my $firstTime = 1;
my $field;
my $details;
foreach $grBed(@geneReviews) {
$details = "";
my @col = split(/\t/, $grBed);
#print "Processing name: ", $col[3], "
";
my @grShort = split('\n',
`hgsql -N -e 'select grShort, diseaseID, diseaseName from geneReviewsRefGene where geneSymbol="$col[3]"' $ARGV[0]`);
$firstTime = 1;
my $count = scalar(@grShort);
my $i;
my $j;
for ($i=0; $i < $count; $i++) {
my @f5 = split(/\t/, $grShort[$i]);
if ($firstTime == 1) {
$firstTime=0;
$details = "";
$details = "
GeneReview available for " . $col[3] . ": " . $clickMsg . "
";
$details .= "
"; $details .= "Short name Disease ID GeneTests disease name
"; $details .= "-----------------------------------------------------------"; $details .= "-----------------------------------------------------------"; $details .= "----------------------------------
"; } - $details .= "" . $f5[0] . ""; + $details .= "" . $f5[0] . ""; if (length($f5[0]) <= 15) { for ($j = 0; $j < 15-length($f5[0]); $j ++ ) { $details .= " "; } } $details .= $f5[1] . " "; - $details .= "" . $f5[2] ."
"; + $details .= "" . $f5[2] ."
"; } $details .= "