3
# Script for checking output of executable src/long_seq_tm_test
4
# To make long_seq_tm_test go to ../src, and do
5
# make long_seq_tm_test
7
my ($seq, $start, $len) = @ARGV;
8
my $s = substr($seq, $start, $len);
9
my $GC_count = $s =~ tr/gGcC/gGcC/;
11
$r = log(50/1000.0)/log(10);
12
$tm = 81.5 + 16.6*$r + 41*$GC_count/$len - 600./$len;