diff --git a/change_quality.py b/change_quality.py index 15562d8..1b5e3da 100644 --- a/change_quality.py +++ b/change_quality.py @@ -6,16 +6,16 @@ import os ### FORWARD def generate_quality_first_bases_fw(read_length): - return ''.join(random.choice('?@ABCD') for _ in range(int(read_length * 0.02))) + return ''.join(random.choice('?@ABCD') for _ in range(int(read_length * 0.01))) def generate_quality_next_bases_fw(read_length): - return ''.join(random.choice('ABCDEG') for _ in range(int(read_length * 0.02))) + return ''.join(random.choice('ABCDEG') for _ in range(int(read_length * 0.01))) def generate_quality_middle_bases_fw(read_length): - return ''.join(random.choice('FGH') for _ in range(int(read_length * 0.10))) + return ''.join(random.choice('FGH') for _ in range(int(read_length * 0.05))) def generate_quality_last_bases_fw(read_length): - return ''.join(random.choice('HI') for _ in range((int(read_length * 0.60)))) + return ''.join(random.choice('HI') for _ in range((int(read_length * 0.70)))) def generate_quality_decreasing_200_220_fw(read_length): return ''.join(random.choice('CDEFGH') for _ in range(int(read_length * 0.10))) @@ -24,10 +24,10 @@ def generate_quality_decreasing_220_260_fw(read_length): return ''.join(random.choice('ABCDEG') for _ in range(int(read_length * 0.10))) def generate_quality_decreasing_260_290_fw(read_length): - return ''.join(random.choice('?@ABCD') for _ in range(int(read_length * 0.04))) + return ''.join(random.choice('?@ABCD') for _ in range(int(read_length * 0.02))) def generate_quality_decreasing_290_300_fw(read_length): - return ''.join(random.choice('?@AB') for _ in range(int(read_length * 0.02))) + return ''.join(random.choice('?@AB') for _ in range(int(read_length * 0.01))) def change_quality_fw(input_file, output_file): with open(input_file, 'r') as f_in, open(output_file, 'w') as f_out: @@ -58,7 +58,7 @@ def generate_quality_first_bases_rv(read_length): return ''.join(random.choice('?@ABCD') for _ in range(int(read_length * 0.01))) def generate_quality_next_bases_rv(read_length): - return ''.join(random.choice('ABCDEG') for _ in range(int(read_length * 0.02))) + return ''.join(random.choice('ABCDEG') for _ in range(int(read_length * 0.01))) def generate_quality_middle_bases_rv(read_length): return ''.join(random.choice('EFGH') for _ in range(int(read_length * 0.02))) @@ -79,7 +79,7 @@ def generate_quality_decreasing_260_290_rv(read_length): return ''.join(random.choice('?@ABC') for _ in range(int(read_length * 0.04))) def generate_quality_decreasing_290_300_rv(read_length): - return ''.join(random.choice('?@AB') for _ in range(int(read_length * 0.01))) + return ''.join(random.choice('?@AB') for _ in range(int(read_length * 0.02))) def change_quality_rv(input_file, output_file): with open(input_file, 'r') as f_in, open(output_file, 'w') as f_out: