A Fine-Tuned Universe

NCBI Assembly에서 다운로드 받은 파일이름 바꾸기 참고용으로 남겨두고 더 자세한 내용은 http://blog.genoglobe.com/2019/02/bash-cut.html [BASH] cut 명령을 이용하여 파일 이름을 간략하게 정리하기 cut - remove sections from each line of files (print selected parts of lines from each file to standard out) 커맨드 라인에서 man cut이라고 치면 나오는 설명... blog.genoglobe.com $ ls | while read f > do > cp $f $(echo $f | cut -d'_' -f 1,2).fna > done

세균 유전체 annotation을 위해 가장 많이 사용하는 방법은 NCBI의 prokaryotic genome annotation pipeline (PGAP, https://www.ncbi.nlm.nih.gov/genome/annotation_prok/) 일 것이다. FASTA 파일을 업로드하면서 PGAP를 이용해 annotation을 하겠다는 옵션에 체크만 해주면 신경쓸 필요없이 쉽게 annotation을 할 수 있고 여기서 부여된 locus tag, accession number 등은 전세계 공통으로 사용할 수 있어서 편리하다 하지만 별로도 유전체 annotation을 할 일도 생기는 데 (유전체 정보를 미리 공개하고 싶지 않은 경우 등) 이 때 Prokka를 사용할 수 있다. Prokka는 세균 유..

KEGG pathway mapper는 내가 분석하고자 하는 유전체 또는 유전자 세트가 KEGG pathway 상에서 어디에 해당하는지 색깔로 표시해주는 도구이다 https://www.genome.jp/kegg/mapper/ 1. amino acid fasta 준비 GenBank FTP에 들어가서 faa 파일을 다운로드 받는다 2. BlastKOALA에 faa를 업로드하고 이메일로 전송되는 결과를 다운로드 받는다 3. KEGG pathwaty mapper - reconstruct 에서 BlastKOALA의 결과파일을 붙여넣고 exec을 클릭한다

1. Upload data to galaxy (1) FTP를 이용하는 방법 Filezilla 설치 후 galaxy ftp에 접속, fastq.gz 파일을 업로드한다 Host: usegalaxy.org Username, Password 입력 usegalaxy.org 접속 후 FTP에 올려놓은 파일을 가져온다 2. Quality check FastQC를 default 값으로 수행 Short read data from your current history: input file 선택 Execute 클릭하면 output으로 raw data, HTML file 생성 3. Trimming Sickle sickle로 low-quality reads 제거 Single-end or paired-end reads: Sing..