INNER CODE UNIT · Python
copy_num
BirolLab/straglr · straglr_compare.py:74
copy_num = cols[-6]
if old_version:
size = cols[-3]
copy_num = cols[-4]
if use_size:
alleles[locus][allele].append(size)
else:
alleles[locus][allele].append(copy_num)
return create_straglr_bed(alleles)
def vs_each_control(test_bed, control_bed, pval_cutoff, min_expansion=0, min_support=0, label=None):
expanded_loci = {}
new_loci = []
common_loci = []
for cols in test_bed.intersect(control_bed, loj=True, wao=True, f=0.9):
if cols[-2] == '.':