From 61a8793d24c8c2e5261443cee952a153bab82655 Mon Sep 17 00:00:00 2001
From: Heng Li <lh3@live.co.uk>
Date: Wed, 3 Jun 2009 08:40:40 +0000
Subject: [PATCH]  * in snpFilter, suppress non-SNP sites

---
 misc/samtools.pl | 1 +
 1 file changed, 1 insertion(+)

diff --git a/misc/samtools.pl b/misc/samtools.pl
index a26ac10..ac5f34a 100755
--- a/misc/samtools.pl
+++ b/misc/samtools.pl
@@ -118,6 +118,7 @@ Options: -d INT        minimum depth to call a SNP [$opts{d}]
   while (<>) {
 	my @t = split;
 	next if ($t[2] eq '*' || $hash{$t[0],$t[1]});
+	next if ($t[2] eq $t[3]);
 	my $is_good = ($t[7] >= $opts{d} && $t[7] <= $opts{D} && $t[6] >= $opts{Q} && $t[5] >= $opts{q})? 1 : 0;
 	next unless ($is_good); # drop
 	if ($t[0] ne $last_chr) { # a different chr, print
-- 
2.39.2