《emboss使用》PPT课件.ppt

《《emboss使用》PPT课件.ppt》由会员分享，可在线阅读，更多相关《《emboss使用》PPT课件.ppt（29页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、EMBOSS软件包方林2006-03-06EMBOSS的特点EuropeanMolecularBiologyOpenSourceSuite软件的分类软件数量统一化输入输出扩展方便图形接口丰富EMBOSS软件分类比对序列编辑结果显示酶动力学分析序列特征分析核酸分析蛋白分析系统发育分析信息分析实用工具EMBOSS的输入序列格式abi,acedb,clustaln,codata,cbid,embl/em,experiment,fasta/pearson,gcg/gcg8,genbank/gb/ddbj,gff,hennig86,ig,jackknifer,jackknifernon,mega,meg

2、anon,msf,nbrf,ncbi,nexus/paup,nexusnon/paupnon,pfam/stockholm,phylip/phylipnon,raw,selex,staden,strider,swissprot/swiss/sw,text,treecon,asis共32种EMBOSS的输出文件格式acedb,asn1,clustal/aln,codata,debug,embl/em,fasta/pearson,fitch,gcg/gcg8,genbank/gb,gff,hennig86,ig,jackknifer,jackknifernon,mega,meganon,msf,n

3、brf,ncbi,nexus/paup,nexusnon/paupnon,phylip,phyplipnon,selex,staden,strider,swiss/sw,text/plain/raw,treecon共30种USA统一序列地址format:file:entrydbname:entrylist几种常用的序列格式FASTAEMBLGenBankSwissprotClustalALNFASTA格式gi|8392890|ref|NM_000477.3|Homo sapiens albumin(ALB),mRNAAGCTTTTCTCTTCTGTCAACCCCACACGCCTTTGGCACA

4、ATGAAGTGGGTAACCTTTATTTCCCTTCTTTTTCTCTTTAGCTCGGCTTATTCCAGGGGTGTGTTTCGTCGAGATGCACACAAGAGTGAGGTTGCTCATCGGTTTAAAGATTTGGGAGAAGAAAATTTCAAAGCCTTGGTGTTGATTGCCTTTGCTCAGTATCTTCAGCAGTGTCCATTTGAAGATCATGTAAAATTAGTGAATGAAGTAACTGAATTTGCAAAAACATGTGTTGCTGATGAGTCAGCTGAAAATTGTGACAAATCACTTCATACCCTTTTTGGAGACAAATTATGCACA

5、GTTGCAACTCTTCGTGAAACCTATGGTGAAATGGCTGACTGCTGTGCAAAACAAGAACCTGAGAGAAATGAATGCTTCTTGCAACACAAAGATGACAACCCAAACCTCCCCCGATTGGTGAGACCAGAGGTTGATGTGATGTGCACTGCTTTTCATGACAATGAAGAGACATTTTTGAAAAAATACTTATATGAAATTGCCAGAAGACATCCTTACTTTTATGCCCCGGAACTCCTTTTCTTTGCTAAAAGGTATAAAGCTGCTTTTACAGAATGTTGCCAAGCTGCTGATEMBL格式Gen

6、Bank格式Swissprot格式ClustALN常用的一些EMBOSS程序wossnametfmseqretrevseqtranseqgetorfcoderetbacktranseqpepstatswossname-search要找的关键字-outfilef输出结果到所给的文件-groups只输出程序所在的组的名字-alphabetic只输出程序名和简单的描述，而不显示程序所在组的名字-noemboss不要在EMBOSS主程序库中搜索-noembassay不要在EMBOSS的附加包中搜索-colon用冒号分隔程序所在的类和父类erestmlRestrictionsiteMaximumLike

7、lihoodmethodfrestbootBootstrappedrestrictionsitesalgorithmfrestdistDistancematrixfromrestrictionsitesorfragmentsfrestmlRestrictionsitemaximumLikelihoodmethodrecoderRemoverestrictionsitesbutmaintainsametranslationremapDisplaysequencewithrestrictionsites,translationetc.tfmtfm是另一个有用的程序，一般配合wossname使用用，

