Accept to deal with sequence in lower case
git-svn-id: https://www.grenoble.prabi.fr/svn/LECASofts/ecoPrimers/trunk@214 60f365c0-8329-0410-b2a4-ec073aeeaa1d
This commit is contained in:
@ -100,6 +100,8 @@ ecoseq_t *readnext_ecoseq(FILE *f)
|
|||||||
int32_t comp_status;
|
int32_t comp_status;
|
||||||
unsigned long int seqlength;
|
unsigned long int seqlength;
|
||||||
int32_t rs;
|
int32_t rs;
|
||||||
|
char *c;
|
||||||
|
int32_t i;
|
||||||
|
|
||||||
raw = read_ecorecord(f,&rs);
|
raw = read_ecorecord(f,&rs);
|
||||||
|
|
||||||
@ -144,6 +146,9 @@ ecoseq_t *readnext_ecoseq(FILE *f)
|
|||||||
if (comp_status != Z_OK)
|
if (comp_status != Z_OK)
|
||||||
ECOERROR(ECO_IO_ERROR,"I cannot uncompress sequence data");
|
ECOERROR(ECO_IO_ERROR,"I cannot uncompress sequence data");
|
||||||
|
|
||||||
|
for (c=seq->SQ,i=0;i<seqlength;c++,i++)
|
||||||
|
*c=toupper(*c);
|
||||||
|
|
||||||
// fprintf(stderr,"seq name : %30s seq size : %d\n",seq->DE,seq->SQ_length);
|
// fprintf(stderr,"seq name : %30s seq size : %d\n",seq->DE,seq->SQ_length);
|
||||||
return seq;
|
return seq;
|
||||||
}
|
}
|
||||||
|
Reference in New Issue
Block a user