《emboss使用》PPT课件

EMBOSS软件包,方林2006-03-06,EMBOSS的特点,EuropeanMolecularBiologyOpenSourceSuite软件的分类软件数量统一化输入输出扩展方便图形接口丰富,EMBOSS软件分类,比对序列编辑结果显示酶动力学分析序列特征分析核酸分析蛋白分析系统发育分析信息分析实用工具,EMBOSS的输入序列格式,abi,acedb,clustaln,codata,cbid,embl/em,experiment,fasta/pearson,gcg/gcg8,genbank/gb/ddbj,gff,hennig86,ig,jackknifer,jackknifernon,mega,meganon,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,nbrf,ncbi,nexus/paup,nexusnon/paupnon,phylip,phyplipnon,selex,staden,strider,swiss/sw,text/plain/raw,treecon共30种,USA统一序列地址,format:file:entrydbname:entrylist,几种常用的序列格式,FASTAEMBLGenBankSwissprotClustalALN,FASTA格式,gi|8392890|ref|NM_000477.3|Homosapiensalbumin(ALB),mRNAAGCTTTTCTCTTCTGTCAACCCCACACGCCTTTGGCACAATGAAGTGGGTAACCTTTATTTCCCTTCTTTTTCTCTTTAGCTCGGCTTATTCCAGGGGTGTGTTTCGTCGAGATGCACACAAGAGTGAGGTTGCTCATCGGTTTAAAGATTTGGGAGAAGAAAATTTCAAAGCCTTGGTGTTGATTGCCTTTGCTCAGTATCTTCAGCAGTGTCCATTTGAAGATCATGTAAAATTAGTGAATGAAGTAACTGAATTTGCAAAAACATGTGTTGCTGATGAGTCAGCTGAAAATTGTGACAAATCACTTCATACCCTTTTTGGAGACAAATTATGCACAGTTGCAACTCTTCGTGAAACCTATGGTGAAATGGCTGACTGCTGTGCAAAACAAGAACCTGAGAGAAATGAATGCTTCTTGCAACACAAAGATGACAACCCAAACCTCCCCCGATTGGTGAGACCAGAGGTTGATGTGATGTGCACTGCTTTTCATGACAATGAAGAGACATTTTTGAAAAAATACTTATATGAAATTGCCAGAAGACATCCTTACTTTTATGCCCCGGAACTCCTTTTCTTTGCTAAAAGGTATAAAGCTGCTTTTACAGAATGTTGCCAAGCTGCTGAT,EMBL格式,GenBank格式,Swissprot格式,ClustALN,常用的一些EMBOSS程序,wossnametfmseqretrevseqtranseqgetorfcoderetbacktranseqpepstats,wossname,-search要找的关键字-outfilef输出结果到所给的文件-groups只输出程序所在的组的名字-alphabetic只输出程序名和简单的描述，而不显示程序所在组的名字-noemboss不要在EMBOSS主程序库中搜索-noembassay不要在EMBOSS的附加包中搜索-colon用冒号分隔程序所在的类和父类,erestmlRestrictionsiteMaximumLikelihoodmethodfrestbootBootstrappedrestrictionsitesalgorithmfrestdistDistancematrixfromrestrictionsitesorfragmentsfrestmlRestrictionsitemaximumLikelihoodmethodrecoderRemoverestrictionsitesbutmaintainsametranslationremapDisplaysequencewithrestrictionsites,translationetc.,tfm,tfm是另一个有用的程序，一般配合wossname使用用，它只能显示给定程序的帮助信息。但它显示的帮助信息较为详细，不像通过-help或-v那样只显示参数信息。用法：tfmoptionsprogram,seqret,使用平率最高的一个程序序列的格式转换特征提取序列的拆分序列操作等,seqret主要参数,-feature用序列的特征信息-firstonly只输出第一条序列-sbegin1n从n位置开始截取序列-send1n截取到n位置-sreverse1对核酸序列反向取补-sask1提问截取序列的起始，终止和是否反向取补-snucleotide1只对核酸进行操作-sprotein1只对蛋白进行操作-slower1序列都表示成小写的形式-supper1序列表示成大写的形式-sformat1s输入文件的格式-sdbname1s数据库名-sid1s序列名-ufo1s统一特征名-fformat1s特征的格式-fopenfile1f特征文件名-osformat2s输出文件格式-ossingle2s将序列分割成一条序列一个文件的形式-ofdirectory2f输出文件的目录-auto关掉提示-stdout输出到标准输出-filter从标准输入读，输出到标准输出,getorf,用于预测序列的开放阅读框常用参数-table用于翻译的表密码子表-minsizen最小的ORF-maxsizen最大的ORF-nomethionine起始密码子是否编码甲硫氨酸-circular序列为环状的-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是否不显示结束密码子，缺省为*表示,transeq支持的密码子表,0(Standard)1(Standard(withalternativeinitiationcodons)2(VertebrateMitochondrial)3(YeastMitochondrial)4(Mold,Protozoan,CoelenterateMitochondrialandMycoplasma/Spiroplasma)5(InvertebrateMitochondrial)6(CiliateMacronuclearandDasycladacean)9(EchinodermMitochondrial)10(EuplotidNuclear)11(Bacterial)12(AlternativeYeastNuclear)13(AscidianMitochondrial)14(FlatwormMitochondrial)15(BlepharismaMacronuclear)16(ChlorophyceanMitochondrial)21(TrematodeMitochondrial)22(Scenedesmusobliquus)23(ThraustochytriumMitochondrial),coderet,用于提取序列的CDS，mRNA或蛋白序列常用的参数有：-nocds不提取CDS序列-nomrna不提取mRNA序列-notranslation不提取translation序列,将蛋白反向翻译成核酸常用参数-cfilef密码子频率文件，格式如下#CodonAAFractionFrequencyNumberGCAA0.08810.82828723GCCA0.42352.293138720GCGA0.42852.999140592GCTA0.0627.61520201TGCC0.8547.85320832TGTC0.1461.3443564GACD0.66536.22596097GATD0.33518.25948437GAAE0.29915.64341497GAGE0.70136.75097488TTCF0.82530.67781378TTTF0.1756.52917320GGAG0.0857.08418792GGCG0.70058.228154464GGGG0.1159.55625350.,backtranseq