1、先用fastaDeal.pl,统计id和len的长度;然后用awk按照300000000,形成切分id;根据这些id用fishInWinter.pl脚本进行调取
路径:/project/nextomics/liujzh/DJX.TMP fastaDeal.pl -attr id:len gy.correctedReads.fasta > gy.correctedReads.fasta.len awk 'BEGIN{i=0;sum=0}{print $1>"p."i".id";sum+=$2}sum>300000000{sum=0;i+=1}' gy.correctedReads.fasta.len ls p.*.id | while read i;do /export/Script/fishInWinter.pl -bf table -ff fasta $i gy.correctedReads.fasta > $i.fasta;done