for i in `ls ./raw_data/*/*.fq.gz`
do
sample_name=`echo $i | awk -F "/" '{print $NF}' | awk -F "." '{print $1}'`
sample=`echo $i | awk -F "/" '{print $NF}'`
cd ${path}
mkdir ${sample_name}
cd ${sample_name}
ln -s $i ./
zcat ${sample} | split -16000000
count=0
for z in `ls x*`
do
count=$((${count}+1))
/bin/mv ${z} ${sample_name}_00${count}.fastq
done
done
创建目录
链接
解压分割
改名
网友评论