#
grep -i '>' yb_4001-8000.orf.fasta > yb_4001-8000.orfs
echo '
########
 ' >> yb_4001-8000.orfs
grep -i '>' yb_4001-8000.orf.fasta | sed s/' '/" stdout | blastall -p blastp -d entero -e 0.01 | grep -h 'Number of sequences better than 1.0e-02:' >> yb_4001-8000.orfs # "/ | sed s/'>'/'seqret yb_4001-8000.orf.fasta:'/ > 1  
chmod +x 1
./1
rm 1