| Previous changeset 19:8938f339ed37 (2014-09-29) Next changeset 21:36e1ba3b4d1b (2016-01-07) |
|
Commit message:
Fixed loading of dbsnp after bcftools-tabix switch. |
|
modified:
DC_Genotyper.pl |
| b |
| 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/"; |