matlab工具箱中的生物信息学工具箱(bioinformatics toolbox)功能还算全面,虽然没有什么突出的优点
。包含了:蛋白和核酸分析,系统发育分析以及基因芯片分析等功能。
1、序列分析:以下,带有>>的为命令行!!
例如,我们要查询
RefSeq的一条序列NM_000799,
>> a =
getgenbank('NM_000799','ToFile','testbank');
%利用getgenbank函数,可以设置所需信息的各种参数。
%getgenbank('accession',TOFILE',FILENAME):可以将信息写入文件;
%getgenbank('accession','SEQUENCEONLY',true);可以只提取序列信息;
%getgenbank('accession','PARTIALSEQ',SEQPARAMS);可以只提取一部分序列,SEQPARAMS为[N,M];
%getgenbank('accession','fileformat',fmt);fmt可以为genbank或者fasta。
>>ntseq =
a.Sequence %这个是结构数组,就是在数组变量名后面加一个圆点,例如array.substructure
>>ntdensity(ntseq)
%图示核酸序列各种核苷酸含量&#