pep and cdna of 5 Bacillus subtilis strains were downloaded from Ensembl


formatdb -i Bacillus_subtilis_subsp_subtilis_str_168.GCA_000009045.1.17.pep.all.fa -p T
formatdb -i Bacillus_subtilis_subsp_spizizenii_tu_b_10.GCA_000227465.1.17.pep.all.fa -p T
formatdb -i Bacillus_subtilis_bsn5.GCA_000186745.1.17.pep.all.fa -p T
formatdb -i Bacillus_subtilis_subsp_subtilis_str_ro_nn_1.GCA_000227485.1.17.pep.all.fa -p T
formatdb -i Bacillus_subtilis_subsp_spizizenii_str_w23.GCA_000146565.1.17.pep.all.fa -p T

blastall -p blastp -d Bacillus_subtilis_subsp_subtilis_str_168.GCA_000009045.1.17.pep.all.fa -i Bacillus_subtilis_subsp_spizizenii_tu_b_10.GCA_000227465.1.17.pep.all.fa -o 99.01168_tu.txt -m 8 -e 10
blastall -p blastp -d Bacillus_subtilis_subsp_subtilis_str_168.GCA_000009045.1.17.pep.all.fa -i Bacillus_subtilis_bsn5.GCA_000186745.1.17.pep.all.fa -o 99.01168_bsn5.txt -m 8 -e 10
blastall -p blastp -d Bacillus_subtilis_subsp_subtilis_str_168.GCA_000009045.1.17.pep.all.fa -i Bacillus_subtilis_subsp_subtilis_str_ro_nn_1.GCA_000227485.1.17.pep.all.fa -o 99.01168_ro.txt -m 8 -e 10
blastall -p blastp -d Bacillus_subtilis_subsp_subtilis_str_168.GCA_000009045.1.17.pep.all.fa -i Bacillus_subtilis_subsp_spizizenii_str_w23.GCA_000146565.1.17.pep.all.fa -o 99.01168_w23.txt -m 8 -e 10

blastall -p blastp -d Bacillus_subtilis_subsp_spizizenii_tu_b_10.GCA_000227465.1.17.pep.all.fa -i Bacillus_subtilis_subsp_subtilis_str_168.GCA_000009045.1.17.pep.all.fa -o 99.01tu_168.txt -m 8 -e 10
blastall -p blastp -d Bacillus_subtilis_bsn5.GCA_000186745.1.17.pep.all.fa -i Bacillus_subtilis_subsp_subtilis_str_168.GCA_000009045.1.17.pep.all.fa -o 99.01bsn5_168.txt -m 8 -e 10
blastall -p blastp -d Bacillus_subtilis_subsp_subtilis_str_ro_nn_1.GCA_000227485.1.17.pep.all.fa -i Bacillus_subtilis_subsp_subtilis_str_168.GCA_000009045.1.17.pep.all.fa -o 99.01ro_168.txt -m 8 -e 10
blastall -p blastp -d Bacillus_subtilis_subsp_spizizenii_str_w23.GCA_000146565.1.17.pep.all.fa -i Bacillus_subtilis_subsp_subtilis_str_168.GCA_000009045.1.17.pep.all.fa -o 99.01w23_168.txt -m 8 -e 10

perl 02.orth.pl
