Mercurial > repos > nilesh > somatic_sniper
view somatic_sniper_wrapper.pl @ 16:1c0abb798dcb draft
Uploaded
author | nilesh |
---|---|
date | Mon, 08 Jul 2013 13:35:51 -0400 |
parents | 2e66a79fe57f |
children |
line wrap: on
line source
use strict; use warnings; use File::Basename; use Cwd; use File::Path qw(make_path remove_tree); die qq( Bad numbr of inputs ) if(!@ARGV); my $options =""; my $normal=""; my $tumor=""; my $output=""; foreach my $input (@ARGV) { my @tmp = split "::", $input; if($tmp[0] eq "NORMAL") { $normal = $tmp[1]; } elsif($tmp[0] eq "TUMOR") { $tumor = $tmp[1]; } elsif($tmp[0] eq "OUTPUT") { $output = $tmp[1]; } elsif($tmp[0] eq "OPTION") { $options = "$options ${tmp[1]}"; } else { die("Unknown Input: $input\n"); } } my $working_dir = cwd(); system ("ln -s $normal $working_dir/normal.bam"); system ("samtools index $working_dir/normal.bam"); system ("ln -s $tumor $working_dir/tumor.bam"); system ("samtools index $working_dir/tumor.bam"); system ("bam-somaticsniper $options $working_dir/tumor.bam $working_dir/normal.bam $output");