from Bio import SeqIO

def print_sequence_lengths_and_gc(fasta_file):
    """Print the ID, length, and GC content of each sequence in a FASTA file."""
    for record in SeqIO.parse(fasta_file, "fasta"):
        seq = record.seq
        gc_count = seq.count('G') + seq.count('C')
        gc_content = (gc_count / len(seq)) * 100
        print(f"Sequence ID: {record.id}\tLength: {len(seq)}\tGC Content: {gc_content:.2f}%")

# Usage
fasta_file = "assembly.fasta"  # Replace with your file
print_sequence_lengths_and_gc(fasta_file)

