In MATLAB, write a function seq_assemble_greedy( fastqfile ) that takes in a fastqfile containing short reads, implements the greedy algorithm of genome assembly, and returns the assembled contigs as a cell array of texts (in python, return as a list of texts) merge these two reads into a single longer read. Repeat the previous step so long as at least two reads can be found sharing one or more residues overlap. Return the collection of extended merged reads you end up with.

Write a script seq_assemble_greedy_test that contains two tests for your function on two separate fastqfiles you created and prints the result of each test case.

short read 1 = ['DGH','GHIABCD','ABCDEFGHI','FGHIABCDE']

short read 2 = ['ABA','BBA','AAB','BAA','AAA']