“Cómo cortar secuencias grandes en Linux” Código de respuesta

Cómo cortar secuencias grandes en Linux

awk 'BEGIN {n_seq=0;} /^>/ {if(n_seq%1000==0){file=sprintf("myseq%d.fa",n_seq);} print >> file; n_seq++; next;} { print >> file; }' < sequences.fa