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))