seqkit split -p 20 ../../test_data/protein.chr4.s.fasta --out-dir protein.split
for i in {01..20};
do
echo "exonerate --model protein2genome --percent 50 -t ../genome.hardmasked \
-q protein.split/*part_0${i}.fasta --showtargetgff yes --showalignment no \
--score 100 --minintron 20 --maxintron 50000 >exonerate${i}.out ">>cmds.sh
done
parallel -j 10 <cmds.sh
parallel -j 10 <w.sh
