probability of two samples sharing a SNP by error

I’m trying to calculate the probability of the same snp occurring in two different samples by chance. I’m pretty lost here, below are the parameters I have, I’m calling SNPs at ≥1% of the reads in targeted, ultradeep sequence data.

gene_length = 500
error_rate = 0.01
read_depth = 100000