Скрипт для получения процентов идентичности генов:

needle g1.fasta g2.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" > percent_id_gene
needle g1.fasta g3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g1.fasta g4.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g1.fasta g5.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g1.fasta g6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g1.fasta g7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
#
needle g1.fasta g3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_gene
#
needle g2.fasta g3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g2.fasta g4.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g2.fasta g5.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g2.fasta g6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g2.fasta g7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
#
needle g1.fasta g3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_gene
#
needle g3.fasta g4.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g3.fasta g5.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g3.fasta g6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g3.fasta g7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
#
needle g1.fasta g3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_gene
#
needle g4.fasta g5.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g4.fasta g6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g4.fasta g7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
#
needle g1.fasta g3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_gene
#
needle g5.fasta g6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
needle g5.fasta g7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene
#
needle g1.fasta g3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_gene
#
needle g6.fasta g7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_gene

Скрипт для получения процентов идентичности белков:

needle p1.fasta p2.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" > percent_id_prot
needle p1.fasta p3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p1.fasta p4.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p1.fasta p5.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p1.fasta p6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p1.fasta p7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
#
needle p1.fasta p3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_prot
#
needle p2.fasta p3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p2.fasta p4.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p2.fasta p5.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p2.fasta p6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p2.fasta p7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
#
needle p1.fasta p3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_prot
#
needle p3.fasta p4.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p3.fasta p5.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p3.fasta p6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p3.fasta p7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
#
needle p1.fasta p3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_prot
#
needle p4.fasta p5.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p4.fasta p6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p4.fasta p7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
#
needle p1.fasta p3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_prot
#
needle p5.fasta p6.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
needle p5.fasta p7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot
#
needle p1.fasta p3.fasta -gapopen 10 -gapextend 0.5 stdout | grep "####" >> percent_id_prot
#
needle p6.fasta p7.fasta -gapopen 10 -gapextend 0.5 stdout | grep "Identity" >> percent_id_prot