Пример #1
0
    def test_pair_direction_and_name(self):
        "it test the pair_name parser"
        title = "seq8:136:FC706VJ:2:2104:15343:197393 1:Y:18:ATCACG"
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == "seq8:136:FC706VJ:2:2104:15343:197393"
        assert dir_ == FWD

        title = "seq8:136:FC706VJ:2:2104:15343:197393/1"
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == "seq8:136:FC706VJ:2:2104:15343:197393"
        assert dir_ == FWD

        title = "seq8:136:FC706VJ:2:2104:15343:197393.f"
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == "seq8:136:FC706VJ:2:2104:15343:197393"
        assert dir_ == FWD

        title = "seq8:136:FC706VJ:2:2104:15343:197393.mp12"
        try:
            name, dir_ = _parse_pair_direction_and_name_from_title(title)
            self.fail()
        except PairDirectionError:
            pass

        title = r"seq8:136:FC706VJ:2:2104:15343:197393\1"
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == "seq8:136:FC706VJ:2:2104:15343:197393"
        assert dir_ == FWD

        # With SeqRecord
        seq = SeqRecord(id=r"seq8:136:FC706VJ:2:2104:15343:197393\1", seq=Seq("ACT"))
        name, dir_ = _parse_pair_direction_and_name(SeqWrapper(SEQRECORD, seq, None))
        assert name == "seq8:136:FC706VJ:2:2104:15343:197393"
        assert dir_ == FWD
    def test_pair_direction_and_name():
        'it test the pair_name parser'
        title = '@seq8:136:FC706VJ:2:2104:15343:197393 1:Y:18:ATCACG'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == '@seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        title = '@seq8:136:FC706VJ:2:2104:15343:197393/1'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == '@seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        title = '@seq8:136:FC706VJ:2:2104:15343:197393.f'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == '@seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD
Пример #3
0
    def test_pair_direction_and_name(self):
        'it test the pair_name parser'
        title = 'seq8:136:FC706VJ:2:2104:15343:197393 1:Y:18:ATCACG'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        title = 'seq8:136:FC706VJ:2:2104:15343:197393/1'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        title = 'seq8:136:FC706VJ:2:2104:15343:197393.f'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        title = 'seq8:136:FC706VJ:2:2104:15343:197393.mp12'
        try:
            name, dir_ = _parse_pair_direction_and_name_from_title(title)
            self.fail()
        except PairDirectionError:
            pass

        title = r'seq8:136:FC706VJ:2:2104:15343:197393\1'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        # With SeqRecord
        seq = SeqRecord(id=r'seq8:136:FC706VJ:2:2104:15343:197393\1',
                        seq=Seq('ACT'))
        name, dir_ = _parse_pair_direction_and_name(
            SeqWrapper(SEQRECORD, seq, None))
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD
Пример #4
0
    def test_pair_direction_and_name(self):
        'it test the pair_name parser'
        title = 'seq8:136:FC706VJ:2:2104:15343:197393 1:Y:18:ATCACG'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        title = 'seq8:136:FC706VJ:2:2104:15343:197393/1'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        title = 'seq8:136:FC706VJ:2:2104:15343:197393.f'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        title = 'seq8:136:FC706VJ:2:2104:15343:197393.mp12'
        try:
            name, dir_ = _parse_pair_direction_and_name_from_title(title)
            self.fail()
        except PairDirectionError:
            pass

        title = r'seq8:136:FC706VJ:2:2104:15343:197393\1'
        name, dir_ = _parse_pair_direction_and_name_from_title(title)
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD

        # With SeqRecord
        seq = SeqRecord(id=r'seq8:136:FC706VJ:2:2104:15343:197393\1',
                        seq=Seq('ACT'))
        name, dir_ = _parse_pair_direction_and_name(SeqWrapper(SEQRECORD, seq,
                                                               None))
        assert name == 'seq8:136:FC706VJ:2:2104:15343:197393'
        assert dir_ == FWD