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
|
#include "btllib/util.hpp"
#include "helpers.hpp"
#include <iostream>
int
main()
{
std::string teststring(" AC TG ");
btllib::trim(teststring);
TEST_ASSERT_EQ(teststring, "AC TG");
btllib::CString testcstring(" TG AC ");
btllib::trim(testcstring);
TEST_ASSERT_EQ(std::string(testcstring), "TG AC");
auto actg_split = btllib::split("A|C|T|G", "|");
TEST_ASSERT_EQ(actg_split.size(), 4);
TEST_ASSERT_EQ(actg_split[0], "A");
TEST_ASSERT_EQ(actg_split[1], "C");
TEST_ASSERT_EQ(actg_split[2], "T");
TEST_ASSERT_EQ(actg_split[3], "G");
auto actg_join = btllib::join(actg_split, "/");
TEST_ASSERT_EQ(actg_join, "A/C/T/G");
TEST_ASSERT(btllib::startswith(actg_join, "A/C"));
TEST_ASSERT(btllib::endswith(actg_join, "/T/G"));
std::string qual = "$$%%)*0)'%%&$$%&$&'''*)(((((()55561--.12356577-++**++,////.*))((()+))**010/..--+**++*+++)++++78883";
double avg = btllib::calc_phred_avg(qual, 0, 10);
double avg1 = btllib::calc_phred_avg(qual);
double avg2 = btllib::calc_phred_avg(qual, 0, 4);
double avg3 = btllib::calc_phred_avg(qual, 5, 20);
TEST_ASSERT_LT(std::abs(avg - 5.34264), 1e-4);
TEST_ASSERT_LT(std::abs(avg1 - 8.54327), 1e-4);
TEST_ASSERT_LT(std::abs(avg2 - 3.47128), 1e-4);
TEST_ASSERT_LT(std::abs(avg3 - 5.48923), 1e-4);
return 0;
}
|
