blast -m1
当database是10个物种(A、B、C、E、F、G、H、J、I、K)时,进行all vs all 比对结果是:
此时reference是物种A的第一个基因:即用10个物种的genome中的所有基因比对物种A的基因1
eg:如下
$ cat A_R00001_A.fa.70.blastp.m1|grep '1_0' -A 2 > r.sh Query= A_R00001_A
(423 letters) Database: all.10
165,229 sequences; 72,505,122 total letters Searching..................................................done Score E
Sequences producing significant alignments: (bits) Value A_R00001_A 855 0.0
B_R00001_B 766 0.0 1_0 1 PITLWLFYTHWLPGSPTTHHRVQGTTAMKDMVKCRVLKAHPCQECGKSFSKKGNLKRHQR 60
A_R00001_A 1 ............................................................ 60
B_R00001_B 4 ..l.......w..dAA....L..rv......................... 53
-- ... 1_0 361 GKCFPHKRHLIKHQLLHSRGGAYKCGVCGKRYRLKKYLRRHQKIHMREGTAPCSKPGYTT 420
A_R00001_A 361 ............................................................ 420
B_R00001_B 357 ..........................d..................t.........wee.. 416
-- 1_0 421 QTA 423
A_R00001_A 421 ... 423
B_R00001_B 417 R.G 419
1.Sequences producing significant alignments:大致情况,并计算score值,E值
关于score值:
he bit score gives an indication of how good the alignment is; the higher the score, the better the alignment. In general terms, this score is calculated from a formula that takes into account the alignment of similar or identical residues, as well as any gaps introduced to align the sequences.
2.具体情况,
1_0:表示reference的的第一个基因,即A物种的基因1
B_R00001_B:表示物种B的基因1与该reference(A物种的基因1) 的比对情况,从第4个AA开始到第53个AA结束
..l.......w..dAA....L..rv........................长度为53-4+1=50个;.是比对成功;空白是B物种基因1没有序列;60个为一组
如果物种A有10个基因,则将会把10个基因全部比对一遍:
$cat A_R00001_A.fa.70.blastp.m1|grep 'Query' >> r.sh Query= A_R00001_A
Query= A_R00002_A
Query= A_R00003_A
Query= A_R00004_A
Query= A_R00005_A
Query= A_R00006_A
Query= A_R00007_A
Query= A_R00008_A
Query= A_R00009_A
Query= A_R00010_A
blast -m1的更多相关文章
- RFID-RC522、FM1702SL、M1卡初探
catalogue . 引言 . RC522芯片(读卡器)简介 . FM1702SL芯片(读卡器)简介 . RFID M1卡简介 . 读取ID/序列号(arduino uno.MFRC522芯片 Ba ...
- M1/M2个人总结
软件工程整个学期结束了,很开心学了这门课,在学到知识的同时也提高了自己的动手实践的能力,感觉自己在整个软件工程的各个环节中都能有所把握,可以将学到的知识运用到设计.实践更多的项目中去. M1阶段个人总 ...
- CPU卡与M1卡的区别
简单来讲CPU卡比M1卡更加安全.扩展性更好.支持更多应用 CPU卡 M1 操作系统 带有COS系统 无COS系统 硬件加密模块 硬件DES运算模块 无实现算法的硬件加密模块 算法支持 标准DES ...
- the comment lines of the blast tabular format
进行Blast比对,用参数-m 6 可以以列表的方式输出结果,结果中从左到右每一列的意义分别是: [00] Query id [01] Subject id [02] % identity [03] ...
- Blast本地化
转载]Blast本地化:使用Blastall进行数据库比对 (2012-02-13 21:25:31) 用blastall进行序列比对 blastall是最常用的blast程序之一,其功能非常强大 ...
- 本地blast用法
格式化数据库: makeblastdb -in db.fasta -dbtype prot -parse_seqids -out dbname 参数说明: -in:待格式化的序列文件 -dbtype: ...
- M1卡介绍
本文整理自网络. M1卡是指菲利浦下属子公司恩智浦出品的芯片缩写,全称为NXP Mifare1系列,常用的有S50及S70两种型号,目前都有国产芯片与其兼容,属于非接触式IC卡.最为重要的优点是可读可 ...
- BLAST - 序列数据库搜索
我生信入门,老师就要求我学好blast比对,说得也确实是很有道理,是个人都知道比对是最基本的东西,现在再想想那老师的建议,也只能呵呵一笑. 北大生物信息公开课有一章专门讲得序列数据库搜索,可以好好看看 ...
- 关于求解不定方程的n(n-1)=2m(m-1)的解法的总结
主要参考下面两篇论文 1.<不定方程x(x-1)=Dy(y-1)的解法> 2.<用递推公式求一个不定方程的正整数解> 原有题目意思是 记得有一次全班去唱K, 其中有个活动是情歌 ...
随机推荐
- xhell ctrl+s 假死
如题. ctrl+q 即可接触 假死状态. 记录下
- 2019中国大学生程序设计竞赛(CCPC) - 网络选拔赛 A题
A - ^&^ Bit operation is a common computing method in computer science ,Now we have two positive ...
- Linux(CENTOS7) Redis安装
1.下载redis 在disk目录下,输入以下命令进行下载: wget http://download.redis.io/releases/redis-2.8.3.tar.gz 2.解 ...
- WIFI无线协议802.11a/b/g/n/ac的演变以及区别
摘自:https://blog.csdn.net/Brouce__Lee/article/details/80956945 毫无疑问,WiFi的出现普及带给我们巨大的上网便利,所以了解一下WiFi对应 ...
- 面向对象 part2 属性的特性
6.1.1理解对象 创建自定义对象最简单的方式就是创建一个object实例.然后添加方法和实例 var person = new Object() person.name = "hi&quo ...
- 洛谷 P3811 【模板】乘法逆元(欧拉定理&&线性求逆元)
题目传送门 逆元定义 逆元和我们平时所说的倒数是有一定的区别的,我们平时所说的倒数是指:a*(1/a) = 1,那么逆元和倒数之间的区别就是:假设x是a的逆元,那么 a * x = 1(mod p), ...
- dubbo的本地存根
在消费者创建存根类 修改消费者XML 也可以修改消费者注解
- Python——Pandas 时间序列数据处理
介绍 Pandas 是非常著名的开源数据处理库,我们可以通过它完成对数据集进行快速读取.转换.过滤.分析等一系列操作.同样,Pandas 已经被证明为是非常强大的用于处理时间序列数据的工具.本节将介绍 ...
- column命令
grep -E "car_flag|feaname" s_35926_uid_psi_table_20170407.csv | column -t -s, -o'|'
- mysql自定义函数初始化数据:init_data()
DELIMITER $$ USE `local_hnyz`$$ DROP FUNCTION IF EXISTS `init_data`$$ CREATE DEFINER=`root`@`localho ...