rust-spoa 0.2.4

/*!
 * @file spoa_test.cpp
 *
 * @brief Spoa unit test source file
 */

#include "spoa_test_config.h"

#include "sequence.hpp"

#include "spoa/spoa.hpp"
#include "bioparser/bioparser.hpp"
#include "gtest/gtest.h"

class SpoaAlignmentTest: public ::testing::Test {
public:
    void SetUp(const std::string& file_name, spoa::AlignmentType alignment_type,
        int8_t match, int8_t mismatch, int8_t gap_open, int8_t gap_extend) {

        parser = bioparser::createParser<bioparser::FastqParser, spoa::Sequence>(
            file_name);
        alignment_engine = spoa::createAlignmentEngine(alignment_type, match,
            mismatch, gap_open, gap_extend);
        graph = spoa::createGraph();
    }

    void TearDown() {}

    void initialize() {
        parser->parse_objects(sequences, -1);

        size_t max_sequence_size = 0;
        for (const auto& it: sequences) {
            max_sequence_size = std::max(max_sequence_size, it->data().size());
        }
        alignment_engine->prealloc(max_sequence_size, 4);
    }

    void construct_partial_order_graph(bool use_qualities) {
        for (const auto& it: sequences) {
            auto alignment = (*alignment_engine)(it->data(), graph);

            if (use_qualities) {
                graph->add_alignment(alignment, it->data(), it->quality());
            } else {
                graph->add_alignment(alignment, it->data());
            }
        }
    }

    std::vector<std::unique_ptr<spoa::Sequence>> sequences;
    std::unique_ptr<bioparser::Parser<spoa::Sequence>> parser;
    std::unique_ptr<spoa::AlignmentEngine> alignment_engine;
    std::unique_ptr<spoa::Graph> graph;
};

TEST(SpoaTest, AlignmentTypeError) {
    try {
        auto alignment_engine = spoa::createAlignmentEngine(
            static_cast<spoa::AlignmentType>(4), 1, 1, -1, -1);
    } catch(std::invalid_argument& exception) {
        EXPECT_STREQ(exception.what(), "[spoa::createAlignmentEngine] error: "
            "invalid alignment type!");
    }
}

TEST(SpoaTest, EmptyInputError) {
    auto alignment_engine = spoa::createAlignmentEngine(
        static_cast<spoa::AlignmentType>(0), 1, 1, -1, -1);
    auto graph = spoa::createGraph();
    auto alignment = (*alignment_engine)("", graph);

    EXPECT_TRUE(alignment.empty());
}

TEST_F(SpoaAlignmentTest, GraphClear) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kSW,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(false);

    auto consensus = graph->generate_consensus();

    graph->clear();
    construct_partial_order_graph(false);

    auto test = graph->generate_consensus();

    EXPECT_TRUE(consensus.compare(test) == 0);
}

TEST_F(SpoaAlignmentTest, LocalConsensus) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kSW,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(false);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "AATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCGA"
        "CCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGG"
        "GAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGC"
        "AGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGT"
        "ACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGC"
        "ACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTT"
        "GAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATA"
        "CGCTTTACACGCGCAACCAAGGATTTCGG";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, LocalAffineConsensus) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kSW,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(false);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "AATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCGA"
        "CCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGG"
        "GAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGC"
        "AGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGT"
        "ACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGC"
        "ACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTT"
        "GAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATA"
        "CGCTTTACACGCGCAACCAAGGATTTCGG";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, LocalConsensusWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kSW,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(true);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "AATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCGA"
        "CCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGG"
        "GAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGC"
        "AGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGT"
        "ACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGC"
        "ACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTT"
        "GAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATA"
        "CGCTTTACACGCGCAACCAAGGATTTCGG";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, LocalAffineConsensusWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kSW,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(true);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "AATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCGA"
        "CCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGG"
        "GAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGC"
        "AGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGT"
        "ACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGC"
        "ACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTT"
        "GAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATA"
        "CGCTTTACACGCGCAACCAAGGATTTCGG";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, GlobalConsensus) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kNW,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(false);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "ATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCGAC"
        "CTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGG"
        "AGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCA"
        "GGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTA"
        "CTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCA"
        "CAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTG"
        "AGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATAC"
        "GC";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, GlobalAffineConsensus) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kNW,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(false);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "ATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCGAC"
        "CTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGG"
        "AGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCA"
        "GGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTA"
        "CTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCA"
        "CAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTG"
        "AGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATAC"
        "GC";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, GlobalConsensusWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kNW,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(true);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "ATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCGAC"
        "CTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGG"
        "AGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCA"
        "GGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTA"
        "CTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCA"
        "CAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTG"
        "AGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATAC"
        "GC";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, GlobalAffineConsensusWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kNW,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(true);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "ATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCGAC"
        "CTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAGGG"
        "AGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGGCA"
        "GGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCGTA"
        "CTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCGCA"
        "CAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGTTG"
        "AGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGATAC"
        "GC";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}


