替代matlab中的fastaread函数

时间:2012-09-25 19:31:04

标签: matlab bioinformatics

我正在尝试运行matlab代码,该代码使用我缺乏的生物信息学工具箱中提供的函数fastaread。我想知道我是否可以创建自己的函数,从.fasta文件函数读取数据以运行此代码。

我很抱歉这种非常具体的天真问题。但任何帮助将非常感激。先感谢您。


感谢您的回复。我将readfasta.m从mbetoolbox重命名为fastaread.m。但是我收到以下错误:

??? Error: File: /home/users/fastaread.m Line: 85 Column: 1
Expression or statement is incorrect--possibly unbalanced (, {, or [.

Error in ==> calculate_cons>read_alignment_fasta at 111
full_alignment = fastaread(msa_fasta_filename);

Error in ==> calculate_cons>read_alignment at 67
[encoded_focus_alignment, focus_index_of_interest, focus_to_uniprot_offset_map]
= read_alignment_fasta(msa_fasta_filename, seqid_of_interest);

Error in ==> calculate_cons at 41
[Pij_true, Pi_true, alignment_width, q, encoded_seq_of_interest,
focus_to_uniprot_offset_map] = read_alignment(msa_fasta_filename,
seqid_of_interest, theta);

感谢您的帮助。提前谢谢。

1 个答案:

答案 0 :(得分:1)

只需使用mbetoolbox中的readfasta功能即可。