def pairwise_snp_counts(fas, out_tsv):
    """ calculate pairwise snp overlap amongst a list of fasta using multiple
    alignment

    Parameters
    ----------
    fas: list
        a list of fasta files
    out_tsv: str
        output path of the pairwise SNP matrix
    """
