diff --git a/images/Ahmed_mahmoud_git_tutorial.PNG b/images/Ahmed_mahmoud_git_tutorial.PNG
new file mode 100644
index 0000000000000000000000000000000000000000..de3936f17494e54ec1e2899e0cd1f2a7b557619b
Binary files /dev/null and b/images/Ahmed_mahmoud_git_tutorial.PNG differ
diff --git a/images/markdown-tutorial-AM.PNG b/images/markdown-tutorial-AM.PNG
new file mode 100644
index 0000000000000000000000000000000000000000..f23c78a1c79d99e264298cb4fdf6850781fbd63e
Binary files /dev/null and b/images/markdown-tutorial-AM.PNG differ
diff --git a/sequence_extractor/exon_concatenation.py b/sequence_extractor/exon_concatenation.py
index 9b2465f1d16bab3bbda22a7b2cbc201a26d156d3..4a2f96055f03a0a250b1183217fc948042c3b472 100644
--- a/sequence_extractor/exon_concatenation.py
+++ b/sequence_extractor/exon_concatenation.py
@@ -1,19 +1,15 @@
-fa = open("fasta.fa",'r')
-lines = fa.readlines()
-for x in range(int(len(lines)/2)):
-    if x == 0:
-        annotation = lines[0]
-        read = lines[1]
-    if x >= 1:
-        if lines[2*x] == lines[2*(x-1)]:
-            read+= lines[(2*x)+1]
-        else:
-            annotation = lines[2*x]
-            read = lines[(2*x)+1]
-
-# Function for random addition of polyA to sequences
-
-import numpy as np
-
-listA = ['A','U','G','C']
-''.join(np.random.choice(listA,250,p=[0.9,0.040,0.020,0.020]))
+def exon_concatenation:
+	fa = open("fasta.fa",'r')
+	lines = fa.readlines()
+	for x in range(int(len(lines)/2)):
+    		if x == 0:
+        		annotation = lines[0]
+        		read = lines[1]
+    		if x >= 1:
+        		if lines[2*x] == lines[2*(x-1)]:
+            			read+= lines[(2*x)+1]
+        		else:
+				return annotation
+				return read
+            			annotation = lines[2*x]
+            			read = lines[(2*x)+1]
diff --git a/sequence_extractor/polyA.py b/sequence_extractor/polyA.py
new file mode 100644
index 0000000000000000000000000000000000000000..2c88e791ec06ef9024ffc45904a6f0c39d2258e8
--- /dev/null
+++ b/sequence_extractor/polyA.py
@@ -0,0 +1,6 @@
+import numpy as np
+
+def PolyA_generator(exon):
+	listA = ['A','T','G','C']
+	polyA = ''.join(np.random.choice(listA,250,p=[0.9,0.040,0.020,0.020]))
+	return (exon+polyA)