반응형
read 서열에서 error를 1이하로 허용하는 내에 BESTMATCH를 찾아 시작과 종료지점 그리고 매치되는 서열을 확인하는 코드. error는 mismatch, insertion, deletion을 의미한다.
import regex
primer_seq, read_seq
regex_primer_seq = fr'({primer_seq}{{e<=1}})'
match_object = regex.search(regex_primer_seq, read_seq, regex.BESTMATCH)
match_start, match_end = match_object.span()
match_seq = match_object.captures()반응형
'Computer Science > python' 카테고리의 다른 글
| pandas 활용하기 (0) | 2022.02.18 |
|---|---|
| logging 모듈 사용하기 (0) | 2022.02.17 |
| Primer 서열 분석을 위한 python 코드 (0) | 2021.08.17 |
| String Format으로 길이 고정하기 (0) | 2020.06.24 |
| python multi-level argparse (0) | 2019.07.12 |