borrageiros/script_ainhoinha
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

SAVE

This commit is contained in:
borrageiros 2024-07-08 17:03:14 +02:00
parent f688641848
commit 70e1ba82a4
1 changed files with 8 additions and 8 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

16
change_quality.py
View File

@ -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: