Exemplo n.º 1
0
def test_minbasequality():
	filter = MinBaseQuality(30)
	seq_len_a = SeqRecord (Seq ('ACGT'), letter_annotations={'phred_quality': [10, 20, 30, 40]})
	assert filter (seq_len_a) == False
	seq_len_b = SeqRecord (Seq ('ACGT'), letter_annotations={'phred_quality': [30, 30, 30, 40]})
	assert filter (seq_len_b) == True
	seq_len_c = SeqRecord (Seq ('ACGT'), letter_annotations={'phred_quality': [40, 40, 40, 40]})
	assert filter (seq_len_c) == True
Exemplo n.º 2
0
def test_minlength():
	filter = MinLength(4)
	seq_len3 = SeqRecord (Seq ('ACG'))
	assert filter (seq_len3) == False
	seq_len4 = SeqRecord (Seq ('ACGT'))
	assert filter (seq_len4) == True
	seq_len5 = SeqRecord (Seq ('ACGTA'))
	assert filter (seq_len5) == True