#!/bin/bash if [ ! -d $1 ]; then echo "$1 is not a directory" exit 1 fi echo "processing $1 ..." rm -f $1/bond_analyze.txt for i in $1/*.save; do bai="`./bond_analyze $i | grep '\[bond' | grep quality`" time=`echo $bai | awk '{ print $4 }'` quality4=`echo $bai | awk '{ print $8 }'` quality3=`echo $bai | awk '{ print $10 }'` count=`echo $bai | awk '{ print $12 }'` echo "$time $quality4 $quality3 $count" >> $1/bond_analyze.txt done echo "done"