计算基因组GC含量（每10kb）

from Bio import SeqIO

def calculate_gc_content_from_seq(seq):
    gc_count = seq.count('G') + seq.count('C')
    return gc_count / len(seq)

def calculate_gc_content(filename, window_size=10000):
    with open(filename, "r") as handle:
        for record in SeqIO.parse(handle, "fasta"):
            chr_name = record.id
            sequence = str(record.seq)
            for start in range(0, len(sequence), window_size):
                end = min(start + window_size, len(sequence))
                subseq = sequence[start:end]
                gc_content = calculate_gc_content_from_seq(subseq)
                print(chr_name, start, end, gc_content)

filename = "C:/Users/fordata/Downloads/Circos/zy130.fasta"
calculate_gc_content(filename)