7 static void replace_unusual( int n, char **seq, char *usual, char unknown, int (*uporlow)( int ) )
16 if( !strchr( usual, *pt ) ) *pt = unknown;
17 else *pt = uporlow( *pt );
24 void arguments( int argc, char *argv[] )
31 while( --argc > 0 && (*++argv)[0] == '-' )
33 while ( (c = *++argv[0]) )
38 seedoffset = atoi( *++argv );
39 fprintf( stderr, "seedoffset = %d\n", seedoffset );
44 fprintf( stderr, "inputfile = %s\n", inputfile );
48 fprintf( stderr, "illegal option %c\n", c );
58 fprintf( stderr, "options: Check source file !\n" );
65 int main( int argc, char *argv[] )
75 arguments( argc, argv );
79 infp = fopen( inputfile, "r" );
82 fprintf( stderr, "Cannot open %s\n", inputfile );
91 getnumlen_casepreserve( infp, &nlenmin );
93 fprintf( stderr, "%d x %d - %d %c\n", njob, nlenmax, nlenmin, dorp );
95 seq = AllocateCharMtx( njob, nlenmax+1 );
96 name = AllocateCharMtx( njob, B+1 );
97 nlen = AllocateIntVec( njob );
99 readData_pointer_casepreserve( infp, name, nlen, seq );
101 // for( i=0; i<njob; i++ ) gappick_samestring( seq[i] );
105 origfp = fopen( "_original", "w" );
108 fprintf( stderr, "Cannot open _original\n" );
111 for( i=0; i<njob; i++ )
113 nlen[i] = strlen( seq[i] );
114 fprintf( origfp, ">%s\n", name[i]+1 );
115 if( seq[i][nlen[i]-1] == '\n' ) seq[i][nlen[i]-1] = 0;
116 fprintf( origfp, "%s\n", seq[i] );
123 usual = "ARNDCQEGHILKMFPSTWYVarndcqeghilkmfpstwyv-.";
124 replace_unusual( njob, seq, usual, 'X', toupper );
128 usual = "ATGCUatgcu-";
129 replace_unusual( njob, seq, usual, 'n', tolower );
134 for( i=0; i<njob; i++ )
136 fprintf( stdout, ">%d_oripos_%s\n", i+seedoffset, name[i] );
137 fprintf( stdout, "%s\n", seq[i] );