From 46e37c3e89e9d9bb245674d8bb8036004417ca92 Mon Sep 17 00:00:00 2001 From: Tanxiao Huang Date: Wed, 4 Nov 2015 17:13:17 +0800 Subject: [PATCH] Maybe raw authir forget to add .so I add . --- pileup2baseindel.pl | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/pileup2baseindel.pl b/pileup2baseindel.pl index 7dadecb..b2d07c0 100644 --- a/pileup2baseindel.pl +++ b/pileup2baseindel.pl @@ -72,7 +72,7 @@ my @region=(3*($i-1),3*($i-1)+1,3*($i-1)+2); my ($dp,$bases,$bq) = @dp_bases_bq[@region]; - my $str = parsePileup($bases,$bq,$BQcut,$offset); + my $str = parsePileup($ref,$bases,$bq,$BQcut,$offset); if($str ne "*"){ print $fh join "\t",($chr,$loc,$ref,$str); @@ -88,7 +88,7 @@ my $fh = $files{$i}; my @region=(3*($i-1),3*($i-1)+1,3*($i-1)+2); my ($dp,$bases,$bq) = @dp_bases_bq[@region]; - my $str = parsePileup($bases,$bq,$BQcut,$offset); + my $str = parsePileup($ref, $bases,$bq,$BQcut,$offset); if($str ne "*"){ print $fh join "\t",($chr,$loc,$ref,$str); } @@ -103,7 +103,7 @@ sub parsePileup{ - my ($bases,$bq,$BQcut,$offset) = @_; + my ($ref,$bases,$bq,$BQcut,$offset) = @_; if($bases eq "*"){ return "*";