使用原理在网上很多,基本基于Blastn,
tblastx基本通用。
参考:
https://bi.biopapyrus.jp/seq/blast/tblastx.html
为了达成all vesus all的blast,首先建立本地数据库。建立本地数据库之前,设置一下电脑的blast path路径,不然无法调用。具体请查询网络。
C:\Users\yemao>cd D:\blast\db
C:\Users\yemao>d:
D:\blast\db>makeblastdb.exe -in new1.fasta -parse_seqids -hash_index -dbtype nucl -out D:\blast\db\new1
调用makeblastdb函数,将new1.fasta建立成本地数据库,(剩下的就是一些参数),序列选择是核苷酸(nucl),如果是蛋白的话再查一下,out出库。建立好数据库后显示:
Building a new DB, current time: 01/04/2021 16:29:38
New DB name: D:\blast\db\new1
New DB title: new1.fasta
Sequence type: Nucleotide
Keep MBits: T
Maximum file size: 1000000000B
Adding sequences from FASTA; added 416 sequences in 0.245911 seconds.
就是建立成功。每一个序列的名称不要太长,不然会提示你 too long。
而后就是tblastx序列比对:
D:\blast\db>tblastx -db new1 -query new1.fasta -matrix BLOSUM45 -outfmt 6 -out result.outfm6 -evalue 1e-3
-matrix 默认是BLOSUM62 -outfmt的格式1-10,可以查询https://blog.csdn.net/dpsh_6/article/details/88382137?utm_medium=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-4.control&depth_1-utm_source=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-4.control
或者自己调用 tblastx -help查看也行