8、它只能显示给定程序的帮助信息。但它显示的帮助信息较为详细，不像通过-help或-v那样只显示参数信息。用法：tfmoptionsprogramseqret使用平率最高的一个程序序列的格式转换特征提取序列的拆分序列操作等seqret主要参数-feature用序列的特征信息-firstonly只输出第一条序列-sbegin1n从n位置开始截取序列-send1n截取到n位置-sreverse1对核酸序列反向取补-sask1提问截取序列的起始，终止和是否反向取补-snucleotide1只对核酸进行操作-sprotein1只对蛋白进行操作-slower1序列都表示成小写的形式-supper1序列表示成

9、大写的形式-sformat1s输入文件的格式-sdbname1s数据库名-sid1s序列名-ufo1s统一特征名-fformat1s特征的格式-fopenfile1f特征文件名-osformat2s输出文件格式-ossingle2s将序列分割成一条序列一个文件的形式-ofdirectory2f输出文件的目录-auto关掉提示-stdout输出到标准输出-filter从标准输入读，输出到标准输出getorf用于预测序列的开放阅读框常用参数-table用于翻译的表密码子表-minsizen最小的ORF-maxsizen最大的ORF-nomethionine起始密码子是否编码甲硫氨酸-circula

10、r序列为环状的-noreverse在或不在反链中预测ORF-flankingn追加n长度的侧翼序列revseq用于对序列进行翻转和取补常用参数-noreverse是否取反-nocomplement是否取补transeq用于蛋白的翻译常用的参数-frame翻译那几个相位，对于正链可以为1,2,3，对于反链，为-1,-2,-3-table密码子表，可以为上面的任何一个-regions要翻译的范围，可以表示为：24-45,56-781:45,67=99;765.8881,5,8,10,23,45,57,99数字可以用任何非数字，非字母的符号表示。-trim是否不显示结束密码子，缺省为*表示trans

11、eq支持的密码子表0(Standard)1(Standard(with alternative initiation codons)2(Vertebrate Mitochondrial)3(Yeast Mitochondrial)4(Mold,Protozoan,Coelenterate Mitochondrial andMycoplasma/Spiroplasma)5(Invertebrate Mitochondrial)6(Ciliate Macronuclear and Dasycladacean)9(Echinoderm Mitochondrial)10(Euplotid Nuclea

12、r)11(Bacterial)12(Alternative Yeast Nuclear)13(Ascidian Mitochondrial)14(Flatworm Mitochondrial)15(Blepharisma Macronuclear)16(Chlorophycean Mitochondrial)21(Trematode Mitochondrial)22(Scenedesmus obliquus)23(Thraustochytrium Mitochondrial)coderet用于提取序列的CDS，mRNA或蛋白序列常用的参数有：-nocds不提取CDS序列-nomrna不提取mR

13、NA序列-notranslation不提取translation序列将蛋白反向翻译成核酸常用参数-cfilef密码子频率文件，格式如下#Codon AA Fraction Frequency NumberGCA A 0.088 10.828 28723GCC A 0.423 52.293 138720GCG A 0.428 52.999 140592GCT A 0.062 7.615 20201TGC C 0.854 7.853 20832TGT C 0.146 1.344 3564GAC D 0.665 36.225 96097GAT D 0.335 18.259 48437GAA E 0.

14、299 15.643 41497GAG E 0.701 36.750 97488TTC F 0.825 30.677 81378TTT F 0.175 6.529 17320GGA G 0.085 7.084 18792GGC G 0.700 58.228 154464GGG G 0.115 9.556 25350.backtranseqpepstats计算蛋白的分子量氨基酸数氨基酸使用频率电荷等电点物理化学性质pepstats常用参数和类似程序常用参数：-notermini是否包含蛋白的N和C端的电荷-aadatas氨基酸分子量信息的文件，缺省下用标准的氨基酸分子量类似程序有：pepinfo，iep，charge和freak等EMBOSS的用户接口JembossEmnuEMBOSS-GUIJembossEMNUEMBOSS-GUI