# 读取fasta文件并删除空白行
def read_fasta(file_path):
    with open(file_path, 'r') as f:
        lines = f.readlines()
    
    cleaned_lines = []
    sequence = ""
    for line in lines:
        line = line.strip()  # 去除行末的换行符和空格
        if line:
            if line.startswith('>'):
                cleaned_lines.append(sequence)
                cleaned_lines.append(line)
                sequence = ""
            else:
                sequence += line
    cleaned_lines.append(sequence)

    return cleaned_lines

分割fasta文件成15份

def split_fasta(file_path, num_parts):
    fasta_lines = read_fasta(file_path)
    total_lines = len(fasta_lines)
    lines_per_part = total_lines // num_parts

    for i in range(num_parts):
        start_idx = i * lines_per_part
        end_idx = start_idx + lines_per_part
        part_lines = fasta_lines[start_idx:end_idx]

        # 创建一个新的FASTA部分文件
        with open(f'part_{i + 1}.fasta', 'w') as f:
            f.write('\n'.join(part_lines))

if __name__ == "__main__":
    input_file = '/Users/yaozhaoqun/Desktop/original.fasta'  # 替换为你的fasta文件路径
    num_parts = 15  # 分割成15份

    # 调用split_fasta函数来分割FASTA文件
    split_fasta(input_file, num_parts)