A Fine-Tuned Universe

FastANI 결과를 보면 아래와 같이 나온다 즉 genome1 genome2 ANI값 이런 식으로 나온다 하지만 보통은 matrix 형태로 정리된 아래 형태에 익숙할 것이다. 아니면 아래와 같은 형식으로 바꿔서 결과를 정리해야 할 때가 있다 그런데 비교하려는 유전체가 많아질 수록 결과를 정리하기가 힘들다. 그래서 주말에 아래와 같이 쉘스크립트를 써서 결과를 정리하였다 먼저 유전체 목록을 아래와 같이 파일로 만들었다. 사실 둘 다 같은 내용인데 가로 세로를 구분해서 생각하려고 따로 만들었다 ls *.gz > query_list ls *.gz > reference_list 그리고 아래와 같이 스크립트를 짜서 실행시켰다. #!/bin/bash while read query do awk -v genome=$q..

Polysaccharide utilization loci (PUL)은 세균 유전체 상에서 다당류 물질을 이용하기 위해 필요한 유전자가 모여있는 부분을 말한다. 대부분 Bacteroidetes phylum의 세균에서 주로 발견된다. Pseudomonadota (전에는 Proteobacteria 였는데 이름이 바뀌었다. 적응이 안된다.) 세균들은 다당류를 이용하더라도 PUL 형태의 gene cluster가 거의 없는 것 같다. PUL에는 다당류 분해에 필요한 glycoside hydrolase, polysaccharide lyase 와 같은 carbohydrate-active enzymes (CAZy) 뿐만 아니라 transcriptional regulator, transporter 등 유전자가 함께 존재하..

Bioinformatics를 잘 하시는 분들은 스크립트로 여러가지 일을 뚝딱뚝딱 하시겠지만 많은 생물학 기반의 연구자와 대학원생들은 갑자기 유전체나 시퀀스 분석을 해야하는 상황을 만나면 당황할 수 밖에 없다. 많은 bioinformatics tool이 리눅스에서 돌아가는데 리눅스도 모르고 설치하는 것도 잘 모르겠고 각종 프로그래밍 언어도 모르고... 아무것도 모르는 상태에서 배워가면서 사용하면서 부딪히는 문제들을 두서없이 메모해서 업데이트하고 있다. 나 같은 사람이 또 있다면 도움이 되길... Official Roary sites https://sanger-pathogens.github.io/Roary/ https://github.com/microgenomics/tutorials/blob/master/p..

지난 번에 nucmer header 에 대해 정리했는데 이번엔 mummer 결과파일의 header에 정리해놓으려고 한다. 좋은 프로그램 같긴 한데 쓰려면 알아야 할 게 많은 것 같다. (결과 예시) > ID1 4655667 1 31 4655699 33 319 4656019 353 520 4656540 874 20 > ID1 Reverse 741743 22 872 (설명) > 뒤에는 시퀀스ID 가 나온다 (inpu..