Commit 03bc5707 authored by Armin Sobhani's avatar Armin Sobhani
Browse files

add commented sample for partial reverse complement of a sequence

parent ddaa97bb
......@@ -25,6 +25,7 @@
#include <boost/range/algorithm/copy.hpp>
#include <boost/range/adaptor/reversed.hpp>
#include <boost/range/adaptor/filtered.hpp>
#include <boost/range/adaptor/replaced.hpp>
#include <boost/bioseq/type_traits.hpp>
#include <boost/bioseq/sequence.hpp>
......@@ -166,6 +167,21 @@ BOOST_AUTO_TEST_CASE_TEMPLATE(read_multi_fasta, T, seq_types)
// BOOST_CHECK_EQUAL(22, reads.size());
// BOOST_CHECK_EQUAL(reads.back().name(), "gi|379009272|ref|NC_016893.1|");
// for (auto read : reads | boost)
// std::cout << read.name() << '\n' << read.data() << std::endl;
// for (auto read : reads)
// std::cout << read.name() << '\t' << read.size() << std::endl;
// {
// std::cout << read.name() << '\n'
// << read.data()
// << std::endl;
// boost::copy(
// read | boost::adaptors::replaced('A', 'T')
// // | boost::adaptors::replaced('T', 'A')
// | boost::adaptors::replaced('C', 'G')
// // | boost::adaptors::replaced('G', 'C')
// | boost::adaptors::reversed
// , std::ostream_iterator<char>(std::cout) );
// std::cout << std::endl;
// }
}
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment