From 38af80acbbf2187e802d298b113e1f2c37865a75 Mon Sep 17 00:00:00 2001 From: clara <clara@dyn-35-97.mobile.unibas.ch> Date: Thu, 10 Nov 2022 14:53:07 +0100 Subject: [PATCH] new file: read_sequencer_package/generate_sequences.py --- read_sequencer_package/generate_sequences.py | 25 ++++++++++++++++++++ 1 file changed, 25 insertions(+) create mode 100644 read_sequencer_package/generate_sequences.py diff --git a/read_sequencer_package/generate_sequences.py b/read_sequencer_package/generate_sequences.py new file mode 100644 index 0000000..aa9932f --- /dev/null +++ b/read_sequencer_package/generate_sequences.py @@ -0,0 +1,25 @@ + +import random + +def generate_sequences(n, mean, sd): + """Summary line. + + Generates random sequences. + + Args: + n (int): Amount of sequences to generate. + mean (int): mean length of sequence (gaussian distribution). + sd (float): standart deviation of length of sequence (gaussian distribution). + + Returns: + list: of n sequences + """ + l1 = [] + for i in range(n): + seq = "" + nt = ["A", "T", "C", "G"] + for pos in range(round(random.gauss(mean, sd))): + seq = seq + random.choice(nt) + l1.append(seq) + return l1 + -- GitLab