count hits from FOAM hmmsearch

foamcount.sh

for FILE in *.tsv; do 
  LI=`grep -v '^#' $FILE | wc -l`; 
  SEQS=`gunzip -c ./../NPRB/$F1.faa.gz | grep -c '^>'`; 
  BP=`gunzip -c ./../NPRB/$F1.faa.gz | grep -v '^>' | wc -b`; 
  echo -e "$FILE\t$LI" >> target_hits.txt; 
done

rm target_hits2.txt;

for FILE in *.tsv; do 
  echo $FILE; 
  F1=`basename $FILE .tsv | cut -f2- -d'.'`; 
  LI=`grep -v '^#' $FILE | wc -l`; 
  SEQS=`gunzip -c ./../*/spades-assembly/$F1.faa.gz | grep -c '^>'`; 
  BP=`gunzip -c ./../*/spades-assembly/$F1.faa.gz | grep -v '^>' | wc -m`; 
  echo -e "$FILE\t$LI\t$SEQS\t$BP" >> target_hits2.txt; 
done