生物信息学编程题目

《生物信息学编程题目》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、生物信息学编程题目1,编程实现将DNA序列文件exam.fa中，去掉每段序列的第一行（＞及其后面的信息以及回车苻信息），将序列信息合并成基因组后统计该基因组的A、C、T、G的个数，并生成把A、C、T、G统计信息写入文件取名stats.fa；同时把拼成基因组写入文件取名genome.fa2,编程实现读入上面完成基因组文件genome.fa,统计酶切位点AC（N5）CTCC个数,并切割模式（N8）AC（N5）CTCC（N6）的片段生成例1所示的Fasta文件Tag.fa,即每段序列第一行为＞Tag-i,第二行为所切基因组序列。例1:>Tag-lTATTGCAC.GAAG]CTCCCA>

2、Tag-2gtattataHatctcctcQccttac>Tag-3attatgac・agttactc(

3、ccttct1、题目写《生物信息学编程试题》,不要写成生物信息学作业，不可以《生物信息学综述》为题目。2、交作业时，写明专业、年级、姓名和学号。姓名请手写签名。3、禁止抄袭和代答，雷同卷低于60分，综述不高于75分。4、7月1日前提交、过期不候。（注：将文件exam.fa存放在D：中，运行“生物信息学1”、“生物信息学2”源程序（见附件），即可得到相应输出文件genoma.fa>stats.fa、Tag.fa。查看・fa文件，可用UltraEdit软件打开，软件安装包见附件）

4、编程1#include"st:ring.h"#include"stdio.h"#includc"stdlib.h"^defineMAXLEN90intmain(void)inti=0;intAcount=0,Tcount二0,Ccount=0,Gcount=0;charbuff[MAXLEN]={0};*F；*fp；F=freopen(，ZD:\exam.faVx，stdin);//文件指针*F指向打开的文件exam.⑴freopen(z/l):Wgcnome.faz，,stdout);//建立1)A序列文件genome,fd,并进行写操作,输LB缓冲区的内容写入到文件当+F二

5、二NULL)printf(z，can，tfineesam・fa〃)；while(!feof(E))fgets(buff,MAXLENT,F);i++；if((i%2)==0)printf(^sW,buff);//输出偶数行Acount++;elseif(num==，)Tcount++;elseif(num=='C)Ccount++;elseif(num==，G*)Geount++;}charnum=0;while(!feof(fp))〃若文件指针*5不为空，贝0进行读取字符操作■num=fgetc(fp);〃读取genome,fa中的字符并输出到字符变量num屮if(num=-

6、A*)//对读取的字符变晁进行判断Ifp=fopenC'D:Wgenome.faz/,"r");//读取genome,fa屮的数^streamstream二fopen("D:\slals・C,"w")

7、//文件指针*stream指向新建立的文件stats,fafprintf(stream,"A二%d,",Acount);//将统计的A、T、C、G的个数通过指针输入到文件ats.fa中fprintf(stream,z/,Tcount)fprintf(stream,fprintf(stream,"G=%d",Gcountfprintf(stream,"共％d〃」count+Tco

8、unt+Ccount+Gcount);//A、T、C、G个数的总共fclose(stream);//关闭文件指针return()二」nop」s-、、M、、"、、ej.胃匚

9、j一(oll丄7「启jnq密・v・==〔「〕JJnq)J)(++「・『><吏〉「o丄)Jo」二VNdH—N5nxvwyjnq)S4①<2?((vw)40£io二言printf(，z>Tag-%ld/z,++n);printf(，，%s，/,enzymecutting)}fclose(stdout);return0