#!/usr/bin/python
import IO

def mosaic_similarity(mosaic):
    genomes = list(set([row[0] for row in mosaic]))
    genomes.sort()
    pairs = [(g1, g2) for g1 in genomes for g2 in genomes[1:]]
    


mosaicfile = "~/DPhil/Recombination/Results2/Mosaics/MaskedTransposons/mosaic10.lowcov.raw.txt"
mosaic = IO.readtable(mosaicfile, ignoreheader=T)
mosaic = [row for row in mosaic if "MAGIC." in row[0]]