切割

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

猜你喜欢

转载自www.cnblogs.com/djx571/p/10298243.html
今日推荐