def test_split_paired_end(self):
merged, r1, r2, r1_unpaired, r2_unpaired = [
get_data_path(i) + '.fq.gz' for i in
['merged', 'r1', 'r2',
'r1_unpaired', 'r2_unpaired']]
with TemporaryDirectory() as d, gzip.open(merged, 'rt') as fh:
fs = split_paired_end(fh, prefix=join(d, 'test'))
for i, j in zip([r1, r2, r1_unpaired, r2_unpaired], fs):
obs = list(skbio.io.read(i, format='fastq'))
exp = list(skbio.io.read(j, format='fastq'))
self.assertTrue(equal_seqs(obs, exp))
def test_equals_seqs(self):
seq1 = [DNA('ATGC'), DNA('A')]
seq2 = [DNA('A'), DNA('ATGC')]
self.assertTrue(equal_seqs(seq1, seq2))
