basic rule to remove ancestral SNPs from st398

[Perugolate]

ancestor sample is:
52,A-C,Ancestor CC398,population

[Perugolate]

Originally thought to mask these (it's a 60-bp region that's causing the problem) using --mask but it isn't available in the version that was used for original analysis. Will try altering the bases directly in the reference.