最经典的Bioedit使用说明书课件

综合序列分析软件BioEdit2003级高芳銮BioEdit简介BioEdit是一个性能优良的免费的生物序列编辑器，可在Windows95/98/NT/2000中运行，它的基本功能是提供蛋白质、核酸序列的编辑、排列、处理和分析。与DNAMAN相比，其分析内容相对丰富一些，而且提供了很多网络程序的分析界面和接口，与DNAMAN等软件配合使用更好。尤其值得一提是利用BioEdit能够十分方面地根据指定的核酸序列绘制相应的质粒图谱。序列的常规操作：序列输入：多种序列输入方式；序列分类：按标题、位置、定义、参数、注释等分类；成对排列：两序列的最佳排列及计算同一性和类似性；序列屏蔽：仅采用联配中部分区域进行分析而排除其他。核酸分析：组成、互补、反转、翻译、质粒、限制性内切酶；蛋白质分析：氨基酸成分、疏水性轮廓、疏水力矩平均数翻译或反翻译：把DNA或RNA翻译成蛋白质；切换翻译：在核酸和编码蛋白质序列中切换核苷酸序列；点图[成对比较]：相互比较两序列的矩阵，生成一个点图。主要内容绘制质粒图限制性内切酶图蛋白质分析组成分析熵图疏水性轮廓联配中搜寻保守区根据密码子的使用翻译核苷酸RNA比较分析共变潜在配对互交信息分析一、绘制质粒图（Plasminddrawing）　使用BioEdit质粒绘图功能，序列可以通过自动的位置标记，自动修改成环形质粒。特征、多连接位点和限制性位点可以通过使用对话框增加。当将一个序列进入质粒图时，在背景上出现一个限制性内切酶图谱，所以可以通过对话框选择可以增加限制性位点。它们自动增加到当前的位点。质粒功能提供简单的绘制和标记工具。标签和绘图可以通过鼠标移动和缩放。想要编辑目标性质，双击目标。　　想要从一个DNA序列产生一个质粒，从“Sequence”菜单中“NucleicAcid”子菜单中选择“CreatePlasmidfromSequence”选项。选择这个选项时，限制性内切酶图谱将会使用通常商业化的，储存在存储器中的限制性内切酶。质粒第一次产生时，它显示成有10个位点标记的圆圈，中央是标题。1.Restrictionsites:(限制性位点)　　想要增加限制性位点，从“Vector”菜单中选择“RestrictionSites”选项。将会显示一下对话框：2.Positionalmarks(位置标记):　　点击“Vector”菜单中的“PositionalMarks”选项，可以出现以下对话框： 　　可以通过移动位置标记到“Show”中，单独增加位置标记，或者设定应用的分割标记数量。想要没有标记，选择“Divideinto:”中的下拉菜单顶端的“None”。4.GeneralVectorproperties　　载体属性可通过选“Vector”菜单中的“Properties”来更改：二、RestrictionMaps(限制性内切酶图)　BioEdit提供两种方法产生核苷酸序列的限制性内切酶图。一种内在的限制性内切酶图功能允许产生序列最多为65,536个核苷酸的限制性内切酶图。实际上，只能检测大约35Kb，而且在速度慢的计算机上会要消耗很长的时间。　你也可以通过万维网直接链接到WebCutter限制性内切酶图上。1.WebCutter：点亮你想要图谱的序列标题，从“WorldWideWeb”菜单中选择“Auto-fedWebCutterRestrictionMapping”2.BioEdit:点亮你想要图谱的序列标题，从“Sequence”菜单选择“RestrictionMap”。以下选项将会显示在一个界面窗口：—显示图谱：显示或省略序列的全图谱,互补链显示每个酶的酶切位点.默认值：yes—按照字母顺序排列名称：显示关于所有内切酶、它们的识别序列、切割频率和所有位置(5’末端开始是1)的列.默认值：yes—位置数：关于酶切位点的列表.默认值：no—唯一位点列表：在全部序列中只有一个酶切位点的内切酶列表.默认值：no—切割5次或更少的酶.默认值：yes—频率汇总表：关于所有正确选择的内切酶和它们切割序列的次数。默认值：no—不能切割的内切酶。默认值：yes—4-碱基内切酶：想要包括这些酶,必须点击这个选项.默认值：no(不包括本身)—5-碱基内切酶：与4-basecutters相同.—非严格识别序列的酶：有时你可排除它们.默认值：yes—大的识别位点：通常用于克隆,只有共同的6-碱基识别酶被使用.—同裂酶：若只显示一个特殊识别位点的一个内切酶,不选(默认值=不选择).—翻　　译：显示沿着排列中的序列翻译(5’端到3’端的由左到右的翻译)—互补翻译：互补链的翻译方向相反.—编号方式：是酶切位点的核酸的号码,而不是识别位点的起点.3.RestrictionEnzymeBrowser(限制性内切酶浏览器)　　从核酸序列中得到内切酶谱时，显示酶的生产公司是很有用的。通过在内切酶图谱中选择制造厂商和按下按钮，可以手动浏览内切酶。你也可以通过选择“Options”菜单中的“ViewRestrictionEnzymesbyManufacturer”选择，在任何时候检查内切酶。显示如右对话框：1.氨基酸的组成　　从“Sequence”菜单下进入“Protein”,再进入“aminaacidcomposition”,可对序列的氨基酸组成分析，结果以摘要和图例的形式给出。　　图例中的柱形条表示每种氨基酸在序列中的摩尔比，如下图：三、蛋白质分析以RGDV的minoroutercapsidprotein－AAS66885为例：3.疏水性轮廓(profile)　　平均疏水性轮廓采用Kyte&Doolittle的方法，平均分值(总和/窗口大小)作为序列中各个位置的疏水性值，并以窗口中中间残基的疏水性值作图。4.瞬间疏水性轮廓(hydrophobicmomentprofile)6.在联配中搜寻保守区　　有时，即使序列之间的变化很大时，在几个序列中搜寻保守区是有用的。例如，根据一系列同源序列发现通用的PCR引物。BioEdiot查找的是低平均“熵”的区域。　　首先选择你的序列，从“Aligment”->“FindConservedRegion”，对话框中各选项的内容：BioEditversion5.0.9ConservedregionsearchAlignmentfile:Q:\Ribosomal_RNA\some_methanos.bio5/10/048:57:33PM Minimumsegmentlength(actualforeachsequence):15Maximumaverageentropy:0.2Maximumentropyperposition:0.2Gapslimitedto2persegmentContiguousgapslimitedto1inanysegment 2conservedregionsfound Region1:Position755to774Consensus:755AUUAGAUACCCGGGUAGUCC774SegmentLength:20Averageentropy(Hx):0.0155Position755:0.0000Position756:0.0000Position757:0.0000Position758:0.0708Position759:0.0000Position760:0.0000Position761:0.0000Position762:0.0000Position763:0.0000Position764:0.0708Position765:0.0000Position766:0.1679Position767:0.0000Position768:0.0000Position769:0.0000Position770:0.0000Position771:0.0000Position772:0.0000Position773:0.0000Position774:0.0000BioEditversion5.0.9ConservedregionsearchAlignmentfile:G:\Ribosomal_RNA\some_methanos.bio5/10/999:34:06PM Minimumsegmentlength(actualforeachsequence):10Maximumaverageentropy:0.4Maximumentropyperposition:0.4with2exceptionsallowedGapslimitedto2persegmentContiguousgapslimitedto1inanysegment 36conservedregionsfound结果：7.根据密码子的使用翻译核苷酸　　核苷酸序列可根据三联体密码翻译预测的蛋白序列。　　从“Sequence”->“Protein”->“Translation”,选择要按何种读框翻译。　　例如，以下是一个假设的Methanobacterium(甲烷细菌)的ORF(开放阅读框架)。>MTH671codingregionATGGTTGCAGTACCCGGCAGTGAGATACTGAGCGGTGCACTACACGTTGTCTCCCAGAGCCTCCTCATACCGGTTATAGCAGGTCTACTGTTATTCATGGTATACGCCATAGTGACCCTCGGAGGGCTCATATCAGAGTACTCTGGAAGGATAAGGACTGATGTTAAGGAACTTGAATCGGCAATAAAATCAATTTCAAACCCAGGAACCCCTGAAAAGATAATTGAGGTCGTCGATTCGATGGACATACCACAGAGCCAGAAGGCCGTGCTCACTGATATCGCAGGGACAGCTGAACTCGGACCAAAATCAAGGGAGGCCCTCGCAAGGAAGTTGATAGAGAATGAGGAACTCAGGGCTGCCAAGAGCCTTGAGAAGACAGACATTGTAACCAGACTCGGCCCAACCCTTGGACTGATGGGGACACTCATACCCATGGGTCCAGGACTCGCAGCCCTCGGGGCAGGTGACATCAATACACTGGCCCAGGCCATCATCATAGCCTTCGATACAACAGTTGTGGGACTTGCATCAGGGGGTATAGCATACATCATCTCCAAGGTCAGGAGAAGATGGTATGAGGAGTACCTCTCAAATCTTGAGACAATGGCCGAGGCAGTGCTGGAGGTGATGGATAATGCCACTCAGACGCCGGCGAAGGCTCCTCTCGGATCAAAAAframe1ofthissequenceisdisplayedasfollowsintheBioEdittexteditor:>MTH671codingregion 1ATGGTTGCAGTACCCGGCAGTGAGATACTGAGCGGTGCACTACAC451MetValAlaValProGlySerGluIleLeuSerGlyAlaLeuHis15 46GTTGTCTCCCAGAGCCTCCTCATACCGGTTATAGCAGGTCTACTG9016ValValSerGlnSerLeuLeuIleProValIleAlaGlyLeuLeu30 91TTATTCATGGTATACGCCATAGTGACCCTCGGAGGGCTCATATCA13531LeuPheMetValTyrAlaIleValThrLeuGlyGlyLeuIleSer45 136GAGTACTCTGGAAGGATAAGGACTGATGTTAAGGAACTTGAATCG18046GluTyrSerGlyArgIleArgThrAspValLysGluLeuGluSer60 181GCAATAAAATCAATTTCAAACCCAGGAACCCCTGAAAAGATAATT22561AlaIleLysSerIleSerAsnProGlyThrProGluLysIleIle75 226GAGGTCGTCGATTCGATGGACATACCACAGAGCCAGAAGGCCGTG27076GluValValAspSerMetAspIleProGlnSerGlnLysAlaVal90 …|ACGT|-----------------------------A|37313|A|0.760.120.040.07||LysThrArgIle|-----------------------------A|1446|C|0.610.430.270.46||AsnThrSerIle|-----------------------------A|8167|G|0.240.230.031||LysThrArgMet|-----------------------------A|4313|T|0.390.210.130.47||AsnThrSerIle|-----------------------------………四、RNA的比较分析　　RNA的结构定义为核苷酸的碱基的相互作用。最简单情况下，即螺旋中的碱基对之间的Waltson-Crick碱基配对。RNA结构的系统发育比较分析方法建立在如下假定上，即在进化中核苷酸改变，但重要的RNA二级和三级结构保持不变。一个可能破坏结构的碱基变化可以由序列中另一处的变化补偿以保持结构稳定。所以不同物种的同源RNA中将包含“补偿碱基变化”或“共变化，协变（covariation）”。所以通过检查来自各个不同生物的同源RNA，确定这些“补偿碱基变化”，从而阐明结构。　　例如，一给定的序列，GAAGA将可能与序列中任一UCUUC配对，而后者可能在序列中出现数次。如何确定到底是和哪一个配对呢？可以检查不同生物的同源RNA序列，找出“补偿碱基变化”。organism#1‑‑‑‑‑GAAGA‑‑‑‑‑‑‑‑‑UCUUC‑‑‑‑‑‑‑‑UCUUC‑‑‑‑‑‑‑‑‑UCUUC‑‑‑‑‑‑‑organism#2‑‑‑‑‑GAUGA‑‑‑‑‑‑‑‑‑UCUUC‑‑‑‑‑‑‑‑UCUGC‑‑‑‑‑‑‑‑‑UCAUC‑‑‑‑‑‑‑organism#2‑‑‑‑‑GAUGA‑‑‑‑‑‑‑‑‑GCUUC‑‑‑‑‑‑‑‑UCUAC‑‑‑‑‑‑‑‑‑UCAUC‑‑‑‑‑‑‑organism#2‑‑‑‑‑GACGA‑‑‑‑‑‑‑‑‑UCUUC‑‑‑‑‑‑‑‑UCUGC‑‑‑‑‑‑‑‑‑UCGUC‑‑‑‑‑‑‑　　在此例中，只有最后一个UCUUC才可和GAAGA配对。象这样在序列中2个位置出现“补偿碱基变化”，被认为是螺旋存在的证据。两条序列不能形成互补，表明不存在配对。在“系统发育比较分析”中关键是序列联配，同源序列必须适当联配。此处同源性是严格意义的：同源的核苷酸来自一个共同的祖先。所以开始时，先使用关系紧密的序列进行联配，这样在序列相似性基础上联配，不需要加入许多联配的空位。联配后互补序列的“协变”可被立即发现，从而开始构建二级结构,然后差异大的序列可以添进联配中。这样持续添加新序列，进行“协变”分析，直到联配和二级结构模型出现此过程的完全描述。一旦一个完整的二级结构模型形成，“协变”分析可以鉴定非螺旋区的核苷酸之间的相互作用以及不规则的相互作用。之所以可以被鉴定,是因为涉及的核苷酸即使不形成规则的碱基配对或是一个螺旋的一部分，也仍一致的变化。1.共变化(Covariation)　　共变化指序列中两个残基步调一致地变化。严格地讲即每当联配序列中x变化时，y也变化，两者是一致的。（例如，当x变为A，y变为T。每次x变为A，y一定变为T）。残基间的共变化表明，它们之间一定有重要的相互作用，当重要结构残基突变时，自然选择保留了那些有补偿突变的序列。共变化的例子　　假设我们现有一个联配序列，它表示了几种物种共有的一个特定的RNA的保守的结构。我们希望从联配中包含的信息推测出RNA二级结构。　　　　　　　....|....|....|....|....|....1020sample1CCGGAUACGAUCGUCGGGUACGUAUCCGGsample2CCGGAUACUAUCUUGGCGAAAGUAUCUGGsample3CGGGAUACGAUCGACGCGUACGUAUCCCGsample4CGCGGUACCAUCCACCCCUAGGUACCGCGsample5CCGGAUACGAUCGUCCCGUUCGUAUCCGGsample6CCGGAUACGAUCGUCGGGUACGUAUCCGGsample7CCGGACACGAUCGUCGGGUACGUAUCCGGsample8CCAGAUACGAUCGAAACUUUCGUAUCUGGsample9CCGGUUACCAUCGUCGGGUAGGUAACCGGsample9CCGGAUACGAUCGACAGGAACGUAUCCGGsample10CCGGAUACGAUCGUCCCGUACGUAUCCGGsample11CCGGAUACGAUCGUCGGGUACGUAUCCGGsample12CCUGAUACUAUCGUCGCCUAAGUAUCGGGsample13CGGGGUACGAUCGAGGCCUACGUACCCCGsample14CCCGCUACGAUCGAGGCCUUCGUAGCGGGsample15CCGGAUACGAUCGAGGCCUUCGUAUCCGG下面是一个联配的例子CovariationanalysisInputfile:I:\BioEdit\help\samples.gbPositionnumberingisrelativetothealignmentnumbering.Nomaskwasused. 1CCCCCCCCCCCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ Position2:‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑2CCGGCCCCCCCCCGCC28GGCCGGGGGGGGGCGGAllpotentialWatsonCrickorG‑Upairs‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑3GGGCGGGAGGGGUGCG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑4GGGGGGGGGGGGGGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ Position5:‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑5AAAGAAAAUAAAAGCA25UUUCUUUUAUUUUCGUAllpotentialWatsonCrickorG‑Upairs‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑6UUUUUUCUUUUUUUUU‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑7AAAAAAAAAAAAAAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑8CCCCCCCCCCCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position9:‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑9GUGCGGGGCGGGUGGG21CACGCCCCGCCCACCCAllpotentialWatsonCrickorG‑Upairs‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑10AAAAAAAAAAAAAAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑11UUUUUUUUUUUUUUUU‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑12CCCCCCCCCCCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑13GUGCGGGGGGGGGGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑14UUAAUUUAUAUUUAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑15CGCCCCCACCCCCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑16GGGCCGGAGACGGGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑17GCCCCGGCGGCGCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑18GGGCGGGUGGGGCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑19UAUUUUUUUAUUUUUU‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑20AAAAUAAUAAAAAAUU‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position21:‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑21CACGCCCCGCCCACCC9GUGCGGGGCGGGUGGGAllpotentialWatsonCrickorG‑Upairs‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑22GGGGGGGGGGGGGGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑23UUUUUUUUUUUUUUUU‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑24AAAAAAAAAAAAAAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ Position25:‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑25UUUCUUUUAUUUUCGU5AAAGAAAAUAAAAGCAAllpotentialWatsonCrickorG‑Upairs‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑26CCCCCCCCCCCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑27CUCGCCCUCCCCGCGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position28:‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑28GGCCGGGGGGGGGCGG2CCGGCCCCCCCCCGCCAllpotentialWatsonCrickorG‑Upairs‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑29GGGGGGGGGGGGGGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑　　在上述联配中共有3对“共变化”的位置点：2/28，5/25，9/21。两个碱基共变表明它们很可能相互作用。如果一个突变发生在与其他碱基有重要作用的碱基上（常是碱基对），选择压力可能会只保留在另一处碱基上发生补偿突变的碱基。事实上，上述的碱基共变化都发生在规则的碱基对（Watson-Crick碱基对或在RNA中G-U）表明它们可能是碱基配对。共变化碱基对2/5分别和5/25的距离相同，而5/25分别和9/21的距离也相同，而且界于它们之间的碱基也可形成碱基互补，这都表明联配序列的两端可能闭合形成螺旋如下是“Sample1”形成的结构。　　　　　　　　　　UCAG--CCGGATACGU--GGCCTATGCCAGUGG2.潜在配对分析potentialpairing　　当RNA分子中两个核苷酸之间存在配对碱基的相互作用力。一个碱基发生突变，另一个碱基为了补偿这一突变，可能不仅仅是某一特定核苷酸突变（例如原来的A-T配对可能在一序列中转换为G-C,而另一序列中为G-U,）这在共变化分析中将被忽略。因为此种改变并不遵循完全相同的模式。要鉴定这种情况，可以在潜在配对中选定碱基配对的规则。　　仍用上例中的序列（sample1－sample15略）　　BioEdit中并不要求有位置变化，所以未改变的位置上只要可以形成碱基对，也能被发现同时也可在“preference”中设置以滤出未改变的位置之间的碱基配对。以下是一个联配序列它和在共变化分析中使用的相同。设置允许A-U/G-C/G-U碱基配对规则以及1个错配，产生下列的结果（以清单格式，滤除了未变化位置的潜在配对）比较这一结果和共变化的结果，发现位置3/27有一潜在的配对，而共变化的结果未检出。潜在配对的数据也可以按允许的配对出现的频率或原始允许配对的数目列出一个（二维矩阵）表。PotentialPairingsListInputFile:I:\BioEdit\help\samples.gbAllowedMispairings=116totalsequences,29nucleotidespersequence.Axesreflectnumberingoftheentirealignment.NoMaskwasused.Hitsoninvariantpairshavebeenfilteredout. ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑1CCCCCCCCCCCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position:2‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑2CCGGCCCCCCCCCGCC28GGCCGGGGGGGGGCGG0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position:3‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑3GGGCGGGAGGGGUGCG27CUCGCCCUCCCCGCGC0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ Position:4‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑4GGGGGGGGGGGGGGGG6UUUUUUCUUUUUUUUU0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ Position:5‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑5AAAGAAAAUAAAAGCA25UUUCUUUUAUUUUCGU0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position:6‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑6UUUUUUCUUUUUUUUU4GGGGGGGGGGGGGGGG0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑6UUUUUUCUUUUUUUUU7AAAAAAAAAAAAAAAA1mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑6UUUUUUCUUUUUUUUU10AAAAAAAAAAAAAAAA1mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑6UUUUUUCUUUUUUUUU22GGGGGGGGGGGGGGGG0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑6UUUUUUCUUUUUUUUU24AAAAAAAAAAAAAAAA1mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑6UUUUUUCUUUUUUUUU29GGGGGGGGGGGGGGGG0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position:7‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑7AAAAAAAAAAAAAAAA6UUUUUUCUUUUUUUUU1mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ 8CCCCCCCCCCCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position:9‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑9GUGCGGGGCGGGUGGG21CACGCCCCGCCCACCC0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ Position:10‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑10AAAAAAAAAAAAAAAA6UUUUUUCUUUUUUUUU1mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑11UUUUUUUUUUUUUUUU‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑12CCCCCCCCCCCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑13GUGCGGGGGGGGGGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑14UUAAUUUAUAUUUAAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑15CGCCCCCACCCCCGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑16GGGCCGGAGACGGGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑17GCCCCGGCGGCGCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑18GGGCGGGUGGGGCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑19UAUUUUUUUAUUUUUU‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑20AAAAUAAUAAAAAAUU‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position:22‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑22GGGGGGGGGGGGGGGG6UUUUUUCUUUUUUUUU0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ 23UUUUUUUUUUUUUUUU‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position:24‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑24AAAAAAAAAAAAAAAA6UUUUUUCUUUUUUUUU1mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ Position:25‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑25UUUCUUUUAUUUUCGU5AAAGAAAAUAAAAGCA0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ 26CCCCCCCCCCCCCCCC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position:27‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑27CUCGCCCUCCCCGCGC3GGGCGGGAGGGGUGCG0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑Position:28‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑28GGCCGGGGGGGGGCGG2CCGGCCCCCCCCCGCC0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ Position:29‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑29GGGGGGGGGGGGGGGG6UUUUUUCUUUUUUUUU0mis‑matches‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑3.交互信息分析（MutualInformationAnalysis）概述　　交互信息，象在系统发育比较分析中的应用一样，主要是衡量在一个适当联配中两个位置共有信息的信息量。符号是M(x,y)(位置x,y的相互信息)。M(x,y)表明两个位置相关的紧密程度。此相关程度显示了两位置的直接相互作用，如碱基配对。BioEdit另外计算R1和R2两个参数，它们分别表示位置x,y对M(x,y)的贡献。什么是交互信息　　交互信息分析是以下思想的拓展--即对某个特定位置的不确定性表示是信息含量的下降。在预先对某位置一无所知的情况下（如RNA中核苷酸），不确定性最大。但一旦确定了某位置是什么核苷酸时,不确定性消除了,此位置的信息量达到最大。现在考虑有多条序列，在某位置均含有一个同源核苷酸。知道第一条序列上此位置上的核苷酸并不能为确定第二条或随机的一条序列中此序列的核苷酸提供多少信息。但是如果已知此位置在许多乃至几乎所有序列中均为某一特定碱基（如C），而不是其它的碱基（如G），则我们积累了相当多的“信息”，可预测另一个未检测的序列中，在此位置某核苷酸出现的可能性。即在另一未检测的序列中，此位置核苷酸的不确定性下降了。　　交互信息进一步拓展了这一思想，对配对位置的信息量进行检查，此信息量依赖于并联系每个位置单独的信息量，但不能将两者混淆。总的讲，它衡量不确定性的下降，此不确定性指两种事物相互影响相互作用的程度。RobinGutell发展了用交互信息预测RNA结构的方法，也很适合系统发育比较分析，因为两个位置交互信息高也提示这2个残基直接相互作用。1234ACGUACGUAGCUAUAUAUAUAAUUAAUUAGCU如左图总共8个序列，其中位置1，4是不改变的,信息量最大。位置2，3中C/G/U/A各出现了2次，信息量为0，我们无法预测下一个序列中这两个位置的核苷酸，但位置2，3都含有它们之间是如何影响彼此的共有信息。我们不能猜出新一序列中位置2的核苷酸，但如果告诉我位置3是C，我们可以推断出位置2是G，这即建立在“交互信息分析”（它们遵循共同的配对模式）交互信息也表明这些碱基可能相互作用。交互信息示例　　以下是分析细菌RNasePRNA的部分序列的一个例子。　　点击（Aligment）可以观察此联配。设置输出是全部列表（fulltable）显示M(x,y)的数值。Nbest列出各个位置5个高分值。序列和编号mask都是根据E.coli.。序列的编号是根据E.coli的mask序列。此序列中包含了一个RNasePRNA结构区域的“cruciformregion”(十字型区域)。由于矩阵文件太大，不能在此说明文件中打开观察。但可通过打点作图方便地观察。在BioEdit矩阵作图程序中，数据既可以数字也可以图形的方式被动态的检查。　　其中交互信息分析的“cruciformregion”(十字型区域)在此输出中是环型的。此图象及全部最新的细菌和古细菌的RNasePRNA结构和序列均可在RNaseP数据库中找到。交互信息的作图示例　细菌RNasePRNA联配的一部分,共有包含极丰富信息的138条序列。序列包含“cruciformregion”(十字型区域)。[具体数据详见文本]。使用矩阵打点作图器分析交互信息数据　　例如使用E.coli作为序列和编号masks，对146条细菌RNasePRNA序列进行M(x,y)分析得到的矩阵进行使用矩阵打点作图器分析交互信息数据。　　作图完成后打开“数据检查器”，用鼠标在图中各点移动观察各数据点，可直接用点击某点在顶部的工具条中将出现数据值。　　通过设置数据点的阈值在矩阵作图中遮蔽（shading）某些数据点，当仅需要显示出高分区域时，此选择可能有用。　　下面的细菌RNasePRNA联配的交互信息数据作图（E.coli作为mask）。在图中很难挑出某一核苷酸三联体，虽然碱基配对的位置94和104是明显可见的，但很难挑出配对核苷酸94-104及第三个核苷酸316组成的三联体。　　此图是联配的全部M(x,y)列式的部分图示，设置固定数据点大小为3*3象素。鼠标箭头指向位置94-104，右边的小红色框中心是位置94-316。　此图中位置316和位置94，104的相互作用并不明显，从作图窗口下“Plot”“LineGraphofRows”进入一维行向作图。　　下图显示第316行的作图用“Row”旁边的上下箭头选择要观察的行或直接在框中输入要观察的行。可移动图中的蓝十字点击任何位置，将在顶部工具条的左上方列出位置x,y及数据值。其中数据值是指图中位置的高峰对应的数值。交互信息检查器(mutualinformationexaminer)　　如果希望在联配窗口中观察任意两个位置的交互信息，从“View”菜单下进入“mutualinformationexaminer”，下面是控制条的格式。　　在“x”,”y”旁边输入要分析的位置。如果希望某一特定序列上此位置的信息（序列需无“gap”）,就把此序列设为“numberingmask”,再输入x、y的位置。图中的x、y位置（x=261,y=289）对应的是前面在矩阵作图时选择的数据点。一定要将需要分析的序列全部选中。最后点击“calculate”,将出现以下窗口。　　在上图窗口中点击“Text”,将出现如下文本编辑窗口，内容可复制和粘贴。　　如果采用此位置输入是X=a-b,Y=c-d,同时计算几对位置，BioEdit假设你要分析是螺旋区域，即位置是反向平行的，所以c-d和d-c是不区分的。如果你要特别地指明某对的顺序，要采用X=a,b,c,dY=e,f,g,h.的形式的输出。-------------------------------------------------------------------附：BioEditV5.09：官方下载红意下载　　　　　　　汉化补丁BioEditV6.0.5:　分卷1　分卷2谢谢！