perl -e '$file = "seq001"; for (@ARGV) { if (10 <= $_ && $_ <= 40) { read STDIN, my $seq, $_; open my $FH, ">", $file++; print $FH $seq } }' 10 18 28 < fasta.txt