PYTHON๐Ÿ 4

[๋ฐ”์ด์˜คํŒŒ์ด์ฌ] biopython ์„ค์น˜ ๋ฐฉ๋ฒ•

NGS ๋ฐ์ดํ„ฐ๋ฅผ ๋‹ค๋ฃจ๋‹ค๋ณด๋ฉด Fastq, BAM, VCF ๋“ฑ ๋‹ค์–‘ํ•œ ํ˜•์‹์˜ ํŒŒ์ผ์„ ์ ‘ํ•˜๊ฒŒ ๋ฉ๋‹ˆ๋‹ค.์˜ค๋Š˜์€ ๊ทธ์ค‘์—์„œ Fastq ํŒŒ์ผ์„  biopython์„ ์ด์šฉํ•˜์—ฌ ๋‹ค๋ค„๋ณด๊ฒ ์Šต๋‹ˆ๋‹ค. 1. Biopython ์ด๋ž€? ์ฐธ๊ณ  ๋งํฌ : https://biopython.org/ ๋ฐ”์ดํŒŒ์ด์ฌ์ด๋ž€ ํŒŒ์ด์ฌ ๋ผ์ด๋ธŒ๋Ÿฌ๋ฆฌ์ค‘ ํ•˜๋‚˜๋กœ, ๋ถ„์ž์ƒ๋ฌผํ•™ ๊ณ„์‚ฐ์— ์ž์œ ๋กญ๊ฒŒ ์‚ฌ์šฉํ•  ์ˆ˜ ์žˆ์Šต๋‹ˆ๋‹ค. ์œ ์ „์ฒด ๋ฐ์ดํ„ฐ ์„œ์—ด์„ ์ฝ๊ฑฐ๋‚˜ ์ƒ๋ฌผ ์ •๋ณด ๋ฐ์ดํ„ฐ๋ฅผ ๋‹ค๋ฃจ๋Š”๋ฐ ์ฃผ๋กœ ์‚ฌ์šฉํ•ฉ๋‹ˆ๋‹ค. ์•„๋ž˜ ํ˜•์‹๋“ค์˜ ํŒŒ์ผ๋“ค์€ Biopython์—์„œ ์ง€์›ํ•˜๋Š” ํŒŒ์ผ ํ˜•์‹๋“ค ์ž…๋‹ˆ๋‹ค.- Blast output- Clustalw- FASTA- GeneBank- PubMed- ExPASy files- SCOP- UniGene- SwissProt 2. Biopython ์„ค์น˜๋ฐฉ๋ฒ• (LINUX, MA..

PYTHON๐Ÿ 2024.05.07

[Tool] ArtDeCo tool ์ด์šฉํ•ด contamination ์•Œ์•„๋‚ด๊ธฐ

Sample contamination ์ด๋ž€ NGS ์‹คํ—˜ํ•˜๋Š” ์‚ฌ๋žŒ์ด๋ผ๋ฉด ๋ˆ„๊ตฌ๋‚˜ ํ•œ๋ฒˆ์ฏค ์ƒ๊ฐํ•ด๋ณผ๋ฒ•ํ•œ ์ด์Šˆ์ด๋‹ค. ์‹คํ—˜ํ•˜๋‹ค๊ฐ€ ์˜†์— well์— contam ๋˜๊ธฐ๋Š” ๋งค์šฐ ์‰ฝ์ง€๋งŒ, ์–ด๋–ค ์ƒ˜ํ”Œ์— ์˜ํ•ด contam์ด ๋˜์—ˆ๋Š”์ง€ ์•Œ์•„๋‚ด๊ธฐ๋Š” ์–ด๋ ต๋‹ค. ์ด๋Ÿฌํ•œ contamination issue๋ฅผ ํ•ด๊ฒฐํ•˜๊ธฐ ์œ„ํ•ด ๋ช‡๊ฐ€์ง€ tool์„ ์ฐพ์•„๋ณธ ๊ฒฐ๊ณผ ArtDeCo๋ผ๋Š” tool์„ ๋ฐœ๊ฒฌํ•˜์˜€๋‹ค. 1. ArtDeCo ๋ž€? NGS๋ฅผ ์ด์šฉํ•œ DNA sample ์‚ฌ์ด์˜ cross-contamination์„ ์•Œ์•„๋‚ผ ์ˆ˜ ์žˆ๋„๋ก ๊ฐœ๋ฐœํ•œ tool ์ด๋‹ค. ์•„๋ž˜ ๋…ผ๋ฌธ์„ ๋ณด๋ฉด batch ๋‚ด ์ƒ˜ํ”Œ๋“ค์—๊ฒŒ์„œ ์ค‘๋ณต๋˜๋Š” SNV์˜ AR(Alleic ratio) ํ™•์ธํ•˜์—ฌ ์ƒ˜ํ”Œ๋“ค์˜ Contamination์„ ํ™•์ธํ•œ๋‹ค. https://bmcbioinformatics.biomed..

