參考一
參考二

測序reads數(shù)兩個例子

#singel end
folder_path="/path_to/donor2"
output_file="/path_to/donor2_nanopore_coverage_2.txt"

# 遍歷文件夾中的所有fastq文件
for fastq_file in ${folder_path}/*.fastq; do
    # 使用SAMtools命令計算reads數(shù)量
    read_count=$(expr $(cat ${fastq_file} | wc -l) / 4)
    # 打印每個文件的reads數(shù)量
    echo "${fastq_file}: ${read_count} reads" >>$output_file
done

#pair end
folder_path="/path_to/short_reads/donor1"
output_file="/path_to/short_reads/donor1_short_reads_coverage.txt"

for i in 1 2 3 4 5 6 7 8
do
    read_count1=$(zcat "${folder_path}/s${i}_R1_val_1.fq.gz" | grep -c '^+$')
    read_count2=$(zcat "${folder_path}/s${i}_R2_val_2.fq.gz" | grep -c '^+$')
    read_count=$(expr($read_count1+$read_count2))
    echo "s${i}: ${read_count} reads" >>$output_file
done