コード例 #1
0
ファイル: test_seq.py プロジェクト: Darlingfuer/fanc
    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()
コード例 #2
0
ファイル: test_seq.py プロジェクト: Darlingfuer/fanc
    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()