diff --git a/heRho/heRho_tally_pairwise_counts_vcf.py b/heRho/heRho_tally_pairwise_counts_vcf.py
index 6c18ea0..09c1d71 100644
--- a/heRho/heRho_tally_pairwise_counts_vcf.py
+++ b/heRho/heRho_tally_pairwise_counts_vcf.py
@@ -265,7 +265,7 @@ def __init__(
 
     def variant_bed_intersect(self, bed_intervals=None):
         self.bed_variant_dict = {}
-        print("Sample: %s, Heterozygosity: %s" % (self.name, int(len(snp_array))/int(snp_array[-1])))
+        print("Sample: %s, Heterozygosity: %s" % (self.name, int(len(self.snp_array))/int(self.snp_array[-1])))
         for interval_index, interval in enumerate(bed_intervals):
             mask_array = (self.snp_array >= interval.start) & (
                 self.snp_array < interval.end