def parse_genes(prod_genes):
    """
    Parse the the Biopython SeqRecord representation of the Prodigal GFF output. Specifically, retrieve the
    start positions, end positions, strand orientation and sequences of each gene
    :param prod_genes: A SeqRecord object with genes found by Prodigal
    :return: A dict with the retrieved data
    """
    starts = []
    ends = []
    strands = []
    seqs = []
    for record in prod_genes:
        data = prod_genes[record].description.split(' # ')
        starts.append(int(data[1]))
        ends.append(int(data[2]))
        strands.append(int(data[3]))
        seqs.append(prod_genes[record].seq)
    samp_id = data[0]
    return {'start': starts, 'end': ends, 'strand': strands, 'sequence':
        seqs, 'id': samp_id}
