“wie man große Sequenzen unter Linux hackt” Code-Antworten

wie man große Sequenzen unter Linux hackt

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