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
|
// -----------------------------------------------------------------------------------------------------
// Copyright (c) 2006-2020, Knut Reinert & Freie Universität Berlin
// Copyright (c) 2016-2020, Knut Reinert & MPI für molekulare Genetik
// This file may be used, modified and/or redistributed under the terms of the 3-clause BSD-License
// shipped with this file and also available at: https://github.com/seqan/seqan3/blob/master/LICENSE.md
// -----------------------------------------------------------------------------------------------------
#include <gtest/gtest.h>
#include <utility>
#include <vector>
#include <seqan3/alignment/aligned_sequence/debug_stream_alignment.hpp>
#include <seqan3/alignment/pairwise/alignment_result.hpp>
#include <seqan3/alphabet/gap/gapped.hpp>
#include <seqan3/alphabet/nucleotide/dna4.hpp>
#include <seqan3/core/detail/debug_stream_tuple.hpp>
TEST(alignment_result_test, debug_streamable)
{
using coordinate_t = std::pair<int, int>;
using aligned_seq_type = std::vector<seqan3::gapped<seqan3::dna4>>;
using seqan3::operator""_dna4;
size_t id = 3;
int score = -15;
coordinate_t begin_coordinate{4, 6};
coordinate_t end_coordinate{23, 35};
aligned_seq_type gapped_seq1{'A'_dna4, 'T'_dna4, seqan3::gap{}, 'A'_dna4};
aligned_seq_type gapped_seq2{'A'_dna4, 'T'_dna4, 'C'_dna4, seqan3::gap{}};
std::pair<aligned_seq_type, aligned_seq_type> alignment{gapped_seq1, gapped_seq2};
std::ostringstream ostream{};
seqan3::debug_stream_type debug_stream{ostream};
{ // Print id and score
seqan3::detail::alignment_result_value_type result_value{id, id, score};
seqan3::alignment_result result{result_value};
debug_stream << result;
EXPECT_EQ(ostream.str(), "{sequence1 id: 3, sequence2 id: 3, score: -15}");
}
{ // Print id and score and back coordinate
ostream.str("");
seqan3::detail::alignment_result_value_type result_value{id, id, score, end_coordinate};
seqan3::alignment_result result{result_value};
EXPECT_TRUE((seqan3::detail::is_type_specialisation_of_v<decltype(result), seqan3::alignment_result>));
debug_stream << result;
EXPECT_EQ(ostream.str(), "{sequence1 id: 3, sequence2 id: 3, score: -15, end: (23,35)}");
}
{ // Print id and score and back and front coordinate
ostream.str("");
seqan3::detail::alignment_result_value_type result_value{id, id, score, end_coordinate, begin_coordinate};
seqan3::alignment_result result{result_value};
debug_stream << result;
EXPECT_EQ(ostream.str(), "{sequence1 id: 3, sequence2 id: 3, score: -15, begin: (4,6), end: (23,35)}");
}
{ // Print id and score and back and front coordinate and alignment
ostream.str("");
seqan3::detail::alignment_result_value_type result_value{id,
id,
score,
end_coordinate,
begin_coordinate,
alignment};
seqan3::alignment_result result{result_value};
debug_stream << result;
EXPECT_EQ(ostream.str(), "{sequence1 id: 3, sequence2 id: 3, score: -15, begin: (4,6), end: (23,35), \n"
"alignment:\n"
" 0 \n"
" AT-A\n"
" || \n"
" ATC-\n"
"}");
}
}
|
