# HG changeset patch
# User geert-vandeweyer
# Date 1411989431 14400
# Node ID 8262299f8f3c49fdf20a2ce485a3d65ff2530b33
# Parent  8938f339ed37996ae5573f8502cd06b11b03ae65
Fixed loading of dbsnp after bcftools-tabix switch.

diff -r 8938f339ed37 -r 8262299f8f3c DC_Genotyper.pl
--- a/DC_Genotyper.pl	Mon Sep 29 05:19:48 2014 -0400
+++ b/DC_Genotyper.pl	Mon Sep 29 07:17:11 2014 -0400
@@ -207,6 +207,13 @@
 	print SNP "chr0\t1\t.\tA\tT\n";
 	close SNP;
 }
+open IN, "$wd/dbsnp.txt";
+while (<IN>) {
+	chomp;
+	my @p = split(/\t/,$_);
+	$dbsnp{$p[0].'-'.$p[1]} = $p[3].'-'.$p[4].'-'.$p[2];
+}
+close IN;
 
 ## now process the bam file.
 mkdir "$wd/WIGS/";