||
#! /bin/bash #download https://github.com/billzt/readfq #purpose: calculate reads from fastq files #set env readfq="/home/zhanghl/supporting_softwares/readfq/readfq-master/kseq_fastq_base" fq_dir="" cd $fq_dir #fq所在位置 ls *.fq.gz > fq.list #collect fastq files for i in $(cat fq.list) do printf $i "\t" >> fq.reads_num $readfq $i >> fq.reads_num done
Archiver|手机版|科学网 ( 京ICP备07017567号-12 )
GMT+8, 2024-11-9 07:18
Powered by ScienceNet.cn
Copyright © 2007- 中国科学报社