awk '/^>/ {printf("\n%s\t",$0);next;} {printf("%s",$0);} END {printf("\n");}' < assembly.fasta | egrep -v '^$' | tr "\t" "\n" | awk '!/^>/ { next } { getline seq } length(seq) > 200 { print $0 "\n" seq }' > assembly.fas
转载于:https://www.cnblogs.com/liuhui0622/p/4616430.html