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

修改可读取压缩格式文件

 

#include <stdio.h>
#include <stdlib.h>
#include <string.h>
#include <zlib.h>

#define RLEN 1024

static void getfq(char *fq);

int main(int argc, char *argv[]){

    if(argc!=2){
        fprintf(stderr,"\tusage: %s fastq[.gz]\n",argv[0]);
        exit(-1);
    }
    getfq(argv[1]);
    exit(0);
}

static void getfq(char *fq){

    gzFile fp1=gzopen(fq,"r"); 
    char fline[RLEN]; 
    char read[4][RLEN]; 
    char *r=NULL; 

    while(fp1 != NULL){    
        for(int i=0; i<4; i++){
            r=gzgets(fp1,fline,RLEN);
            if(r != NULL){
                strcpy(read[i],fline);
            }
            else{ // EOF
                goto LoopOut;
            }
        }
        for(int i=0; i<4; i++){
            if(i==0){
                memcpy(read[i],">",1);
                fprintf(stdout,"%s",read[i]);
            }
            else if(i==1){
                fprintf(stdout,"%s",read[i]);
            }
            else{
                continue;
            }
        }
    }
    LoopOut:
        free(r); 
        r=NULL;
        gzclose(fp1);

    return;
}