def test_call_consensus_simple(self):
res = assemble.call_consensus(self.ref, [self.aln1, self.aln2],
self.args.cons_cov,
self.args,
strict=True)
exp = 'NNNNNNNNNNNNNAANNNNN'
self.assertEqual(res, exp)
def test_call_consensus_no_aligned_sequences(self):
res = assemble.call_consensus(self.ref, [],
self.args.cons_cov,
self.args,
strict=True)
exp = ''
self.assertEqual(res, exp)
def test_call_consensus_disagreement_coverage_requirement(self):
res = assemble.call_consensus(self.ref,
[self.aln1, self.aln2, self.aln3],
self.args.cons_cov,
self.args,
strict=True)
exp = 'NNNNNNNNNNNNNAANAANN'
self.assertEqual(res, exp)
def test_call_consensus_high_coverage_requirement(self):
self.args.cons_cov = 3
res = assemble.call_consensus(self.ref, [self.aln1, self.aln2],
self.args.cons_cov,
self.args,
strict=True)
exp = 'NNNNNNNNNNNNNNNNNNNN'
self.assertEqual(res, exp)