As we calculate SNPs pairwise anyway, and people like to use SNP distances for various reasons, might be worth adding a flag that saves SNP output in some way. Possible options:
just save nucmer output
write out snp distance for each pair as well as dcj (will be faster as not writing all the nucmer output)
As we calculate SNPs pairwise anyway, and people like to use SNP distances for various reasons, might be worth adding a flag that saves SNP output in some way. Possible options: