28 lines
573 B
C
28 lines
573 B
C
|
#ifndef _GLOBAL_DEFS_H
|
||
|
|
#define _GLOBAL_DEFS_H
|
||
|
|
#include "global_defs.h"
|
||
|
|
#endif
|
||
|
|
|
||
|
|
main()
|
||
|
|
{
|
||
|
|
Sequence tSeq;
|
||
|
|
int ii, jj, rr;
|
||
|
|
char acgt[128];
|
||
|
|
|
||
|
|
for(ii = 0; ii < 128; ii++)
|
||
|
|
acgt[ii] = '0';
|
||
|
|
|
||
|
|
acgt[0x41] = acgt[0x43] = acgt[0x47] = acgt[0x54] = acgt[0x55] = '1';
|
||
|
|
acgt[0x61] = acgt[0x63] = acgt[0x67] = acgt[0x74] = acgt[0x75] = '1';
|
||
|
|
|
||
|
|
while((rr = ReadRecord(stdin, &tSeq)) != -1)
|
||
|
|
{
|
||
|
|
strcpy(tSeq.type, "MASK");
|
||
|
|
for(ii = 0; ii < tSeq.seqlen; ii++)
|
||
|
|
{
|
||
|
|
tSeq.c_elem[ii] = acgt[tSeq.c_elem[ii]];
|
||
|
|
}
|
||
|
|
WriteRecord(stdout, tSeq, NULL, 0);
|
||
|
|
}
|
||
|
|
}
|