解析TMalign结构比对结果

TM-align 将两个蛋白质结构通过旋转和位移对齐后

TMalign 6owd_A.pdb  6owd_B.pdb -m 6owd_AB_mtx.txt 1>output.txt

屏幕输出内容,可以重定向为文件,文件内容为:


 *********************************************************************
 * TM-align (Version 20220412): protein structure alignment          *
 * References: Y Zhang, J Skolnick. Nucl Acids Res 33, 2302-9 (2005) *
 * Please email comments and suggestions to zhanglab@zhanggroup.org   *
 *********************************************************************

Name of Chain_1: 6owd_A.pdb (to be superimposed onto Chain_2)
Name of Chain_2: 6owd_B.pdb
Length of Chain_1: 28 residues
Length of Chain_2: 28 residues

Aligned length= 27, RMSD=   0.80, Seq_ID=n_identical/n_aligned= 0.000
TM-score= 0.85914 (if normalized by length of Chain_1, i.e., LN=28, d0=1.12)
TM-score= 0.85914 (if normalized by length of Chain_2, i.e., LN=28, d0=1.12)
(You should use TM-score normalized by length of the reference structure)

(":" denotes residue pairs of d <  5.0 Angstrom, "." denotes other aligned residues)
-QIARLQRQIRALQRQNARLQRQIRALQW
 ::::::::::::::::::::::::::: 
QIARLQRQIRALQRQNARLQRQIRALQW-

Total CPU time is  0.00 seconds
输出结果解释
  • Aligned length:两个结构的对齐长度为 27 个残基,这意味着虽然每条链都有 28 个残基,但比对过程中只有 27 个残基是成功对齐的。
BWA(Burrows-Wheeler Alignment)是一种常用的基因组比对工具,用于将测序数据与参考基因组进行比对。在使用BWA对样本进行比对后,可以根据比对结果来计算比对率。下面是一种常用的统计方法: 首先,从BWA比对输出的文件中提取比对信息。BWA比对结果一般以SAM(Sequence Alignment/Map)格式保存,可以使用SAMtools等相关工具对比对结果进行解析和提取。 从SAM文件中,可以获得每个比对结果的相关信息,如比对的状态(比对成功或失败)、比对上的read数量等。 比对率通常被定义为成功比对的read数量占总read数量的比例。因此,可以通过计算成功比对的read数量除以总read数量得到比对率。 具体的计算步骤如下: 1. 从SAM文件中按行读取每个比对结果。 2. 对于每行结果,判断比对的状态是否为成功(比对状态通常在每行的FLAG列中标注)。 3. 如果比对成功,则成功比对的read数量加1。 4. 继续读取下一行并重复上述步骤,直到读取完所有比对结果。 5. 计算成功比对的read数量除以总read数量,并乘以100%得到比对率(以百分比表示)。 需要注意的是,统计比对率时,应该使用有效的read数量。在有些情况下,比如双端测序,一对read中的其中一条无法比对成功,但另一条可以成功比对,这种情况下,只统计比对成功的read数量。 总之,通过解析BWA比对结果的SAM文件,并统计成功比对的read数量除以总read数量,就可以得到比对率。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值