待整理

Posted by Icey的博客 on September 19, 2019

https://www.jianshu.com/p/471283080bd6 seqkit seqkit stat seqkit fx2tab -l -g -n -i -H assembly.fasta

按染色体拆分fasta

for i in `cat ref.fa|grep ">"|cut -c 2-`
do
cat ref.fa|seqkit grep -p ${i} > ${i}.fa
done