C语言 fastq文件转换为fasta文件

目前只能处理短序列，若要处理长序列，可按照https://www.cnblogs.com/mmtinfo/p/13036039.html的读取方法。

 1 #include <stdio.h>
 2 #include <stdlib.h>
 3 #include <string.h>
 4 #define RLEN 1024
 5 
 6 static void getfq(char *fq, char *ofq);
 7 
 8 int main(int argc, char *argv[]){
 9 
10     if(argc!=3){
11         printf("infomation:\n");
12         fprintf(stderr,"\tusage: %s <fastq> <fast>\n",argv[0]);
13         exit(-1);
14     }
15     getfq(argv[1],argv[2]);
16     exit(0);
17 }
18 
19 static void getfq(char *fq, char* ofq){
20 
21     FILE *fp1=fopen(fq,"r"); 
22     FILE *fp2=fopen(ofq,"w"); 
23 
24     char fline[RLEN]; 
25     char read[4][RLEN]; 
26     char *r=NULL; 
27 
28     while(fp1 != NULL && fp2 != NULL){
29         
30         for(int i=0; i<4; i++){
31             r=fgets(fline,RLEN,fp1);
32             if(r != NULL){
33                 strcpy(read[i],fline);
34             }
35             else{ // EOF
36                 goto LoopOut;
37             }
38         }
39 
40         for(int i=0; i<4; i++){
41             switch(i)
42             {
43                 case 0:
44                     memcpy(read[i],">",1);
45                     fprintf(fp2,"%s",read[i]);
46                     break;
47                 case 1:
48                     fprintf(fp2,"%s",read[i]);
49                     break;
50                 case 2:
51                     continue;
52                     break;
53                 case 3:
54                     continue;
55                     break;
56                 default:
57                     break;
58             }
59         }
60     }
61     LoopOut:
62         free(r); 
63         r=NULL;
64         free(fp1);
65         free(fp2);
66 
67     return;
68 }

 

posted @ 2020-11-04 10:34  天使不设防  阅读(303)  评论(0)  编辑  收藏  举报