def generate_sequence(haplotype_length, total_number, mutation_rate, recombination_rate): global original_haplotype global mutated_haplotypes global recombinated_haplotypes logger = logging.getLogger('control.generate_sequence') logger.setLevel(logging.INFO) logger.debug("generator.generate_haplotype(%d, %d, %f, %f)" % (haplotype_length, total_number, float(mutation_rate)/100, float(recombination_rate)/100)) original_haplotype, mutated_haplotypes, recombinated_haplotypes = generator.generate_haplotypes(haplotype_length, total_number, float(mutation_rate)/100, float(recombination_rate)/100)
def test_generate_haplotypes(self): haplotype_length = random.randint(1,500) total_length = random.randint(1,1000) mutation_rate = random.random() recombination_rate = 1.0 - mutation_rate original, mutated_list, recombinant_list = generator.generate_haplotypes(haplotype_length, total_length, mutation_rate, recombination_rate) self.assertIsNotNone(original) self.assertIsNotNone(mutated_list) self.assertIsNotNone(recombinant_list)