def test_seqs_string():
''' Here we test it we can initialice a seq with quality
and if we can print it. Seq is going to be a normal string'''
#sequence1 = Seq('aaavvttt')
# First we initialice the quality in the init
seq1 = SeqWithQuality(name='seq1', seq=Seq('aaavvttt'),
qual=[2, 4 , 1, 4, 5, 6, 12, 34])
assert seq1
# Here we add the quality after the initialization
seq2 = SeqWithQuality(name='seq2', seq=Seq('aaavvttt'))
seq2.qual = [2, 4 , 1, 4, 5, 6, 12, 34]
assert seq2
# Let's check the sliceability of the seq
for num in range(len(seq1)):
seq3 = seq1[num]
assert seq3
def test_seq_seq():
''' We are going to check the same tests but with a BioPython seq
class object inside seq'''
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])
assert seq1
# Here we add the quality after the initialization
seq2 = SeqWithQuality(name='seq2', seq=sequence1)
seq2.qual = [2, 4 , 1, 4, 5, 6, 12, 34]
assert seq2
# We check if the seq can be complemented
seq3 = seq2.complement()
assert seq3.seq == 'tttggaaa'
# Let's check the sliceability of the seq
assert seq2[0:2].seq == 'aa'
assert seq2[0:2].qual == [2, 4]