sh 单域与多域蛋白的数量

Posted 2021-05-06

#!/usr/bin/env bash

MULTI=$(awk '{print $3}' $1 | uniq -c | awk '$1 != 1 {print $2}' | wc -l)

SINGLE=$(awk '{print $3}' $1 | uniq -c | awk '$1 == 1 {print $2}' | wc -l)

echo "No. of single domain proteins = $SINGLE"

echo "No. of multidomain proteins = $MULTI"