def test_reverse_complement(): 'It test the reverse complement fucntion' sequence1 = Seq('aaaccttt') # First we initialice the quality in the init seq1 = SeqWithQuality(name='seq1', seq=sequence1, \ qual=[2, 4 , 1, 4, 5, 6, 12, 34]) seq2 = reverse_complement(seq1) assert seq2.seq == sequence1.reverse_complement() assert seq2.qual == seq1.qual[::-1]