我正在尝试运行使用该函数的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);
我会很感激你的帮助。先感谢您。