s.load_fasta(seq_path);
BOOST_CHECK_EQUAL(s.subseq(0, 5), "GGATC"); // first few chars of fasta file
BOOST_CHECK_EQUAL(s.subseq(2, 3), "ATC");
- BOOST_CHECK_EQUAL(s.get_header(), "gi|180579|gb|M21487.1|HUMCKMM1 Human "
+ BOOST_CHECK_EQUAL(s.get_fasta_header(), "gi|180579|gb|M21487.1|HUMCKMM1 Human "
"muscle creatine kinase gene (CKMM), "
"5' flank");
}
/* FIXME: enable this when i find a way of passing storing the motif name
// does our annotation travel?
Sequence motif_seq(m);
- motif_seq.set_header("hi");
+ motif_seq.set_fasta_header("hi");
s1.add_motif(motif_seq);
BOOST_CHECK_EQUAL(s1.motifs().size(), 2);
BOOST_CHECK_EQUAL( s, buf.str() );
}
+BOOST_AUTO_TEST_CASE( get_name )
+{
+ Sequence seq("AAGGCCTT");
+
+ BOOST_CHECK_EQUAL( seq.get_name(), "" );
+ seq.set_species("hooman"); // anyone remember tradewars?
+ BOOST_CHECK_EQUAL( seq.get_name(), "hooman");
+ seq.set_fasta_header("fasta human");
+ BOOST_CHECK_EQUAL( seq.get_name(), "fasta human");
+}