1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30
|
int get_trinucls(char* sequence, int position_start, int position_end, char* trinucl_sens, char* trinucl_anti, int cut_rank) {
int j;
char base_anti, base_sens;
for (j=0; j<3; j++) {
base_anti = sequence[position_end-cut_rank-j-1];
base_sens = sequence[position_start+j+(ENZYME_LEN-cut_rank)-1];
if( (base_anti == 'a' || base_anti == 'c' || base_anti == 'g' || base_anti == 't') && (base_sens == 'a' || base_sens == 'c' || base_sens == 'g' || base_sens == 't') ) {
trinucl_sens[j] = base_sens;
// if (j == 0) printf("%c", trinucl_sens[j]);
switch (base_anti) {
case 'a' : base_anti = 't'; break;
case 'c' : base_anti = 'g'; break;
case 'g' : base_anti = 'c'; break;
case 't' : base_anti = 'a'; break;
}
trinucl_anti[j] = base_anti;
}
else {
return 0;
}
}
trinucl_sens[3]='\0';
trinucl_anti[3]='\0';
// printf("\n");
// et là j'ai exactement les memes caractères !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
printf("%c\t%c\n", trinucl_anti[4], trinucl_sens[1]);
// printf("%s\t%u\t%s\t%u\t%s\n", psub->name, position_start, trinucl_sens, position_end, trinucl_anti);
return 1;
} |
Partager