def test_query_indexed_seqs_fasta(self): """pull out some sequences""" idx = index_fasta(self.seqs_fasta) seqs = query_indexed_seqs(['ABC', '123'], idx) exp = sorted([">ABC 567", "AATTGGCC", ">123", "CCCC"]) obs = sorted(seqs.splitlines()) self.assertEqual(obs, exp)
def test_query_indexed_seqs_fasta(self): """pull out some sequences""" idx = index_fasta(self.seqs_fasta) seqs = query_indexed_seqs(['ABC','123'], idx) exp = sorted([">ABC 567","AATTGGCC",">123","CCCC"]) obs = sorted(seqs.splitlines()) self.assertEqual(obs,exp)
def test_index_fasta(self): """Index a fasta!""" res = index_fasta(self.seqs_fasta) self.assertEqual(res.rsplit('/')[-1], "seqs.fasta.cidx")