+BOOST_AUTO_TEST_CASE( seqspan_global_mutable_start_stop_minus_strand )
+{
+ std::string seq_string("AAAAGCTA");
+ SeqSpanRef s1(new SeqSpan(seq_string));
+
+ SeqSpanRef s2 = s1->subseq(2,3, SeqSpan::MinusStrand);
+ BOOST_CHECK_EQUAL( s2->start(), 2);
+ BOOST_CHECK_EQUAL( s2->stop(), 2+3);
+ BOOST_CHECK_EQUAL( s2->size(), 3);
+ BOOST_CHECK_EQUAL( s2->sequence(), "CTT");
+
+ SeqSpanRef s3 = s2->subseq(1,2, SeqSpan::SameStrand);
+ BOOST_CHECK_EQUAL(s3->sequence(), "TT");
+
+ // Could also argue that it should be CT
+ // if you assume that the locations are all relative to the global sequence
+ // and are then reverse complemented
+
+ s2->setStart(1);
+ BOOST_CHECK_EQUAL( s2->sequence(), "CTTT");
+}
+