如何写shell脚本统计fasta序列的长度和gc含量

可以使用以下shell脚本来统计fasta序列的长度和GC含量：

#!/bin/bash

# 读取fasta文件
fasta_file=$1
fasta=$(cat $fasta_file)

# 统计序列长度
len=$(echo "$fasta" | grep -v "^>" | tr -d "\n" | wc -m)
echo "Sequence length: $len"

# 统计GC含量
gc=$(echo "$fasta" | grep -v "^>" | tr -d "\n" | grep -o "[GCgc]" | wc -l)
gc_perc=$(echo "scale=2;$gc/$len*100" | bc)
echo "GC content: $gc_perc%"

使用该脚本时，需要将fasta文件作为第一个参数传递给脚本，例如：

bash fasta_stats.sh input.fasta

该脚本首先读取fasta文件，然后使用grep和tr命令来统计序列长度和GC含量。最后，使用echo命令打印结果。

相关问题

提取fasta序列中的指定位置的长度序列的脚本

以下是一个 Python 脚本，可用于从 FASTA 格式的文件中提取指定位置的长度序列：

def read_fasta_file(filename):
    """读取 FASTA 格式的文件并返回一个字典，其中键为序列 ID，值为序列字符串"""
    sequences = {}
    with open(filename, 'r') as f:
        current_id = None
        current_sequence = ''
        for line in f:
            line = line.strip()
            if line.startswith('>'):
                if current_id is not None:
                    sequences[current_id] = current_sequence
                current_id = line[1:]
                current_sequence = ''
            else:
                current_sequence += line
        if current_id is not None:
            sequences[current_id] = current_sequence
    return sequences

def extract_sequence(sequence, start, end):
    """从序列中提取指定位置的长度序列"""
    return sequence[start-1:end]

if __name__ == '__main__':
    # 读取 FASTA 文件
    sequences = read_fasta_file('sequences.fasta')

    # 提取指定位置的长度序列
    sequence_id = 'sequence1'
    start = 10
    end = 20
    extracted_sequence = extract_sequence(sequences[sequence_id], start, end)

    # 输出结果
    print(f'>{sequence_id}_{start}_{end}')
    print(extracted_sequence)

你只需要将上述脚本中的 `sequences.fasta` 文件名、`sequence1` 序列 ID、`start` 和 `end` 位置替换为你自己的需要即可。

写一段计算fasta长度的脚本 Python

脚本如下：f = open("fasta.txt", "r")
seqs = f.readlines()for seq in seqs:
print("长度：",len(seq))f.close()