From 70e1ba82a4ae5f639e26365e868e3e1888934be8 Mon Sep 17 00:00:00 2001 From: borrageiros Date: Mon, 8 Jul 2024 17:03:14 +0200 Subject: [PATCH] SAVE --- change_quality.py | 16 ++++++++-------- 1 file changed, 8 insertions(+), 8 deletions(-) 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: