In [ ]:

    

from skbio import DistanceMatrix
from glob import glob
from os import path
import os
import re


    

In [ ]:

    

ddir = "../../coalescent/2017-01-23_data/"
cov="32x"
var="0.1"
seeds = [s for s in os.listdir(ddir) if re.match(r'\d+', s) is not None]


    

In [ ]:

    

for seed in seeds:
    truth_f = "{d}/{s}/all_genomes-{v}.dist".format(d=ddir, s=seed, v=var)
    kwip_f = "{d}/{s}/kwip/1e8sz/{c}-{v}-ip.dist".format(d=ddir, s=seed, c=cov, v=var)
    
    td = DistanceMatrix.read(truth_f)
    print(td.ids)
    kd = DistanceMatrix.read(kwip_f)
    print(kd.ids)


    

In [ ]: