blob: 1f495eb9b984762768e3b029555ea6a9240d62f8 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
|
#!/usr/bin/env bash
if [ -z "$foreground" ] && [ -z "$background" ]; then
if (( $# < 2 )); then
echo "please supply two genomes, foreground and background"
exit 1
fi;
: ${foreground=$1}
: ${background=$2}
fi
if [[ ! -f $foreground ]]; then
echo "Could not open $foreground"
exit 1
fi
if [[ ! -f $background ]]; then
echo "Could not open $background"
exit 1
fi
# output directory
: ${output_directory=`basename $foreground`_`basename $background`}
# temp directory
: ${tmp_directory=$output_directory/.tmp}
# directory to store our counts and sorted counts
: ${counts_directory=$tmp_directory}
# range of mers, min and max
: ${min_mer_range=6}
: ${max_mer_range=10}
# max mer distance, the distance between two mers in our selected outputs
: ${max_mer_distance=5000}
# min/maximum kmer meling point
: ${max_melting_temp=30}
: ${min_melting_temp=0}
# minimum mer count
: ${min_mer_count=0}
# maximum mers to pick
: ${max_select=15}
# mers to specifically IGNORE, space delimited
: ${ignore_mers=''}
export ignore_mers
export min_mer_range
export max_mer_range
export max_select
export min_mer_count
export max_mer_distance
export max_melting_temp
export min_melting_temp
# Make our output directory
if [ ! -d $output_directory ]; then
mkdir $output_directory
fi
# Make our counts directory
if [ ! -d $counts_directory ]; then
mkdir $counts_directory
fi
# Make our temporary directory
if [ ! -d $tmp_directory ]; then
mkdir $tmp_directory
fi
for fasta_file in $foreground $background; do
counts=$counts_directory/$(basename $fasta_file)
tmp=$tmp_directory/$(basename $fasta_file)
echo pre-processing $fasta_file
# check if our preprocessed file exists
if [[ ! -f $tmp ]]; then
echo "> pre processed $fasta_file" >> $tmp
cat $fasta_file | grep -v "^>" | tr -d '\n' >> $tmp
fi
# run counts if they haven't been created
rm $counts-counts
for mer in `seq $min_mer_range $max_mer_range`; do
if [ ! -e $counts-counts-$mer ]; then
echo checking $mer mers for $fasta_file
kmer_total_count -i $tmp -k $mer -l -n >> $counts-counts-$mer
else
echo "$mer mers already done for $fasta_file"
fi
cat $counts-counts-$mer >> $counts-counts
done
done
fg_counts=$counts_directory/$(basename $foreground)-counts
bg_counts=$counts_directory/$(basename $background)-counts
fg_tmp=$tmp_directory/$(basename $foreground)
bg_tmp=$tmp_directory/$(basename $background)
# remove ignored mers
if [ "$ignore_mers" ]; then
echo "removing ignored mers: " + $ignore_mers
for mer in $ignore_mers; do
sed -i '/^'$mer'\t/d' $fg_counts
sed -i '/^'$mer'\t/d' $bg_counts
done
fi
echo "checking if mers are below melting temperature in the foreground"
rm $fg_counts-fg-non-melting
melting_range $min_melting_temp $max_melting_temp < $fg_counts > $fg_counts-fg-non-melting &
echo "checking if mers are below melting temperature in the background"
rm $bg_counts-bg-non-melting
melting_range $min_melting_temp $max_melting_temp < $bg_counts > $bg_counts-bg-non-melting
# echo "scoring mer selectivity"
# python ./mer_selectivity.py $fg_counts-fg-non-melting $bg_counts-bg-non-melting
echo ""
echo "scoring mers"
select_mers.py $fg_counts-fg-non-melting $fg_tmp $bg_counts-bg-non-melting $bg_tmp $output_directory/output_`date +%s`
|
