#include <boost/filesystem/operations.hpp>
namespace fs=boost::filesystem;
+#include <boost/algorithm/string/case_conv.hpp>
+
#include <list>
#include <iostream>
#include <sstream>
BOOST_CHECK_EQUAL(s2.rev_comp(), "GGCCAATT");
BOOST_CHECK_EQUAL(s2.rev_comp(), "ggccaatt"); // we should be case insensitive now
BOOST_CHECK_EQUAL(s2.size(), s2.size());
- BOOST_CHECK_EQUAL(s2.get_sequence(), "aattggcc");
+ //We're currently forcing sequences to uppercase
+ BOOST_CHECK_EQUAL(s2.get_sequence(), "AATTGGCC");
Sequence s3("asdfg", Sequence::reduced_dna_alphabet);
BOOST_CHECK_EQUAL(s3, "ANNNG");
std::string bdv("BDv");
Sequence seq_bdv(bdv, Sequence::nucleic_alphabet);
BOOST_CHECK_EQUAL(seq_bdv.size(), bdv.size());
- BOOST_CHECK_EQUAL(seq_bdv.get_sequence(), bdv);
+ // forcing sequence to upper case
+ BOOST_CHECK_EQUAL(seq_bdv.get_sequence(),
+ boost::algorithm::to_upper_copy(bdv));
}
BOOST_CHECK( seq != seqr );
BOOST_CHECK_EQUAL( seq, seqr.rev_comp() );
- BOOST_CHECK_EQUAL( seq.get_sequence(), iupac_symbols );
+ // forcing sequence to upper case
+ BOOST_CHECK_EQUAL( seq.get_sequence(),
+ boost::algorithm::to_upper_copy(iupac_symbols) );
}
BOOST_AUTO_TEST_CASE( sequence_reverse_complement_dna )