def test_query_indexed_seqs_fastq(self):
"""pull out some sequences"""
idx = index_fastq(self.seqs_fastq)
seqs = query_indexed_seqs(['ABC','123'], idx)
exp = sorted(["@ABC 567","AATTGGCC","+","12345678",
"@123","CCCC","+","9876"])
obs = sorted(seqs.splitlines())
self.assertEqual(obs,exp)
def test_query_indexed_seqs_fastq(self):
"""pull out some sequences"""
idx = index_fastq(self.seqs_fastq)
seqs = query_indexed_seqs(['ABC', '123'], idx)
exp = sorted([
"@ABC 567", "AATTGGCC", "+", "12345678", "@123", "CCCC", "+",
"9876"
])
obs = sorted(seqs.splitlines())
self.assertEqual(obs, exp)
def test_index_fastq(self):
"""Index a fastq!!!"""
res = index_fastq(self.seqs_fastq)
self.assertEqual(res.rsplit('/')[-1], "seqs.fastq.cidx")
def test_index_fastq(self):
"""Index a fastq!!!"""
res = index_fastq(self.seqs_fastq)
self.assertEqual(res.rsplit('/')[-1], "seqs.fastq.cidx")
