Fixed indel calling in pileup wrapper. TO DO: Write a separate wrapper for indels only.

This commit is contained in:
Anton Nekrutenko
2010-04-14 13:52:35 -04:00
parent 96f85263de
commit 254ad33d95
+1 -1
View File
@@ -60,7 +60,7 @@ while (<IN>) {
my @qv = split //, $base_quality;
for my $base ( 0 .. @bases - 1 ) {
if ( ord( $qv[ $base ] ) - 33 >= $quality_cutoff )
if ( ord( $qv[ $base ] ) - 33 >= $quality_cutoff and $bases[ $base ] ne '*')
{
++$above_qv_bases;