统计fastq文件中读段的数量

mycount=`cat SRR108114_new_1.fastq | wc -l`
echo 'Number of clean reads, SRR108114_new_1.fastq: '$((mycount/4)) >> SRR108114.stat
mycount=`cat SRR108114_new_2.fastq | wc -l`
echo 'Number of clean reads, SRR108114_new_2.fastq: '$((mycount/4)) >> SRR108114.stat

posted @ 2015-09-21 18:36  emanlee  阅读(998)  评论(0编辑  收藏  举报