084578b135956fb5fd981acfad2bd929bd8f0114 hiram Fri Oct 19 11:50:10 2012 -0700 update to handle mm10 correctly diff --git src/utils/ncbi/cytoBandVerify.pl src/utils/ncbi/cytoBandVerify.pl index 5c7abc8..86c82a4 100755 --- src/utils/ncbi/cytoBandVerify.pl +++ src/utils/ncbi/cytoBandVerify.pl @@ -29,21 +29,22 @@ while (my $line = <FH>) { chomp $line; my ($chr, $start, $end, $name, $band) = split('\s+', $line); die "start < 0 at $line" if ($start < 0); die "end > $chromInfo{$chr} at $line" if ($end > $chromInfo{$chr}); $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/); die "no zero coordinate on chrom $chr" if (!exists($zeros{$chr})); die "no end coordinate on chrom $chr" if (!exists($ends{$chr})); ++$chrCount; } print "everything checks out OK on $chrCount chroms\n";