From 425a8d43fcb29439856321409da4300f75060b3e Mon Sep 17 00:00:00 2001
From: sandho0000 <michael.sandholzer@unibas.ch>
Date: Wed, 9 Nov 2022 13:48:36 +0100
Subject: [PATCH] 	new file:   read_sequencer_package/read_sequence.py

---
 read_sequencer_package/read_sequence.py | 31 +++++++++++++++++++++++++
 1 file changed, 31 insertions(+)
 create mode 100644 read_sequencer_package/read_sequence.py

diff --git a/read_sequencer_package/read_sequence.py b/read_sequencer_package/read_sequence.py
new file mode 100644
index 0000000..45a5c6a
--- /dev/null
+++ b/read_sequencer_package/read_sequence.py
@@ -0,0 +1,31 @@
+'''
+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
-- 
GitLab