zzh
/
gatk-3.8
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

Updated to reflect Mark's renaming binge 

git-svn-id: file:///humgen/gsa-scr1/gsa-engineering/svn_contents/trunk@3697 348d0f76-0448-11de-a6fe-93d51630548a
This commit is contained in:
ebanks 2010-06-30 22:48:42 +00:00
parent 5f2b2d860e
commit 15a6be5d50
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
perl/liftOverVCF.pl
View File

@ -34,7 +34,7 @@ my $unsorted_vcf = "$tmp_prefix.unsorted.vcf";
# lift over the file
print "Lifting over the vcf...";
my $cmd = "java -jar $gatk/dist/GenomeAnalysisTK.jar -T LiftoverVCF -R $oldRef.fasta -B vcf,vcf,$in -o $unsorted_vcf -chain $chain -dict $newRef.dict";
my $cmd = "java -jar $gatk/dist/GenomeAnalysisTK.jar -T LiftoverVariants -R $oldRef.fasta -B vcf,vcf,$in -o $unsorted_vcf -chain $chain -dict $newRef.dict";
system($cmd);
# we need to sort the lifted over file now
@ -61,7 +61,7 @@ close(SORTED);
# Filter the VCF for bad records
print "\nFixing/removing bad records...\n";
$cmd = "java -jar $gatk/dist/GenomeAnalysisTK.jar -T FilterLiftedVCF -R $newRef.fasta -B vcf,vcf,$sorted_vcf -o $out";
$cmd = "java -jar $gatk/dist/GenomeAnalysisTK.jar -T FilterLiftedVariants -R $newRef.fasta -B vcf,vcf,$sorted_vcf -o $out";
system($cmd);
# clean up