def Simple_Percent_Puller(fastq_dict, num_needed): #TODO this is incomplete needed_range = range(1, num_needed + 1) stripped_dict_count = 0 strip_num = 0 base_location_chance = {} range_dict = {} percent_error_list = [] stripped_dict = name_stripper(mpileup_percent_correct(fastq_dict)) for num in stripped_dict: stripped_dict_count += 1 range_dict[stripped_dict_count] = stripped_dict[num] for num in range_dict: strip_num += 1 float_marker = float(strip_num / stripped_dict_count) base_location_chance[num] = float_marker base_location_chance_clone = base_location_chance for number in needed_range: for base in base_location_chance_clone: rand_float = float(random.random()) if rand_float <= base_location_chance_clone[base]: percent_error_list.append(range_dict[base]) base_location_chance_clone = base_location_chance break return percent_error_list
def Simple_Percent_Puller(fastq_dict, num_needed): #TODO this is incomplete needed_range = range(1, num_needed+1) stripped_dict_count = 0 strip_num=0 base_location_chance = {} range_dict = {} percent_error_list = [] stripped_dict = name_stripper(mpileup_percent_correct(fastq_dict)) for num in stripped_dict: stripped_dict_count += 1 range_dict[stripped_dict_count] = stripped_dict[num] for num in range_dict: strip_num += 1 float_marker = float(strip_num/stripped_dict_count) base_location_chance[num] = float_marker base_location_chance_clone = base_location_chance for number in needed_range: for base in base_location_chance_clone: rand_float = float(random.random()) if rand_float <= base_location_chance_clone[base]: percent_error_list.append(range_dict[base]) base_location_chance_clone = base_location_chance break return percent_error_list
def percent_puller(fastq_dict, num_needed): needed_range = range(1, num_needed + 1) stripped_dict_count = 0 base_location_chance = {} range_dict = {} stripped_dict = name_stripper(mpileup_percent_correct(fastq_dict)) for num in needed_range: float_marker = float(num / num_needed) base_location_chance[num] = float_marker for num in stripped_dict: stripped_dict_count += 1 range_dict[stripped_dict_count] = stripped_dict[num] for base in needed_range: rand_float = float(random.random()) if rand_float <= base_location_chance[base]: print base
def percent_puller(fastq_dict, num_needed): needed_range=range(1, num_needed+1) stripped_dict_count=0 base_location_chance={} range_dict={} stripped_dict=name_stripper(mpileup_percent_correct(fastq_dict)) for num in needed_range: float_marker=float(num/num_needed) base_location_chance[num]=float_marker for num in stripped_dict: stripped_dict_count+=1 range_dict[stripped_dict_count]=stripped_dict[num] for base in needed_range: rand_float=float(random.random()) if rand_float <= base_location_chance[base]: print base
def old_percent_puller(fastq_dict, num_needed): """ Creates a dictionary for use of the error_producer to fit the error probability of a real fastq :param fastq_dict: a dictionary which uses the position of base pairs and it's assinged chance of being a correct base pair :param num_needed: the amount of error percentages needed :return: dictionary with an assigned intiger and a percentage chance for a base pair at a location to be correct or incorrect """ needed_range = range(1, num_needed + 1) range_dict = {} count = 1 base_location_chance = {} percent_error_list = [] stripped_dict = name_stripper(mpileup_percent_correct(fastq_dict)) for num in stripped_dict: base_location_chance[count] = stripped_dict[num] count += 1 for num in base_location_chance: if num in needed_range: range_dict[num] = base_location_chance[num] percent_error_list.append(range_dict[num]) return percent_error_list
def old_percent_puller(fastq_dict, num_needed): """ Creates a dictionary for use of the error_producer to fit the error probability of a real fastq :param fastq_dict: a dictionary which uses the position of base pairs and it's assinged chance of being a correct base pair :param num_needed: the amount of error percentages needed :return: dictionary with an assigned intiger and a percentage chance for a base pair at a location to be correct or incorrect """ needed_range = range(1, num_needed+1) range_dict = {} count = 1 base_location_chance = {} percent_error_list = [] stripped_dict = name_stripper(mpileup_percent_correct(fastq_dict)) for num in stripped_dict: base_location_chance[count] = stripped_dict[num] count += 1 for num in base_location_chance: if num in needed_range: range_dict[num] = base_location_chance[num] percent_error_list.append(range_dict[num]) return percent_error_list