def test_unmapped_filter(self): reads = Reads(self.lambda_sam1_file) l = len(reads) unmapped_filter = UnmappedFilter( reads.add_mask_description("unmapped", "Filter unmapped reads")) reads.filter(unmapped_filter) assert len(reads) < l reads.close()
def test_contaminant_filter(self): reads = Reads(self.lambda_sam1_file) l = len(reads) contaminant = Reads(self.lambda_sam1_contaminant_file) contaminant_filter = ContaminantFilter(contaminant) reads.filter(contaminant_filter) assert len(reads) == l - 9 reads.close() contaminant.close()