資源簡介
這段程序的主要作用是,從一個基因庫中,查找到我們的目的基因,然后將該目的基因連帶其后面的核苷酸序列一起按照fasta格式輸出。為了方便后面的靶基因的預(yù)測,可以減少很大的工作量。本人從事生物信息分析很多年,可以幫忙進行一些個性化的數(shù)據(jù)分析工作。有興趣聯(lián)系我吧

代碼片段和文件信息
clear?all
fid=fopen(‘a(chǎn)ll0.txt‘‘r‘);
fid1=fopen(‘result_selete.txt‘‘wt+‘);?%查找
fid2=fopen(‘target.txt‘‘r‘);
?tline2=fgetl(fid2);
??tline=fgetl(fid);
?while?~feof(fid)?
?????tline0=tline;
?????tline=fgetl(fid);
?????a=size(findstr(tline2tline0(2:15)));
?????if?a(11)==1&&tline0(16)==‘.‘
????????fprintf(fid1‘%s\n‘tline0);
????????fprintf(fid1‘%s\n‘tline);
?????end
??????
?end
?fclose?all
??clear?all
?fid2=fopen(‘result_selete.txt‘‘r‘);
???fid1=fopen(‘部分查找出的序列.txt‘‘wt‘);?%分行顯示查找結(jié)果
???while?~feof(fid2)
????????tline2=fgetl(fid2);
??if?tline2(1)==‘>‘
????????fprintf(fid1‘%s\n‘tline2);?
??end?
??if?tline2(1)~=‘>‘
?????len=length(tline2);
?????e=mod(len60);
?????k=(len-e)/60;???????
????for?i=1:k
??????fprintf(fid1‘%s\n‘tline2(60*(i-1)+1:60*i));?
????end
???if?e~=0
?????fprintf(fid1‘%s\n‘tline2(60*k+1:len));
???end
??end
???end
???
?fclose?all??????????????
???????????
?
?
?
?
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????967??2014-02-25?20:59??基因序列的篩選\select.m
?????文件????????146??2014-04-08?16:40??基因序列的篩選\使用說明.txt
?????文件??????????0??2014-04-08?16:37??基因序列的篩選\使用說明.txt.bak
?????目錄??????????0??2014-04-08?16:40??基因序列的篩選
-----------?---------??----------?-----??----
?????????????????1113????????????????????4
評論
共有 條評論