Prot

#! usr/bin/bash

ALL=$(zgrep "^PE" $1 | wc -l)
TR=$(zgrep "KW" $1 | grep "Transmembrane" | wc -l)
EN=$(zgrep "DE" $1 | grep 'EC=' | wc -l)
VI=$(zgrep "KW" $1 | grep "Virulence" | wc -l)
TRp=$(bc<<<"scale=2;$TR*100/$ALL")
ENp=$(bc<<<"scale=2;$EN*100/$ALL")
VIp=$(bc<<<"scale=2;$VI*100/$ALL")


echo "Доля трансмембранных белков в протеоме равна $TRp%"
echo "Доля ферментов в протеоме равна $ENp%"
echo "Доля белков, регулирующих вирулентность, в протеоме $VIp%"