Skip to content
Snippets Groups Projects

read_sequencer new file: read_sequencer_package/read_sequence.py

Merged Michael Sandholzer requested to merge read_sequencer into main
1 file
+ 31
0
Compare changes
  • Side-by-side
  • Inline
+ 31
0
'''
This function reads sequences
arguments: seq is a list of sequences
padding_probabilities is a number??
returns sequenced element
'''
import random
def read_sequence(seq, padding_probabilities, read_length):
reading_element = random.choice(seq)
bases =["A", "T", "C", "G"]
if read_length > len(reading_element):
for nt in [0:len(reading_element)]:
sequenced += reading_element[nt]
for nt2 in [len(reading_element):read_length]:
sequenced += random.choice(bases)
else:
for nt in [0:read_length]
sequenced += reading_element[nt]
return sequenced
\ No newline at end of file
Loading