TEST_F(SpoaAlignmentTest, SemiGlobalConsensus) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kOV,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(false);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "ACATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCG"
        "ACCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAG"
        "GGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGG"
        "CAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCG"
        "TACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCG"
        "CACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGT"
        "TGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGAT"
        "ACGCGTTTTACACGCGCAACCAAGGATTTCGG";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, SemiGlobalAffineConsensus) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kOV,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(false);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "GTATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCG"
        "ACCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAG"
        "GGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGG"
        "CAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCG"
        "TACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCG"
        "CACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGT"
        "TGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGAT"
        "ACGCGTTTTACACGCGCAACCAAGGATTTCGG";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, SemiGlobalConsensusWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kOV,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(true);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "ACATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCG"
        "ACCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAG"
        "GGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGG"
        "CAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCG"
        "TACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCG"
        "CACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGT"
        "TGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGAT"
        "ACGCGTTTTACACGCGCAACCAAGGATTTCGG";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, SemiGlobalAffineConsensusWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kOV,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(true);

    auto consensus = graph->generate_consensus();

    std::string valid_result = "ACATGATGCGCTTTGTTGGCGCGGTGGCTTGATGCAGGGGCTAATCG"
        "ACCTCTGGCAACCACTTTTCCATGACAGGAGTTGAATATGGCATTCAGTAATCCCTTCGATGATCCGCAG"
        "GGAGCGTTTTACATATTGCGCAATGCGCAGGGGCAATTCAGTCTGTGGCCGCAACAATGCGTCTTACCGG"
        "CAGGCTGGGACATTGTGTGTCAGCCGCAGTCACAGGCGTCCTGCCAGCAGTGGCTGGAAGCCCACTGGCG"
        "TACTCTGACACCGACGAATTTTACCCAGTTGCAGGAGGCACAATGAGCCAGCATTTACCTTTGGTCGCCG"
        "CACAGCCCGGCATCTGGATGGCAGAAAAACTGTCAGAATTACCCTCCGCCTGGAGCGTGGCGCATTACGT"
        "TGAGTTAACCGGAGAGGTTGATTCGCCATTACTGGCCCGCGCGGTGGTTGCCGGACTAGCGCAAGCAGAT"
        "ACGCGTTTTACACGCGCAACCAAGGATTTCGG";

    EXPECT_TRUE(consensus.compare(valid_result) == 0);
}

TEST_F(SpoaAlignmentTest, LocalMSA) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kSW,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(false);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, LocalAffineMSA) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kSW,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(false);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, LocalMSAWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kSW,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(true);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, LocalAffineMSAWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kSW,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(true);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, GlobalMSA) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kNW,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(false);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, GlobalAffineMSA) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kNW,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(false);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, GlobalMSAWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kNW,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(true);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, GlobalAffineMSAWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kNW,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(true);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, SemiGlobalMSA) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kOV,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(false);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, SemiGlobalAffineMSA) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kOV,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(false);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, SemiGlobalMSAWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kOV,
        5, -4, -8, -8);

    initialize();
    construct_partial_order_graph(true);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}

TEST_F(SpoaAlignmentTest, SemiGlobalAffineMSAWithQualities) {
    SetUp(spoa_test_data_path + "sample.fastq", spoa::AlignmentType::kOV,
        5, -4, -8, -6);

    initialize();
    construct_partial_order_graph(true);

    std::vector<std::string> msa;
    graph->generate_multiple_sequence_alignment(msa);

    EXPECT_TRUE(msa.size() == sequences.size());

    for (uint32_t i = 0; i < msa.size(); ++i) {
        std::string tmp = "";
        for (const auto& c: msa[i]) {
            if (c != '-') tmp += c;
        }

        EXPECT_TRUE(tmp.size() == sequences[i]->data().size());
        EXPECT_TRUE(tmp.compare(sequences[i]->data()) == 0);
    }
}