+
+ FILE *fasta;
+
+ if ((fasta=fopen(argv[1],"r"))==NULL) {
+ printf("Could not open %s\n",argv[1]),exit(1);
+ }
+
+/* Creating output files depending on the type */
+ createOutputFiles(type);
+
+/* Read input file sequence by sequence */
+
+ SEQS *fastaseq;
+
+ do {
+
+ fastaseq = gseq_fasta(fasta);
+ if(fastaseq==NULL) {
+ break;
+ }
+ SEQ=malloc(sizeof(SEQ_STR));
+
+ //printf("No: %s\n",fastaseq->id );
+ //printf("L: %d\n",fastaseq->slen);
+
+ SEQ->name=fastaseq->id;
+ SEQ->seq=fastaseq->seq;
+ SEQ->le=fastaseq->slen;
+
+#ifdef DEBUG
+ printf("N: %s\n",SEQ->name);
+ printf("S: %s\n",SEQ->seq);
+ printf("L: %d\n",SEQ->le);
+#endif
+
+