56da8a5745b1ef79c8087ab53e3590e4395d0751
hiram
  Thu May 18 11:34:20 2017 -0700
catch up checkin no redmine

diff --git src/utils/ncbi/cytoBandVerify.pl src/utils/ncbi/cytoBandVerify.pl
index 5a08e6f..0811e8b 100755
--- src/utils/ncbi/cytoBandVerify.pl
+++ src/utils/ncbi/cytoBandVerify.pl
@@ -1,63 +1,63 @@
 #!/usr/bin/env perl
 
 ## run in the directory bed/cytoBand/ where cytoBand.bed has been
 ##	produced, verifies coordinates are OK relative to chrom.sizes
 
 #	$Id: cytoBandVerify.pl,v 1.1 2007/08/15 18:50:11 hiram Exp $
 
 use warnings;
 use strict;
 
 open (FH, "<../../chrom.sizes") or die "can not open ../../chrom.sizes";
 
 my %chromInfo;	# key is chrom name, value is size
 while (my $line = <FH>) {
     chomp $line;
     my ($chr, $size) = split('\s+',$line);
     $chromInfo{$chr} = $size;
 }
 close (FH);
 
 # foreach my $chr (sort (keys %chromInfo)) {
 #   print "$chr\t$chromInfo{$chr}\n"; 
 #}
 
 my %zeros;  # key is chrom, value unimportant, indicates start==0 found
 my %ends;   # key is chrom, value unimportant, indicates end==chromSize found
 
 open (FH, "<cytoBand.bed") or die "can not open cytoBand.bed";
 
 my $errors = 0;
 while (my $line = <FH>) {
     chomp $line;
     my ($chr, $start, $end, $name, $band) = split('\s+', $line);
     if ($start < 0) {
 	printf STDERR "start < 0 at %s\n", $line;
 	++$errors;
     }
     if ($end > $chromInfo{$chr}) {
 	printf STDERR "end > %d at %s\n", $chromInfo{$chr}, $line;
 	++$errors;
     }
     $zeros{$chr} = 1 if (0 == $start);
     $ends{$chr} = 1 if ($chromInfo{$chr} == $end);
 }
 
 my $chrCount = 0;
 #  all zeros and ends covered ?
 foreach my $chr (sort (keys %chromInfo)) {
-    next if ($chr =~ m/random/);
-    next if ($chr =~ m/chrM/);
-    next if ($chr =~ m/chrUn/);
+#    next if ($chr =~ m/random/);
+#    next if ($chr =~ m/chrM/);
+#    next if ($chr =~ m/chrUn/);
     if (!exists($zeros{$chr})) {
 	printf STDERR "no zero coordinate on chrom %s\n", $chr;
 	++$errors;
     }
     if (!exists($ends{$chr})) {
 	printf STDERR "no end coordinate on chrom %s\n", $chr;
 	++$errors;
     }
     ++$chrCount;
 }
 
 print "everything checks out OK on $chrCount chroms\n" if ($errors == 0);