PYTHON๐Ÿ 2024.03.22

[ํŒŒ์ด์ฌ ๊ธฐ์ดˆ] ํŒŒ์ด์ฌ for ๋ฌธ ๊ธฐ์ดˆ ์‚ฌ์šฉ๋ฒ•(range, dictionary, ๋Œ€์†Œ๋ฌธ์ž, break, continue)

ํŒŒ์ด์ฌ for๋ฌธ์„ ์ด์šฉํ•˜์—ฌ ๋ช‡๊ฐ€์ง€ ์œ ์šฉํ•œ ๋ฐฉ๋ฒ•๋“ค์„ ์•Œ์•„๋ณด์ž. 1. range ๋ฅผ ์‚ฌ์šฉํ•œ for ๋ฌธ ์ด์šฉ range๋ฅผ ์ด์šฉํ•˜๋ฉด ์‹œ์ž‘๋ถ€ํ„ฐ ๋๊ฐ’ ๊นŒ์ง€ ์‚ฌ์šฉ์ž๊ฐ€ ๋ฒ”์œ„๋ฅผ ์„ค์ • ๊ฐ€๋Šฅํ•˜๋‹ค. for๋ฌธ์„ ์ด์šฉํ•˜์—ฌ ์›ํ•˜๋Š” ๋ฒ”์œ„์˜ ์ˆซ์ž, ํ™€์ˆ˜, ์ง์ˆ˜ ๋“ฑ ์ถœ๋ ฅ ๊ฐ€๋Šฅํ•˜๋‹ค. #for ์‹ค์Šต for v1 in range(4): print('v1 is:', v1) for v2 in range(1,11): print('v2 is :' , v2) for v3 in range(0,11,2): print('v3 is:', v3) sum=0 for v in range(1,1001): sum += v print(sum) print(type(range(1,11))) 2. Dictionary๋ฅผ ์‚ฌ์šฉํ•œ for ๋ฌธ ์ด์šฉ Dictionary์˜ key ๊ฐ’..

PYTHON๐Ÿ 2024.03.13

[python] ๋ฆฌ๋ˆ…์Šค์—์„œ ํ˜„์žฌ ์‚ฌ์šฉํ•˜๊ณ  ์žˆ๋Š” python ๋ฒ„์ „ ๋ฐ ๊ฒฝ๋กœ ์•Œ์•„๋‚ด๊ธฐ

1. python ๋ฒ„์ „ ์•Œ์•„๋‚ด๊ธฐ python --version 2.python ์„ค์น˜๋œ ๊ฒฝ๋กœ ์•Œ์•„๋‚ด๊ธฐ ์ผ๋‹จ ๋ฆฌ๋ˆ…์Šค์—์„œ $python์„ ์น˜๊ณ  python์„ ์‹คํ–‰ํ•œ๋‹ค. $python import sys sys.executable sys.executable์„ ์น˜๋ฉด ํ˜„์žฌ ์‚ฌ์šฉํ•˜๊ณ  ์žˆ๋Š” ํŒŒ์ด์ฌ์˜ ๊ฒฝ๋กœ๊ฐ€ ๋‚˜์˜จ๋‹ค.

PYTHON๐Ÿ 2023.10.23