The coalescent

       Hello world
    

Parse index numbers from loci


def get_ref_locus_idxs(self):
    idxs = []

    with open(self.data) as indata:
        liter = (indata.read().strip().split("|\n"))

    for idx, loc in enumerate(liter):
        lines = loc.split("\n")
        snpline = loc[-1]
        locidx, chidx, pos = snpline.split("|")[1].split(":")            
        names = [i.split()[0] for i in lines[:-1]]

        ## check coverage
        coverage = 0
        for node in self.imap:
            mincov = self.minmap[node]
            if sum([i in names for i in self.imap[node]]) >= mincov:
                coverage += 1
        
        if coverage == len(self.imap.keys()):
            pos1, pos2 = pos.split('-')
            refinfo = (idx, chidx, pos1, pos2)
            idxs.append(refinfo)
    return idxs
    

def get_ref_locus_idxs(self):
    idxs = []

    with open(self.data) as indata:
        liter = (indata.read().strip().split("|\n"))

    for idx, loc in enumerate(liter):
        lines = loc.split("\n")
        snpline = loc[-1]
        locidx, chidx, pos = snpline.split("|")[1].split(":")            
        names = [i.split()[0] for i in lines[:-1]]

    

Eaton Lab China expedition 2018

Limits on the application of genomics to phylogeny

Hello world

def wat(x, y):
    for i in range(x):
        y += x
        return y
    
Vertical slide 1
Vertical slide 2

Slide 1

2018 sampling in Yunnan and Sichuan provinces, China

P. integrifolia

return

when curried, will become

a chain of functions, each

with a single argument

function(x) {

return function(y) {

return x + y;

};

}

Pedicularis

Fading

hello

heereree

YEEP