Mercurial > repos > pjbriggs > pal_finder

annotate pal_finder_wrapper_utils.sh @ 11:84613d5f1713 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Uploaded fix to printing read IDs with associated with 'bad' ranges.
author pjbriggs
date Mon, 24 Jul 2017 11:04:27 -0400
parents e95e2b7b3c84
children a3af1ff4cad1
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
10
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
1 #!/bin/bash
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
2 #
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
3 # Helper functions for the pal_finder_wrapper.sh script
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
4 #
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
5 # Utility function for terminating on fatal error
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
6 function fatal() {
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
7 echo "FATAL $@" >&2
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
8 exit 1
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
9 }
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
10 #
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
11 # Check that specified program is available
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
12 function have_program() {
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
13 local program=$1
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
14 local got_program=$(which $program 2>&1 | grep "no $(basename $program) in")
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
15 if [ -z "$got_program" ] ; then
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
16 echo yes
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
17 else
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
18 echo no
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
19 fi
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
20 }
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
21 #
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
22 # Set the value for a parameter in the pal_finder config file
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
23 function set_config_value() {
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
24 local key=$1
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
25 local value=$2
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
26 local config_txt=$3
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
27 if [ -z "$value" ] ; then
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
28 echo "No value for $key, left as default"
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
29 else
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
30 echo Setting "$key" to "$value"
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
31 sed -i 's,^'"$key"' .*,'"$key"' '"$value"',' $config_txt
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
32 fi
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
33 }
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
34 #
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
35 # Identify 'bad' PRIMER_PRODUCT_SIZE_RANGE from pr3in.txt file
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
36 function find_bad_primer_ranges() {
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
37 # Parses a pr3in.txt file from pal_finder and reports
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
38 # sequence ids where the PRIMER_PRODUCT_SIZE_RANGE has
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
39 # upper limit which is smaller than lower limit
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
40 local pr3in=$1
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
41 local pattern="^(SEQUENCE_ID|PRIMER_PRODUCT_SIZE_RANGE)"
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
42 for line in $(grep -E "$pattern" $pr3in | sed 's/ /^/' | sed 'N;s/\n/*/')
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
43 do
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
44 # Loop over pairs of SEQUENCE_ID and PRIMER_PRODUCT_SIZE_RANGE
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
45 # keywords in the primer3 input
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
46 if [ ! -z "$(echo $line | grep ^SEQUENCE_ID)" ] ; then
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
47 # Extract the values
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
48 local size_range=$(echo $line | cut -d'*' -f2 | cut -d'=' -f2 | tr '^' ' ')
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
49 local seq_id=$(echo $line | cut -d'*' -f1 | cut -d'=' -f2)
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
50 else
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
51 local size_range=$(echo $line | cut -d'*' -f1 | cut -d'=' -f2)
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
52 local seq_id=$(echo $line | cut -d'*' -f2 | cut -d'=' -f2)
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
53 fi
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
54 seq_id=$(echo $seq_id | cut -d')' -f3)
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
55 # Check the upper and lower limits in each range
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
56 # to see if it's okay
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
57 local bad_range=
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
58 for range in $(echo $size_range) ; do
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
59 local lower=$(echo $range | cut -d'-' -f1)
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
60 local upper=$(echo $range | cut -d'-' -f2)
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
61 if [ $lower -gt $upper ] ; then
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
62 bad_range=yes
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
63 break
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
64 fi
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
65 done
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
66 # Report if the range is wrong
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
67 if [ ! -z "$bad_range" ] ; then
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
68 echo "$seq_id ($size_range)"
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
69 fi
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
70 done
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
71 }