我正在尝试运行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);
感谢您的帮助。提前谢谢。