Blast를 하기 위해서는 query와 database가 있어야 합니다.
그러므로 어떤 유전체에 대한 all-vs-all blast를 하려면 일단 그 유전체를 사용해 database를 만들고 다시 그 유전체를 query로 사용하여 blast를 수행하면 됩니다.
Database를 만드는 방법:
>makeblastdb -in input_database(query)_file.fasta -dbtype prot -out my_db_file
그다음 똑같은 유전체를 query로 blast를 수행하면 됩니다.
>blastp -db my_db_file -query input_database(query)_file.fasta -outfmt 6 -out all-vs-all.tsv -num_threads 4
그럼 완료
출처
http://www.biostars.org/p/6541/
'생물정보학 > Bioinformatics' 카테고리의 다른 글
[논문] 출판물과 유전자 및 단백질 데이터를 서로 엮기 (0) | 2014.02.05 |
---|---|
BLAST에서 BLAST+로 명령어 바꾸기 (0) | 2014.01.22 |
박사 학위 또는 포닥을 Bioinformatics/Computational Biology로 해야하는 탑 N 가지 이유 (0) | 2013.12.20 |
[스크랩] 초보도 할수있느 조립식PC 용어이해 및 견작짜보기. (마지막수정) (0) | 2013.08.07 |
UltraEdit HTML tag removing regexp (0) | 2013.05.15 |