readfq 软件的安装 统计 reads的数目及碱基数目

 

01、下载

git clone https://github.com/billzt/readfq.git

 

02、编译

gcc -o kseq_fastq_base kseq_fastq_base.c -lz

 

03、统计fastq文件的reads数目和碱基数目

#!/bin/bash
for i in $(find SRR*)
do
        printf $i"\t" >> result.txt; kseq_fastq_base $i >> result.txt
done

(base) [b20223040323@admin1 test]$ cat result.txt   ## 统计结果
SRR21814498_1.fastq.gz  Num reads:11671642      Num Bases: 1750746300
SRR21814498_2.fastq.gz  Num reads:11671642      Num Bases: 1750746300
SRR21814509_1.fastq.gz  Num reads:11420702      Num Bases: 1713105300
SRR21814509_2.fastq.gz  Num reads:11420702      Num Bases: 1713105300
SRR21814514_1.fastq.gz  Num reads:14051283      Num Bases: 2107692450
SRR21814514_2.fastq.gz  Num reads:14051283      Num Bases: 2107692450