把几条同一基因的序列放进MacVector，比对后却出现成片差异，甚至原本熟悉的保守区域也对不上，这时很容易怀疑参数设错了。不过，序列方向相反、比较范围不一致、短片段缺少重叠，都可能产生类似结果。调整参数之前，应当看清错配分布在哪里，以及参与比较的序列是否适合放在一起。以下操作采用英文界面，部分选项的位置会随版本有所变化。

　　一、MacVector如何设置多序列比对参数

　　多序列比对适合比较有相关性的核酸或蛋白序列。建立文件时选对序列类型，再保留一份默认参数结果，后面的调整才有可比较的依据。

　　1、建立对应类型的比对文件

　　核酸和蛋白需要分别建立比对文件。同一批材料中，完整基因、编码区和局部测序片段也应先区分，避免把范围差异误当成序列差异。

　　①点击【File】→【New】，按材料类型选择【Nucleic Acid Alignment】或【Protein Alignment】。

　　②点击【Edit】→【Add Sequences From File】，加入准备好的序列文件。

　　③核对各行名称、序列长度及内容，删除误导入的重复记录。

　　④点击【File】→【Save As】，保存原始比对文件。

　　2、选择算法并建立对照结果

　　①在比对窗口中，点击【Analyze】→【Align Multiple Sequences Using】→【ClustalW】。

　　②在【Sequences to Align】中确认参与比对的序列，保留默认参数，点击【OK】。

　　③保存结果副本，再从【Align Multiple Sequences Using】中选择【Muscle】，比较相同材料的排列。

　　两种算法给出相同的保守区域，可以增加对该区域的信心。若差异集中在重复片段或缺口附近，这些位置还需要单独判断，不能仅凭某次结果看起来整齐就接受。

　　3、分别调整缺口相关参数

　　①重新打开【ClustalW】设置窗口，记录当前【Open Gap Penalty】与【Extend Gap Penalty】。

　　②对连续错配但缺口很少的区域，适当降低【Open Gap Penalty】，重新运行并对照原结果。

　　③对零散缺口较多的结果，适当提高【Open Gap Penalty】；需要比较缺口延伸情况时，单独调整【Extend Gap Penalty】。

　　④每次只改一项，点击【OK】运行，并将参数值记入结果文件名。

　　二、MacVector比对结果出现大量错配如何判断

　　缺口开启罚分提高后，算法通常更少引入新缺口；延伸罚分提高后，延长已有缺口的代价会增加。这些设置改变的是排列取舍，无法修正输入材料本身的问题。排查时可以按错配出现的位置缩小范围。

　　1、只有某条核酸序列整行对不上

　　同一基因的序列可能来自不同链方向。多序列比对前，应统一核酸序列方向，单纯倒转字符顺序不能代替反向互补。

　　①单独打开异常序列，进入【Editor】，核对序列来源和链方向。

　　②对需要转换的核酸序列全选，点击【Edit】→【Reverse&Complement】。

　　③使用【Save As】另存转换后的文件，将其加入新的比对副本后重算。

　　2、差异主要集中在两端

　　①打开原始序列的【Editor】，核对两端是否包含引物、载体片段或较长的【N】区段。

　　②对照原始测序记录，选中有依据需要去除的区域，在副本中删除。

　　③通过【Add Sequences From File】导入整理后的序列，再次比对共同覆盖的区域。

　　短序列两端的空白，可能只是没有测到相应位置。缺口、未知碱基和明确的碱基替换应分别看待。没有测序依据时，不要为了减少差异直接修改碱基，也不要把真实存在的长度差异全部裁掉。

　　3、从某个位置开始连续错开

　　一个插入或缺失附近，可能出现较长的错配带。尤其是连续相同碱基或重复片段，缺口放在不同位置都可能得到相近的排列，视觉上却差别很大。

　　①在比对编辑器中找到错配带起点，向前查看相邻的保守位置。

　　②打开对应原始序列的【Editor】，核对该位置的实际字符及测序记录。

　　③返回【ClustalW】，只调整【Open Gap Penalty】，比较缺口位置与两侧保守区域。

　　④对仍无法确定的区域保留原始数据，记录起止位置，暂缓解释其中的变异。

　　4、整组序列都缺少稳定匹配

　　①点击比对窗口的【Views】，选择【Identity/Similarity Matrix】，查看各对序列的比较结果。

　　②从【Sequences to Align】中暂时只保留两条来源清楚、范围对应的序列，重新运行。

　　③逐条加入其余序列，记录保守区域开始明显变化时加入的材料。

　　蛋白序列中的相似性可以包含性质相近的氨基酸替换，相同性则要求字符一致，两项不能混用。某条序列与其他材料普遍差异较大时，应核对基因名称、物种、转录本及截取范围。相同名称也可能对应不同片段，不能直接认定它们覆盖同一段区域。

　　三、MacVector怎样处理不适合直接多序列比对的材料

　　有些错配来自任务选择。用短测序片段核对已知模板，与比较完整蛋白之间的差异，处理方式应有所区别。持续调节缺口罚分，未必能够解决材料之间缺少重叠的问题。

　　1、核对短读段时改用参考序列比对

　　几条测序读段可能分别覆盖模板的不同位置，彼此没有足够重叠。把它们直接放进多序列比对，容易出现被勉强排在同一区域的情况。

　　①打开已确认的参考序列，点击【Analyze】→【Align to Reference】。

　　②在参考比对窗口中点击【+】，加入原始测序文件。

　　③点击【Align】，查看各读段在模板上的覆盖位置，并核对错配处的原始峰图。

　　2、蛋白片段不重叠时比较另一组设置

　　①将全长蛋白和对应片段加入同一个【Protein Alignment】文件。

　　②点击【Analyze】→【Align Multiple Sequences Using】→【Muscle】。

　　③在提供相应选项的设置窗口中，将【Profile】设为【PAM200】，将【Diagonals】设为【On】，点击【OK】。

　　④对照全长蛋白中的已知位置，核查各片段是否落入对应区域。

　　这组设置适用于检查不重叠蛋白片段与全长序列的排列，并非所有材料都应照搬。若当前版本没有这些选项，应查阅对应版本说明。即使片段排列改善，也要核对已知功能区域，避免把算法能够排出的结果直接当成生物学结论。

　　总结

　　MacVector把序列排列、成对比较和原始数据核查放在同一套工作环境中，便于发现值得继续研究的差异。参数调整的价值，在于帮助判断某个位置的排列是否稳定，而非让错配数量看起来更少。对于研究记录和后续交流，说明材料范围、保留判断依据，往往比给出一张整齐的比对图更有用。希望这些经验能帮助大家读懂比对结果，减少对差异的误判。如有MacVector产品选型或版本配置方面的需求，欢迎联系咨询。