diff --git a/read_sequencer_package/read_sequence.py b/read_sequencer_package/read_sequence.py
new file mode 100644
index 0000000000000000000000000000000000000000..45a5c6a6a2f9f4a8c73ba95b8c9d5d458fcc97df
--- /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