Mercurial > repos > pjbriggs > pal_finder

annotate pal_finder_wrapper_utils.sh @ 15:a3af1ff4cad1 draft

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
pal_finder 0.02.04.7 for testing.
author pjbriggs
date Mon, 14 May 2018 11:10:19 -0400
parents e95e2b7b3c84
children
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
15
a3af1ff4cad1 pal_finder 0.02.04.7 for testing.
pjbriggs
parents: 10
diff changeset
41 local outfile=$2
10
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
42 local pattern="^(SEQUENCE_ID|PRIMER_PRODUCT_SIZE_RANGE)"
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
43 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
44 do
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
45 # 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
46 # keywords in the primer3 input
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
47 if [ ! -z "$(echo $line | grep ^SEQUENCE_ID)" ] ; then
15
a3af1ff4cad1 pal_finder 0.02.04.7 for testing.
pjbriggs
parents: 10
diff changeset
48 # Lines look like:
a3af1ff4cad1 pal_finder 0.02.04.7 for testing.
pjbriggs
parents: 10
diff changeset
49 # SEQUENCE_ID=(AT_1_16)(AT_1_16)M00879:99:000000000-AH9KG:1:2107:10006:2535*PRIMER_PRODUCT_SIZE_RANGE=194-329^59-194
10
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
50 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
51 local seq_id=$(echo $line | cut -d'*' -f1 | cut -d'=' -f2)
15
a3af1ff4cad1 pal_finder 0.02.04.7 for testing.
pjbriggs
parents: 10
diff changeset
52 elif [ ! -z "$(echo $line | grep ^PRIMER_PRODUCT_SIZE_RANGE)" ] ; then
a3af1ff4cad1 pal_finder 0.02.04.7 for testing.
pjbriggs
parents: 10
diff changeset
53 # Lines look like:
a3af1ff4cad1 pal_finder 0.02.04.7 for testing.
pjbriggs
parents: 10
diff changeset
54 # PRIMER_PRODUCT_SIZE_RANGE=194-329^59-194*SEQUENCE_ID=(AT_1_16)(AT_1_16)M00879:99:000000000-AH9KG:1:2107:10006:2535
a3af1ff4cad1 pal_finder 0.02.04.7 for testing.
pjbriggs
parents: 10
diff changeset
55 local size_range=$(echo $line | cut -d'*' -f1 | cut -d'=' -f2 | tr '^' ' ')
10
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
56 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
57 fi
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
58 seq_id=$(echo $seq_id | cut -d')' -f3)
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
59 # 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
60 # to see if it's okay
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
61 local bad_range=
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
62 for range in $(echo $size_range) ; do
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
63 local lower=$(echo $range | cut -d'-' -f1)
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
64 local upper=$(echo $range | cut -d'-' -f2)
15
a3af1ff4cad1 pal_finder 0.02.04.7 for testing.
pjbriggs
parents: 10
diff changeset
65 if [ "$lower" -gt "$upper" ] ; then
10
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
66 bad_range=yes
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
67 break
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
68 fi
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
69 done
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
70 # Report if the range is wrong
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
71 if [ ! -z "$bad_range" ] ; then
15
a3af1ff4cad1 pal_finder 0.02.04.7 for testing.
pjbriggs
parents: 10
diff changeset
72 echo "${seq_id}"$'\t'"(${size_range})" >>$outfile
10
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
73 fi
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
74 done
e95e2b7b3c84 Fix previous upload (was missing pal_finder_wrapper_utils.sh file)
pjbriggs
parents:
diff changeset
75 }