diff --git a/pileup2baseindel.pl b/pileup2baseindel.pl index 7dadecb..084ca1a 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($bases,$bq,$BQcut,$offset,$ref); 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($bases,$bq,$BQcut,$offset,$ref); if($str ne "*"){ print $fh join "\t",($chr,$loc,$ref,$str); } @@ -103,7 +103,7 @@ sub parsePileup{ - my ($bases,$bq,$BQcut,$offset) = @_; + my ($bases,$bq,$BQcut,$offset,$ref) = @_; if($bases eq "*"){ return "*"; @@ -158,7 +158,6 @@ sub parsePileup{ my $reverse_T=0; my $reverse_C=0; my $reverse_G=0; - #start the loop for(my $i=0;$i<@base;$i++){ my $ch=$base[$i]; @@ -226,4 +225,3 @@ sub parsePileup{ $str.=$insertion."\t".$deletion."\n"; return $str; } -