# HG changeset patch # User jshay # Date 1565900339 14400 # Node ID ce6806e0539f1644d50eace99497ff4fd1cbff02 # Parent e40ccd3eab4483e90a9207fce146c28fb9c0b3fa Uploaded diff -r e40ccd3eab44 -r ce6806e0539f filterinfo/filterinfo.sh --- a/filterinfo/filterinfo.sh Tue Aug 13 15:34:03 2019 -0400 +++ b/filterinfo/filterinfo.sh Thu Aug 15 16:18:59 2019 -0400 @@ -37,13 +37,13 @@ divideby=4 fi # print % reads passing filter to an outfile - if [[ $FILTERED == *.gz ]]; then + if [ $(file $FILTERED | awk '{print $2}') == "gzip" ]; then pass=$(gunzip -c $FILTERED | wc -l | awk '{print $1}') else pass=$(wc -l $FILTERED | awk '{print $1}') fi pass=$(expr $pass / $divideby) - if [[ $FULL == *.gz ]]; then + if [ $(file $FULL | awk '{print $2}') == "gzip" ]; then total=$(gunzip -c $FULL | wc -l | awk '{print $1}') else total=$(wc -l $FULL | awk '{print $1}') diff -r e40ccd3eab44 -r ce6806e0539f filterinfo/test --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/filterinfo/test Thu Aug 15 16:18:59 2019 -0400 @@ -0,0 +1,2 @@ +Input Reads Reads Passing Filter Fraction Reads Passing Filter +100000 100000 1.00000