From f9620392732a172b74b097b1d8ee175e1af00ecf Mon Sep 17 00:00:00 2001
From: ebanks <ebanks@348d0f76-0448-11de-a6fe-93d51630548a>
Date: Mon, 18 Oct 2010 01:44:54 +0000
Subject: [PATCH] Using James P's patch for the liftover script

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@4515 348d0f76-0448-11de-a6fe-93d51630548a
---
 perl/liftOverVCF.pl | 6 +++---
 1 file changed, 3 insertions(+), 3 deletions(-)

diff --git a/perl/liftOverVCF.pl b/perl/liftOverVCF.pl
index 5a34594d1..6e3b0fccf 100755
--- a/perl/liftOverVCF.pl
+++ b/perl/liftOverVCF.pl
@@ -54,11 +54,11 @@ while ( $inHeader == 1 ) {
     }
 }
 close(UNSORTED);
-
-$cmd = "grep \"^#\" -v $unsorted_vcf | sort -n -k2 -T $tmp | $gatk/perl/sortByRef.pl --tmp $tmp - $newRef.fasta.fai";
-print SORTED `$cmd`;
 close(SORTED);
 
+$cmd = "grep \"^#\" -v $unsorted_vcf | sort -n -k2 -T $tmp | $gatk/perl/sortByRef.pl --tmp $tmp - $newRef.fasta.fai >> $sorted_vcf";
+system($cmd);
+
 # Filter the VCF for bad records
 print "\nFixing/removing bad records...\n";
 $cmd = "java -jar $gatk/dist/GenomeAnalysisTK.jar -T FilterLiftedVariants -R $newRef.fasta -B:variant,vcf $sorted_vcf -o $out";