sequences = {}
current_gene_id = None
current_sequence = ""

# Use context manager to safely open and read the file
with open(filename, 'r') as file:
    for line in file:
        line = line.strip()  # Remove whitespace

        if line.startswith('>'):
            # Save previous sequence if we have one
            if current_gene_id is not None:
                sequences[current_gene_id]['sequence'] = current_sequence

            # Parse new header
            header_parts = line[1:].split(' ', 1)  # Split on first space only
            current_gene_id = header_parts[0]

            # Store header info
            sequences[current_gene_id] = {
                'header': line[1:],  # Full header without >
                'sequence': ""
            }

            current_sequence = ""  # Reset sequence
            print(f"Found sequence: {current_gene_id}")

        else:
            # Add to current sequence (sequences can span multiple lines)
            current_sequence += line.upper()

# Don't forget the last sequence!
if current_gene_id is not None:
    sequences[current_gene_id]['sequence'] = current_sequence

# Display results
for gene_id, data in sequences.items():
    print(f"\nGene: {gene_id}")
    print(f"Header: {data['header']}")
    print(f"Length: {len(data['sequence'])} bases")
    print(f"First 50 bases: {data['sequence'][:50]}...")