“comment hacher de grandes séquences dans Linux” Réponses codées

comment hacher de grandes séquences dans 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